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PREFACE 


This  document  Is  volume  one  (D2- 125286-1)  of  three 
volumes.  It  describes  a  FORTRAN  IV  digital  computer 
program  developed  for  analysis  of  solid  rocket  motor 
Internal  ballistics.  Volume  one,  "Program  Manual," 
explains  the  theory  and  describes  the  mathematical 
model,  program  capabilities  and  Information  necessary 
for  program  maintenance  and  revision.  Volume  two, 
"User’s  Guide,"  (02-125286-2)  describes  program  op¬ 
tions,  preparation  of  program  Input  data  and  program 
output.  Volume  three,  "Sample  Case  Results  and  Program 
Listings,"  (02-125286-3)  contains  the  sample  case  re¬ 
sults  and  complete  program  listings  of  the  computer 
program.  The  document  Is  divided  Into  three  volumes 
for  handling  convenience.  Section  numbering  Is  con¬ 
tinuous  through  the  three  volumes.  A  complete  table 
of  contents  appears  In  each  volume. 
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ABSTRACT 


This  report  describes  a  FORTRAN  IV  digital  computer  pro- 
gram  developed  to  calculate  internal  ballistic  performance 
of  solid  propellant  rocket  motors  with  high  burning  rates, 
short  burning  durations,  and  high  vehicle  accelerations. 
Forked  wagon  wheel,  conventional  wagon  wheel,  standard 
star,  slot ted -cone,  and  circular  port  monolithic  and 
segmented  grain  designs  may  be  considered.  Accurate 
description  of  an  Inert  sliver  in  the  cylindrical  section 
Is  allowed  for  all  but  the  forked  wagon  wheel  grain  design. 
The  effects  of  anlsotroolc  burning  of  the  propellant  may 
be  considered.  The  storage  of  mass  and  momentum  (capac- 
Itance  effects)  and  vehicle  acceleration  are  Included  In 
the  internal  gas  dynamic  equations.  Ignition  transients 
may  be  calculated.  Tabular  Input  of  the  motor  grain 
description  is  available  for  special  grain  configurations 
that  cannot  be  described  by  the  program  geometry  constants. 
Motor  performance  parameters  such  as  delivered  and  vacuum 
thrust  and  total  Impulse,  fore-head  and  sf'"head  total 
pressure,  nozzle  discharge  flow,  fore-heat:  <nd  aft-head 
total  pressure  Integrals,  pitch  and  roll  moments  of 
Inertia,  center  of  gravity  locations,  burn  surface  area, 
and  weight  of  propellant  remaining  are  printed  for  each 
time  Interval. 

This  report  Is  divided  Into  three  volumes.  Volume  one  Is 
the  Program  Manual,  volume  two  Is  the  User1*  Guide,  and 
volume  three  Is  the  Sample  Case  Results  and  Program  Listings. 


KEY  WORDS 


The  following  Key  Words  Identify  the  major  program  capa¬ 
bilities; 

$ 

Internal  Ballistics 

Solid  Propellant  Rocket  Mo*or 

Monolithic  Grain 

Segmented  Grain 

One-D I mens  local  Gas  Dynamics 

Steady  Flow  Gas  Dynamics 

Non-Steady  Flow  Gas  Dynamics 

Isotropic  Propellant  Burning 

Anisotropic  Propellant  Burning 

Vehicle  Acceleration 

ignition  Transient  Interval 

Web -Time  Interval 

Tall -Off  Interval 

Fore-Head  Section 

Aft-Head  Section 

Cylindrical  Section 

Center  of  Gravity 

Moment  of  Inertia 

Grain  Geometry 

Nozzle 


ACKNOWLEDGEMENTS 


Acknowledgement  Is  given  to  the  Thlokol  Chemical  Corpora¬ 
tion,  Redstone  Division,  for  the  original  SCAT  machine 
language  program  from  which  this  program  was  developed 
(References  1  and  2),  and  to  the  Thlokol  Chemical  Corpora¬ 
tion,  Wasatch  Division,  for  the  documentation  of  the  grain 
geometry  calculations  from  which  the  description  of  the 
mathematical  model  was  developed  (Reference  3) « 


DEFINITION  OF  TERMS 


The  following  general  definitions  apply  throughout  this 
document: 

Internal  Ballistics:  Analysis  of  the  burning  characteristics 
and  progression  of  the  propellant  surface,  dynamics  of  the  gas 
flow,  and  the  gas  generation  In  the  Interior  of  solid  propellant 
rocket  motors. 

Performance  Characteristics:  Parameters  that  specify  motcr 
performance,  e.g.,  thrust  vs  time,  maximum  chamber  pressure, 
specific  Impulse,  and  burn  time. 

Gas  Dynamics:  Study  of  the  generation  and  flow  of  combustion 
products  along  the  propellant  grain  and  through  the  nozzle. 

Steady  Flow  Gas  Dynamics:  Mass,  energy,  and  momentum  within 
a  control  volume  are  constant  with  time. 

Non-Steady  Flow  Gas  Dynamics;  Mass,  energy,  and  momentum  within 
a  control  volume  are  not  constant  with  time. 

Grain  Design:  The  cross  sectional  grain  configuration  of  the 
propellant. 

Monolithic  Grain?  The  propellant  grain  Is  one  single  piece. 

Segmented  Grain:  The  propellant  grair^  is  divided  up  Into  a 
number  of  longltudlial  segments. 

Slots:  The  region  between  the  segments  which  does  not  contain 
propellant. 

Web:  The  minimum  distance  between  the  grain  surface  and  the 
case  wall. 

Core:  The  region  occupied  by  the  combustion  gases. 

Reference  Planes:  The  stations  In  the  cylindrical  section  of 
a  motor  where  the  grain  design  Is  specified. 

Increment  Dividing  Planes:  The  stations  tr.  the  cylindrical 
section  of  the  motor  where  the  solution  of  the  Internal  gas 
dynamic  equations  is  obtained. 

Mass  Addition  Region:  The  region  between  Increment  dividing 
planes. 


Ignition  Transient  Interval;  The  time  required  to  obtain  motor 
operating  pressure. 

Web-Tine  Interval;  The  time  required  to  burn  through  the  web. 

Tail-off  Interval?  The  time  Interval  after  web  burn  through. 

Isotropic  Propellant  Burning;  Where  the  burning  rate  character¬ 
istics  are  Independent  of  distance  burned. 

Anisotropic  Propellant  Burning:  Where  the  burning  rate  character 
Istlcs  are  dependent  on  distance  burned. 
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INTRODUCTION 


This  report  describes  a  FORTRAN  IV  digital  computer  program, 
developed  by  The  Boeing  Company  for  operation  on  an  IBM  7094 
computer  with  an  IBSYS  version  13  monitor  system  from  an  earlier 
Thlokol  Chemical  Corporation  program,  to  perform  dynamic  analyses 
of  solid  propellant  rocket  motor  Internal  ballistics.  A  complete 
description  of  the  mathematical  model,  the  engineering  equation 
development,  the  method  of  solution,  and  detailed  programming 
information  is  presented  to  xplaln  program  principles  and  theory 
and  to  facilitate  program  maintenance  and  revision. 

This  program  was  developed  under  contract  from  the  Army  Missile 
Command  Propulsion  Laboratory,  Redstone  Arsenal,  Alabama, 
contract  number  DA-01-021-ANC-15557(Z)  by  the  Engineering 
Digital  Computing  Organization  2-2640,  with  support  from  the 
Klsstle  and  Information  Systems  Division,  Flight  Technology 
Propulsion  Organization  2-5711  of  The  Boeing  Company. 

Generc!  Information 

The  report  Is  divided  Into  three  volumes.  Volume  1,  the  Program 
Manual .  provides  a  technical  explanation  of  the  theory,  mathe¬ 
matical  model,  program  capabilities  and  Information  necessary 
for  program  maintenance  and  revision.  Volume  II,  the  User's 
Guide,  describes  program  options,  preparation  of  Input  data  and 
program  output.  Volume  III,  the  Sample  Case  Results  and  Program 
Listings ,  contains  the  complete  program  listings  of  the  computer 
program.  The  following  paragraphs  deserf briefly  the  sections 
of  the  report  comprising  each  volume. 

Volume  I 


The  Program  Capabilities  and  Limitations,  Section  2.C,  Indicates 
the  capability  the  program  has  to  evaluate  grain  designs  and 
Internal  ballistics  and  the  program  limitations  that  ex.st  In 
these  areas. 

The  Method  of  Solution,  Section  3*0;  describes  the  method  the 
program  uses  to  obtain  the  Internal  ballistic  solutions  and  the 
organization  of  major  program  sections  which  divide  the  solution 
Into  logical  blocks  or  core  loads  that  reside  tn  core  at  separate 
times . 

The  Gas  Dynamic  Equation  Development,  Sect  ton  4.0,  presents  the 
development  of  the  equations  for  the  non-steady  flow  gas  dynamics 
for  both  segmented  and  monolithic  motors,  propellant  description, 
end  for  an  accelerating  reference  system.  In  general,  this 
section  describes  modifications  made  to  the  original  Thlokol 
Chemical  Corporation  program  for  conducting  design  studies  of 
solid  propellant  configurations  (Reference  1)  to  simulate  the 


(Continued) 


Internal  ballistics  of  high  burning  rate  propellants  with  charac¬ 
teristically  short  burn  durations.  These  mod  If l cat  Ions  were 
made  specifically  to  include  the  storage  of  mass  and  momentum 
In  the  gas  dynamic  equations,  to  consider  the  effects  of  very 
high  vehicle  acceleration  on  the  Internal  ballistics,  and  to 
study  anisotropic  propellant  burning.  As  a  result  of  these 
modifications,  ignition  transients  may  be  calculated. 

The  Geometry  Equation  Development,  Section  5.0,  explains  in 
detail  the  setup  and  solution  of  the  grain  geometry  equations 
which  determine  the  perimeter  length,  cross  sectional  area,  burn 
surface  area,  moments  of  Inertia,  and  center  of  gravity  location 
of  the  various  grain  options  and  longitudinal  configurations. 

The  equation  development  and  figures  presented  In  this  section 
were  obtained  from  References  1  and  3* 

The  Detailed  Programming  Information,  Section  6.0,  presents  a 
brief  description  of  all  program  subroutines,  macroscopic  program 
logic  flow  charts,  a  description  of  the  computing  system  and 
program  storage  allocation,  program  diagnostic  aids,  and  a  list 
of  the  program  nomenclature.  Appropriate  comments  are  placed 
throughout  the  program  listings  as  a  supplement  to  the  subroutine 
flow  charts  to  assist  In  program  maintenance  and  revision. 

The  Results,  Section  7.0,  presents  a  comparison  of  a  computer 
prediction  with  three  full  scale  HI BEX  motor  firings.  Dimen¬ 
sionless  fore-head  pressure  traces  are  shown  for  the  computer 
prediction  and  the  motor  firings. 

Volume  ! I 


The  User's  Guide,  Section  8.0,  presents  an  explanation  of  the 
required  program  inputs,  sample  cases  showing  the  available 
program  options,  a  description  of  the  output  format,  and  the 
required  control  cards  to  permit  effective  program  use  and 
operation  without  knowledge  of  the  program  technical  aspects. 

This  section  Is  arranged  to  be  complete  without  reference  to 
the  program  manual  technical  sections  and  may  be  used  Independent 
of  the  program  manual. 

Volume  i 1 1 


The  Listings,  Section  9.0,  contains  the  sample  case  results  and 
program  lti  .  ings  of  the  computer  program. 


HUtory  of  Program  Development 

Jn  March  1360,  work  was  Initiated  at  the  Thtokol  Chemical 
Corporation,  Redstone  Division,  under  the  auspices  of  Systems 
Analysis  Laboratory,  Army  Rocket  and  Guided  Missile  Agency, 
Redstone  Arsenal,  Alabama,  for  the  development  of  a  solid 
propellant  rocket  motor  design  program  (References  1  and  2). 

In  1962,  The  Boeing  Company  received  a  copy  of  this  program, 
and  In  1963  developed  a  segmented  motor  version.  In  1964, 
the  Saturn  Branch  of  the  Aerospace  Division  of  Boeing  at  New 
Orleans  converted  the  SCAT  machine  language  program  to  FORTRAN 
M  for  operation  on  the  IBM  7094  computer.  In  1365,  The 
Boeing  Aerospace  Division  In  Seattle,  Washington,  converted 
the  FORTRAN  li  version  to  FORTRAN  IV  for  operation  on  the  SRU 
1107  and  made  the  modifications  discussed  In  Section  1.1  for 
the  HtBEX  contract  to  add  the  transient  capabtittyo  In  May 
1966,  The  Boeing  Company  proposed  to  the  Army  Missile  Command 
to  segment  the  existing  SRU  1108  program,  perform  the  necessary 
conversion  for  operation  on  the  IBM  7094,  and  completely 
document  the  advanced  program  version.  A  contract  was  received 
In  July  1966  from  the  Army  Missile  Command  for  a  6  month 
development  effort  to  perform  the  required  conversion. 
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PROGRAM  CAPABILITIES  AND  LIMITATIONS 


This  computer  program  was  developed  to  calculate  solid  propellant 
rocket  motor  Internal  ballistics.  Because  of  Its  development 
from  earlier  grain  design  and  ballistic  performance  programs, 
additional  capabilities  are  preser.t.  Throughout  the  development 
effort,  all  prior  program  capabilities  have  been  retained  so 
that  a  general  program  exists  with  both  grain  design  and  internal 
ballistic  evaluation  capability. 

2.1  Program  Capabl lit les 

The  basic  propellant  grain  design  Is  the  forked  wagon  wheel, 
however,  the  grain  design  equations  are  general  so  that  the  con¬ 
ventional  wagon  wheel,  standard  star  and  circular  port  as  well 
as  the  slotted-cone  may  be  described  by  variations  In  the  Input 
data,  Figure  2.1  shows  the  five  grain  design  options.  Other 
more  complicated  grain  designs  may  be  evaluated  by  describing 
the  perimeter  length  and  burn  area  as  a  function  of  distance 
burned  and  Input  to  the  computer  program  as  tabular  data. 

The  propellant  grain  configuration  may  be  either  monolithic  or 
segmented  with  up  to  11  slots.  The  propellant  case  and  port 
cavity  may  be  either  cylindrical  or  tapered.  The  fore-head 
section  configuration  mav  be  a  straight  through  grain  or  may 
contain  a  complete  web.  Figure  2.2  shews  the  various  motor 
configuration  options.  The  aft-head  section  configuration  Is 
a  straight  through  grain. 

The  propellant  characteristics  are  described  by  definitive 
properties  of  the  combustion  gases  and  a  generalized  burning 
rate  equation.  The  propellant  gas  properties  may  be  held 
constant  or  may  be  varied  as  a  function  of  the  static  pressure 
In  the  port  cavity.  The  burning  rate  model  Includes  erosive 
burning  and  will  allow  either  Isotropic,  or  anisotropic  burning 
of  the  propellant  surface. 

Either  steady  or  non-steady  flow  gas  dynamics  are  available  to 
obtain  the  Internal  ballistic  solution.  The  steady  flow  gas 
dynamics  solve  the  momentum  and  continuity  equations  without 
consideration  of  time  dependent  terms  such  that  there  Is  no 
storage  of  mass  or  momentum  (no  capacitance  effects).  The  non¬ 
steady  flew  gas  dynamics  solve  the  momentum  and  continuity 
equations  with  respect  to  time  so  that  the  storage  of  mass  and 
momentum  Is  considered  and  the  start  transient  and  tali -off 
Intervals  may  be  determined. 

The  effects  of  vehicle  acceleration  on  the  Internal  ballistic 
solution  may  be  considered.  The  acceleration  term  Is  included 
in  the  momentum  equation. 

The  effect  of  a  tapered  inert  silver  may  be  considered  In  the 
cylindrical  section  for  all  but  the  forked  wagon  wheel  grain 
design. 
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Program  Capabl lltlec  (Continued) 
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The  following  Is  a  summary  of  the  program  capabilities; 


1.0  Grain  Design 

A.  Circular  Port 

B.  Standard  Star 

C.  Slotted-Cone 

D.  Conventional  Wagon  Wheel 

E.  Forked  Wagon  Wheel 

F.  Tables  of  Perimeter  and  Burn  Area  as  a  Function 
of  Distance  Burned  can  be  Input 

2.0  Kotor  Configuration 

A.  Cylindrical  Section 

1)  Monolithic  grain 

2)  Segmented  grain 

3)  Tapered  Inert  silver 

B.  Fore-head  Section 

1)  Straight  through  grain 

2)  Complete  web  (head -end  with  web) 

C-  Aft-head  Section 

1)  Straight  through  grain 

2)  Burning  on  aft  face 

3.0  Propellant  Characteristics 

A.  isotropic  Propellant  Burning 

B.  Anisotropic  Propellant  Burning 

C.  Erosive  Burning 

D.  Variable  Gas  Properties  (function  of  static  pressure) 
4.0  Internal  Ballistics 

A.  Steady  Gas  Flow 

B.  Non-Steady  Gas  Flow 

1)  Ignition  transient  Interval 

2)  Web-time  Interval 

3)  Tall -off  Interval 

C.  Vehicle  Acceleration 
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Program  Output 

the  program  provides  the  following  output: 

1.0  Motor 

A.  Delivered  Thrust 

B.  Vacuum  Thrust 

C.  Fore-head  Total  Pressure 

D.  Nozzle  Total  Pressure 

E*  Fore-head  Total  Pressure  Time  Integral 

F.  Nozzle  Total  Pressure  Time  Integral 

G.  Delivered  Total  Impulse 

H.  Vacuum  Total  Impulse 

I.  Nozzle  Discharge  Flow  Rate 

J.  Polar  and  Rectangular  Moment  of  inertia 

K.  Center  of  Gravity 

2.0  Propellant  Characteristics 

A.  Weight  of  Propellant  Remaining 

3.  Forward  Tangent  Plane  Propellant  Burning  Rate 

C.  Aft  Tangent  Plane  Propel  lar'-  Burning  Rate 

D.  Total  '  > 1 gh t  of  Propellant  Expended 

3.0  Grain  Geometry 

A.  Cylindrical  Section  Burn  Area 
8.  Fore-head  Section  Burn  Area 

C.  Aft-head  Section  Burn  Are* 

D.  Grain  Segment  Face  Burn  Area 

E.  Total  Motor  Burn  Area 

4.0  Nozzle  Characteristics 

A.  Throat  Area 

B.  Effective  Expansion  Ratio  (flow  separation  accounted  for) 
C«  Pressure  Ratio  Across  Nozzle 

D.  Momentum  Portion  of  Thrust  Coefficient 

Program  Assumptions 

The  following  assumptions  were  made  to  translate  the  physical 
system  Into  the  one-dimers  Iona  1  gas  dynamic  model: 

1.  Propellant  burning  during  Ignition  and  steady  state 
operation  occurs  normal  to  the  grain  surface. 

2.  The  burn  rate  In  the  fore-head  and  aft-head  sections 
Is  assumed  to  be  constant  over  the  burning  surface 
of  the  entire  section. 


Program  Assumptions  (Continued) 


3.  Hass  addition  occurs  Instantaneous^  with  no  velocity 
component  along  the  motor's  longitudinal  axis  (dZ/dt  *  0). 

4.  The  products  of  combustion  obey  the  perfect  gas  law. 

5.  The  gas  flow  Is  one-df:  ansional  and  adiabatic. 

6.  The  combustion  temperature  Is  constant  throughout  the 
motor. 

7.  The  heat  capacity  of  the  combustion  gases  Is  constant. 

8.  The  friction  forces  of  the  combustion  gases  In  the  port 
cavity  a-e  negligible. 

9.  The  moments  oi  Inertia  about  the  pitch  and  yaw  axis  are 
equa 1 . 

Limitations 

The  program  has  the  following  limitations! 

1.  A  maximum  of  11  reference  planes  are  allowed  In  the 
cylindrical  section  to  describe  the  grain  design. 

2.  A  maximum  of  100  Increment  dividing  planes  are  allowed 
In  the  cylindrical  section  to  define  the  mass  addition 
regions.  If  this  restriction  Is  exceeded  by  defining 

a  AZ  too  small,  the  case  execution  will  be  terminated 
and  an  appropriate  comment  will  be  printed. 

3.  A  maximum  of  11  slots  are  allowed  for  segmented  motors. 

4.  The  slotted-cone  grain  design  Is  applicable  only  to  the 
cylindrical  section.  Burn  area  tables  must  be  Input 
for  the  forward  and  aft  domes  when  this  grain  design  Is 
used. 

5.  The  Inert  silver  option  Is  restricted  to  the  cylindrical  { 
section  and  does  not  apply  to  either  end  section. 

6.  The  effects  of  an  accelerating  reference  system  can  be 
determined  only  for  non-steady  flow  gas  dynamics. 
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Forked  Wagon  Wheel 


Conventional  Wagon  Wheel 


Standard  Star 


Slotted-Cone 


Circular  Port 


Figure  2*1.  Grain  Design  Options 


ntinuy  ur  sulunun 


The  program  calculations  are  based  on  the  geometrical  model 
shown  In  Figure  3*1*  The  motor  configuration  Is  divided  Into 
three  sections:  the  head-end  section  (forward  dome),  the 
cylindrical  section,  and  the  aft-head  section  (aft  dome)  or 
nozzle.  The  grain  geometry  Is  described  by  Input  reference 
planes  within  the  cylindrical  section.  The  cylindrical  section, 
which  may  contain  either  a  monolithic  or  a  segmented  grain.  Is 
further  divided  Into  a  number  of  Increments  or  mass  addition 
regions  by  the  location  of  Increment  dividing  planes,  at  each 
reference  plane  and  at  specified  Intervals  (AZ)  from  each  re¬ 
ference  plane  until  either  a  segment  slot  Interface  or  the  next 
reference  plane  Is  passed.  During  the  computer  solution  of  the 
gas  dynamics,  port  perimeter,  port  cross  sectional  area  and  mo¬ 
ments  of  Inertia  are  determined  at  each  Increment  dividing  plane 
b/  llnec~  Interpolation  between  adjacent  reference  planes.  Mass 
addition  Is  assumed  to  occur  as  a  step  process  between  two  In¬ 
crement  dividing  planes. 

The  program  method  of  solution  Is  divided  into  four  separate 
computer  core  loads  linked  together  by  a  main  control  program 
and  a  common  data  region.  Each  core  load  Is  unique,  but  depen¬ 
dent  on  preceding  core  loads  for  generated  data.  The  flow 
chart  shown  In  Figure  3*2  presents  the  macroscopic  program 
order  of  solution  with  the  separate  core  loads  linked  together. 

The  first  computer  core  load  contains  the  subroutines  required 
to  read  the  Input  data,  initialize  the  data  cells,  compute  the 
Input  reference  plane  constants,  locate  the  Increment  dividing 
planes,  check  for  Input  data  errors,  and  print  the  program  Inputs 
and  computed  constants. 

The  second  computer  core  load  contains  the  geometry  subroutines 
required  to  compute  the  initial  propellant  cross  section  area 
and  perimeter  length  for  the  cylindrical  section  reference  planes, 
thn  aft-head  and  straight  through  grain  fore-head  sections  burn 
area  and  Initial  propellant  volume  as  a  function  of  distance 
burned.  The  moments  of  inertia  and  centers  of  gravity  for  the 
a't-head  and  straight  through  grain  fore-head  sections  and  the 
radius  of  gyration  for  the  cylindrical  section  are  also  calculated. 
The  computed  yalues  for  each  section  are  stored  In  tables  for  use 
during  the  solution  of  the  internal  ballistics  in  the  fourth  core 
load. 

The  third  computer  core  load  contains  the  geometry  subroutines 
required  to  compute  the  Initial  propellant  volume,  burn  area, 
moments  of  inertia,  and  CG  location  tables  of  the  head-end  with 
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web.  The  third  core  load  1$  loaded  only  If  a  head-end  with  web 
Is  required  after  the  cylindrical  section  and  aft-head  section 
geometry  calculations  have  been  completed  In  the  second  core  load. 

After  the  grain  geometry  calculations  a*e  reformed  In  the  first, 
second,  and  third  core  loads  and  the  perimeter  and  area  tables 
have  been  established,  the  Internal  ballistic  solution  Is  Initiated 
In  the  fourth  core  load  by  determining  the  geometry  values  of 
each  reference  plane  and  each  end  section  ; rom  a  table  look-up 
procedure  In  the  geometry  tables.  For  steady  flow  conditions,  an 
Initial  estimate  of  the  fore-head  pressure  Is  made,  and  the  burning 
rate  (which  Is  assumed  to  be  constant  over  the  entire  head-end 
section)  Is  determined  as  described  In  Sect  Ion 4JLL  With  this 
burn  rate  and  the  tabular  value  of  the  fore-head  burn  area,  the 
Instantaneous  value  of  mass  addition  Is  determined  for  the  fore¬ 
head  section.  The  state  and  gas  dynamic  properties  of  the 
propellant  at  the  forward  tangent  plane,  the  first  Increment 
dlvld.ng  plane,  are  determined  from  the  simultaneous  solution 
of  the  momentum  and  continuity  equations  assuming  perfect  ga:; 
relationships < 

The  grain  geometry  at  the  first  Increment  dividing  plane  Is  then 
determined  and  stored  In  temporary  locations  for  future  reference. 
The  grain  geometry  of  Increment  dividing  plane  two  Is  determined. 
The  Increment  section  mass  generation  rate  Is  determined  from 
the  perimeter  lengths  of  Increment  dividing  planes  one  and 
two,  the  Increment  length,  and  the  burning  rates  at  the  upstream 
adjacent  Increment  dividing  planes.  The  propellant  properties  and 
mass  flow  at  Increment  dividing  plane  two  are  then  determined  by  a 
simultaneous  solution  of  the  momentum  and  continuity  equations 
for  either  steady  or  non-steady  flow  conditions. 

The  above  procedure  Is  repeated  for  each  cylindrical  section  mass 
addition  region.  Once  the  cylindrical  section  Is  complete,  the 
same  general  calculations  are  performed  for  the  aft-head  section. 
The  port  cross  sectional  area  and  perimeter  length  and  burning 
rate  are  assumed  to  be  constant  In  the  aft-head  section  and 
Identical  to  the  values  at  the  aft  tangent  plane. 

The  nozzle  throat  area  Is  compared  to  the  maximum  value  that  will 
maintain  subsonic  flow  In  the  port.  If  this  maximum  value  Is 
exceeded,  the  program  prints  an  error  comment  and  terminates  the 
case.  If  this  maximum  value  Is  not  exceeded,  the  flow  rate  of 
propellant  discharge  through  the  nozzle  Is  computed  on  the  bcsls 
of  Isentropic  flow.  This  flow  rate  Is  compared  to  the  flow  rate 
of  propellant  discharged  from  the  grain.  If  these  two  values  do 
not  agree  within  .1  percent,  the  fore-head  pressure  Is  adjusted 
and  the  program  returns  to  the  fore-head  and  repeats  until 
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convergence  Is  attained.  Once  equilibrium  Is  reached,  additional 
ballistic  properties  are  computed  and  the  performance  data  Is 
printed. 

Following  the  performance  printout,  the  thickness  burned  In  each 
Increment  dividing  plane  and  slot  interface  Is  then  determined 
from  the  previous  web  thickness.  A  check  Is  then  made  to  see 
If  burnout  has  occurred  at  any  of  the  Increment  dividing  planes. 

If  burnout  occurs,  a  comment  Is  printed  that  the  Increment 
dividing  plane  has  burned  out.  The  progression  of  the  slot  Inter¬ 
faces  for  segmented  motors  Is  Indicated  by  a  printout  of  the  In¬ 
crement  dividing  plane  longitudinal  location. 

The  time  Is  then  Incremented  by  the  computed  time  Interval  and 
the  program  returns  to  the  fore-head  to  compute  new  equilibrium 
conditions  and  determine  new  values  cf  the  perimeter  length  and 
port  cross  section  area  at  each  Increment  dividing  plane  and  burn 
area  for  each  end  section.  This  process  Is  then  repeated  until  the 
termination  option  is  exceeded. 

The  general  program  solution  of  the  Internal  ballistics  outlined 
above  Is  modified  for  non-steady  flow  conditions.  When  the  start 
transient  Interval  Is  computed,  the  fore-head  pressure  Is  defined 
by  tabular  Input  of  the  fore-head  pressure  as  a  function  of  time, 
or  the  burn  rate  coefficient  is  defined  by  tabular  Input  of  the 
burn  rate  coefficient  as  a  function  of  distance  burned.  When  the 
fore-head  pressure  Is  Input  for  the  start  transient,  the  burn 
rate  coefficient  Is  varied  to  obtain  convergence;  and  when  the 
burn  rate  coefficient  Is  Input  for  the  start  transient,  the  fore-  j 
head  pressure  Is  varied,  as  for  steady  flow  conditions,  to  obtain 
convergence.  When  the  fore-head  pressure  Is  Input  for  the  start 
transient,  an  Initial  estimate  of  the  burn  rate  coefficient  Is 
made  by  computing  a  first  guess  of  the  burn  rate  coefficient  from 
the  motor  configuration  parameters  and  the  fore-head  pressure 
variation.  With  this  burn  rate  coefficient,  the  instantaneous 
value  of  mass  addition  end  mass  discharge  Is  determined  for  the 
head-end  section.  The  propellant  gss  properties  for  the  first 
Increment  dividing  plane  are  then  determined  from  a  simultaneous 
solution  of  the  non-steady  gas  flow  equations  as  above.  The 
remainder  of  the  ballistic  solution  Is  unchanged. 


GAS  OYNAMIC  EQUATION  DEVELOPMENT 


The  static  and  total  pressure,  temperature,  velocity,  and  flow 
rate  of  the  gas  along  the  length  of  the  propellant  grain  and 
through  the  nozzle  are  required  to  determine  the  progression  of 
the  propellant  burning  surface  and  motor  performance  parameters 
such  as  chamber  pressure,  thrust,  and  total  Impulse.  These  para¬ 
meters  are  obtained  from  an  Iterative  solution  of  the  perfect 
gas  law  and  the  equations  of  continuity,  momentum,  and  energy  for 
one-d f mens lor.a  1  gas  flow.  The  gas  flow  along  the  propellant 
grain  Is  determined  by  dividing  the  grain  Into  a  number  of  Incre¬ 
ments  which  are  termed  mass  addition  regions.  Two  control  volumes 
are  shown  In  Figures  4.1  and  4.2  which  define  the  mass  addition 
regions  for  the  monolithic  and  segmented  motor  configurations 
shown  In  Figure  4.3.  The  gas  dynamic  equations,  which  are  solved 
for  each  mass  addition  region  along  the  propellant  grain,  are 
described  In  the  following  sections.  Thci  Isotropic  and  erosive 
burning  rate  equation,  anisotropic  mass  generation  model,  and 
propellant  gas  properties  required  for  the  calculations  are  also 
described. 

Gas  Dynamics  for  Incremental  Control  Volumes 

The  solution  of  the  continuity,  momentum,  and  energy  equations 
for  the  mass  addition  regions  along  the  propellant  grain  and 
through  the  nozzle  are  described  In  this  section.  Steady  flow, 
non-steady  flow,  or  non*steady  flow  with  acceleration  may  be 
selected. 

The  following  assumptions  are  basic  to  the  development  of  the  gas 
dynamic  equations: 

1.  Mass  addition  occurs  as  an  Instantaneous  process  with  no 
velocity  component  parallel  to  the  motor  axis  (dZ/dt  =  0). 

2.  The  products  of  combustion  obey  the  perfect  gas  law. 

3.  The  gas  flow  is  one-dimensional  and  adiabatic. 

4.  The  combustion  temperature,  specific  heat  ratio,  and  molecu¬ 
lar  weight;  ere  constant  throughout  the  motor. 

5.  The  friction  forces  of  the  combustion  gases  In  the  port 
cavity  are  negligible. 

6.  The  static  pressure  is  constant  across  the  fore-head  section, 
f.e.,  no  static  pressure  loss  resulting  from  mass  addition 
or  area  change. 

7.  The  port  area  and  perimeter  are  constant  across  the  aft-head 
section  from  the  aft  tangent  plane  to  the  grain  exit,  l«e., 
constant  area  duct. 
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Gas  Dynamics  fo»*  Incremental  Control  Volumes  (Continued) 


The  control  volume  or  mass  addition  region  Inlet  Is  defined  as 
station  1  and  outlet  Is  defined  as  station  2  as  shown  In  Figures 
4.1  and  4.2.  The  analysis  of  the  mass  addition  regions  along  the 
grain  segments  Is  accomplished  In  subroutines  AjBSUB  and  AIB’ST 
to  obtain  the  solution  for  the  discharge  pressure,  temperature, 
and  flow  rate  from  known  Inlet  conditions  at  station  1  and  a 
known  value  of  the  Instantaneous  mass  generation  rate.  The  In¬ 
stantaneous  mass  generation  rate  Is  determined  from  the  port 
perimeters,  upstream  burning  rates,  and  mass  addition  region  Incre¬ 
mental  length  as  follows: 


dW 


Pf 


where 


L 


Kfc 


mass  generation  rate,  lb/sec 

port  perimeter  at  station  1,  In 

port  perimeter  at  station  2,  In 

mass  addition  region  length,  In 

solid  propellant  density,  lb/In 

burning  rate  at  station  1,  In/sec 
(determined  from  conditions  at  adjacent 
upstream  Increment  dividing  plane) 

burning  rate  at  station  2,  In/sec 
(determined  from  conditions  at  adjacent 
upstream  increment  dividing  plane) 


A  mass  balance  for  the  slots  between  the  grain  segments  Is  ob¬ 
tained  In  subroutine  SLOT.  The  Instantaneous  mass  generation 
rate  at  each  slot  Interface  Is  determined  as  follows: 

dW  -=  Af  pf  a  Pn 

where 

2 

A*  =  burning  area  at  slot  Interface,  In 

'  2 
P  =  slot  Interface  static  pressure,  lb/fn 

a  =  slot  burn  rate  coefficient 

n  =  slot  burn  rate  coefficient  pressure  exponent 


The  general  energy  equation  and  the  perfect  gas  law  are  applied 
at  station  2  to  obtain  the  state  properties: 
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I  2  9.  *  7  (T.  -  T)  \m 

u  =  I  - - —  ]  (general  energy) 

P  a  12  o  R  T  (perfect  gas) 


The  continuity  and  momentum  eolations  are  iterated  to  obtain 
the  discharge  pressure,  P*,  and  flow  rate,  *L,  for  each  mass 

addition  region.  For  steady  flow  conditions,  the  basic  equa¬ 
tions  are  as  follows: 


12  Pj  Aj  Uj  +  dW  =  12  p2  A?  U2 


(continuity) 


o  ^  p 

P1  A1  +  2  *A2  "  Al'  "  P2  A2 

12  p.  A.  U.  U,  dW 

*  — - -  (lL  -  U, )  + -  (momentum) 

9a  9a 


For  non-steady  flow  conditions,  the  basic  equations  are  as 
follows: 


W2  =  Wl  +  dW-f 


dW  d 

dt  St 


(continuity) 


V 

12R 


a  f  (o  UA)  .  (p2  U2  A2  *  P1  Ut  V 

St  ./  ■iKs-t<lx  +  — : - s - - 


(perfect  gas) 


1 

Pl+P2 

=  p1a1  -  p2a2  +  -y2  iv*j) 


(momentum) 


where 

2 

P  =  static  pressure,  lb/ In 
U  =  velocity,  ft/sec 
V  »  flow  rate,  Ib/sec 
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4.1  Gas  Dynamics  for  Incremental  Control  Volumes  (Continu'd) 


T  =  static  temperature,  R 

2 

A  =  port  cress  sectional  area,  In 

'< 

V  =  port  volume  of  mass  addition  region,  In" 

3 

0  =  density,  lb/ In 
R  =  gas  constant,  °P./ft 


4.1.1  Steady  Flow  Gas  Dynamics 

The  original  program  developed  by  the  Thlokol  Chemical  Corpora¬ 
tion  (references  1  and  2)  based  the  Internal  gas  dynamics  on 
steady  flow.  Continuity  and  momentum  equations  were  Iterated 
to  obtain  mass  addition  region  discharge  pressure,  temperature, 
and  velocity.  This  capability  Is  retained  In  subroutine  A I BSUB. 
The  steady  flow  solution  for  a  segment  slot  Is  obtained  In  sub¬ 
routine  SLOT  by  neglecting  the  time  dependent  terms  of  the  non¬ 
steady  flow  equations.  Therefore  the  steady  flow  option  can  be 
exercised  for  both  monolithic  and  segmented  grain  motor  designs. 

The  solution  of  the  steady  flow  gas  dynamics  Is  obtained  from  an 
Iterative  solution  of  the  continuity  and  momentum  equations  with 
the  general  energy  equation  and  perfect  gas  la w  applied  at  the 
discharge  section  to  obtain  the  density  and  temperature  of  the 
combustion  gases.  The  discharge  flow  rate  is  obtained  from  the 
sum  of  the  mass  generation  rate  and  Inlet  mass  flow.  An  Initial 
value  of  the  discharge  velocity  Is  then  obtained  from  the  momen¬ 
tum  equation,  assuming  P-  -  P.,  and  a  starting  value  of  the  dis¬ 
charge  density  Is  obtained  from  the  continuity  equation.  With 
the  Initial  value  of  the  discharge  density,  the  following  Itera¬ 
tive  procedure  Is  employed  to  converge  the  steady  flow  equations. 

Initially,  a  guess  of  the  discharge  velocity  U2  Is  obtained  from 

the  momentum  equation  by  assuming  no  pressure  loss,  P.  =  P., 
as  follows: 
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Steady  Flow  Gas  Dynamics  (Continued) 

and 


*1 U1 

=  w 
*1 

U2  +  U2  dW 

=  w  u 

*2  u2 

From  the  perfect 

gas  l*w  and 

continuity  equation 

e 

9 

• 

1! 

CM 

=> 

wi  ui 

.  wi  ui 

Hi  ui 

*2 

PoAjlIj 

-  r2a2u2 

12  RT2 

or 

2 

e 

12  RT2 

U2  ' 

Hi  ui 

?2A2 

Using  the  general  energy  equation: 


=%  R  f  <To  -  V 

7  -  1 


12  RT 
=  vi  ui  r2a2 


Multiplying  ooth  sides  of  the  above  equation  by  -2A-  and  adding 
like  terms  to  eliminate  T2  yields: 


*  9„  R 

r-i 


“»  +  <Al+A2'9oR  To  - 


T2  •  (A1+A2)90RT2  = 


(A1+A2)3oR  To 

-M2  VlRT2 
p2a2 


Rearranging  and  combining  terms  yields: 


-2go7(T  -T) 

-  -V-r~~  «2*«VWtVT2> 


M2  *l¥T2 
P2A2 


<A1+A2>9oRT2 


-  (A^2)goRT0  -  - 


Multiplying  through  by  V#2  and  again  rearranging  terms  yields: 
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Steady  Flow  Gas  Dynamics  (Continued) 


.  2<U'(T -TJ 


V'o  2 
'2  7-i 


*•  «■ 
(7-1)  (A.4A?) 

« 

-2A  +  1  L 

^2  +  27 

+ 

■ 

-1 

w> 

(A1+A2)g0ST, 


„  “*  W1  U1  RT2 


P2  A2 


«2  -  <W«»  RTo  *2  "  0 


The  terms  in  the  second  set  of  brackets  may  be  rearranged  with 
the  perfect  gas  law,  remembering  that  P2  ~  Pj! 


2V<W 


y-i 


*  M2  W1  U1 


{7—1) (Aj+A.) 

‘ 

4- 

w2+  2r 

. 

(Al+A2)g0  ?!A2 


U2  -  <V*2>»o  RTo  “2  ‘  0 


Thus,  from  the  general  energy  equation; 

{  r^.j+Aj)  (Aj+A.,) 


‘  2 

-W  U 
2  2 


2A, 


(  2r  2y  )  j 

4  2A2  w.  Uj  j  U2  -  (a1+a2)tor9o  W2  -  0 


(A1-tA2)goA2?1 


The  above  equation  is  solved  for  IL  using  the  quadratic  formula 
as  follows? 


CAA  =  -W, 


2A, 


j  r(At+A2)  (A^) 


‘2  }  2  7  2  7 

CBS  -  'Ai+A2^9oA2P1  +  ^2  *1U1 

CC  =  -(AjtA.2)T0  R  g0  W2 


'?SWS£P  •  .,.»! 


4.1.1 


Steady  Flow  Gas  dynamics  (Continued) 

where  Is  based  on  the  assumption  of  ^  *5  used  onty 

once  In  calculating  the  Initial  value  of 

The  continuity  and  momentum  equations  are  Iterated  by  converging 
on  the  discharge  density  ai  follows: 


1.  P 


W 

L 


tmp 


12  Aa  V 


(continuity) 


•  •  • 
where  dW 


*•  UTmp  *  V!2A2ptap 

(7-1)  U 

7  T  —  T  «.  — .—ilHE 

‘Tmp  o  2a  R  7 

0 

(general  energy) 

2(W,  UT  -  tf.U.) 

4  p  =  P  -  ...  ..1.  Tmp  .1.11 

Tmp  *1  90(Aj-WV2) 

(momentum) 

-  PTmp 

5*  ?2  '  12RTTmp 

(perfect  gas) 

6.  Temp  =  Ptmp  -  p2 

7«  jf  iTgnpJ.  <  .0001  oo  to  8,  otherwise  set  p..  *  p0 

%>  ~  tap  2 

and  return  to  2. 


8-  P2  =  PTmp 
T2  "  TTir,p 
U2  -  UTmp 


9*  Determine  discharge  mach  number. 
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Steady  Flow  Gas  Dynamics  (Continued) 


Non-Steady  Klow  Gas  Dynamics 

The  non-steaay  flow  gas  dynamic  equations  were  developed  to  pre¬ 
dict  Ignition  and  tall-off  transients.  The  non-steady  gas  flow 
equations  account  for  mass,  momentum,  volume,  and  pressure  with¬ 
in  the  control  volume  varying  with  time.  The  fundamental  equa¬ 
tions  of  momentum  and  continuity  along  with  the  perfect  gas  law 
are  expressed  In  partial  differential  form  and  then  Integrated 
across  the  control  volume  describing  an  Incremental  mass  addi¬ 
tion  region  using  the  technique  of  finite  differences. 

The  gas  dynamic  solution  for  the  mass  addition  regions  Is  ob¬ 
tained  by  Iterating  the  continuity  and  momentum  equations  for 
the  discharge  pressure,  temperature,  and  flow  rate.  When  the 
solution  of  the  above  discharge  parameters  have  converged, 
the  discharge  values  of  mach  number,  density,  and  total  pressure 
are  determined.  The  derivation  of  the  time  dependent  momentum 
equation  and  discharge  flow  rate  and  pressure  equations  are  pre¬ 
sented  In  the  following  sections  for  a  mass  addition  region. 
Figure  4.1,  and  a  segment  slot.  Figure  4.2. 

Mass  Addition  Region 


This  section  develops  gas  dynamic  equations  and  solutions  for  a 
mass  addition  region  as  shown  in  figure  4.1.  The  equation  de¬ 
velopment  from  fundamental  enql  Bering  principles  to  obtain 
the  discharge  conditions  Is  take:*  from  reference  6  and  Is  pre¬ 
sented  below.  The  discharge  pressure. 
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4. 1.2.1  Mass  Addition  Region  (Continued) 

and  discharge  flow  rate, 

...  ...  <?1+P2>  V-V  V(P1+P2-P}-Pi) 

U2  «  Wj  +  dW  -  12R  (Tj+Tj)  At  "  lajTj+TjW 

where 

V  =  current  port  volume,  In^ 

3 

V*  -  one  past  time  Increment  volume,  In 

2 

PI  =  Inlet  past  time  Increment  pressure,  lb/ln 

*  2 
P£  *  discharge  past  time  Increment  pressure,  lb/ln 

At  =  time  Increment,  sec 

2 

a  =  vehicle  longitudinal  acceleration,  ft/sec 


are  Iterated  In  subroutine  AIBST  to  obtain  thr  solution  of  the 
gas  dynamic  equations  for  a  mass  addition  region.  Total  pres- 
sure,  static  temperature,  mach  nunber,  end  density  are  then 
calculated  at  the  discharge  station  2. 

The  derivation  cf  the  above  equations  and  method  of  convergence 
follov/s: 

Starting  with  the  continuity  equation: 

1.  ^  (p  U  A)  +  ^  (p  A)  =  0  (continuity  equation) 

Integrating  with  respect  to  x  between  station  1  and  2  gives? 

,  i 

2.  P2  U2  A2  ‘  P1  U1  A1  +  St  J  <*  *■'  ’  0 

1 

The  first  two  terms  are  the  discharge  and  Inlet  mass  flow 
rates,  Wj  and  tfj,,  respectively.  The  tntegral  ts  the  rate 

of  change  of  mass  between  the  stations  1  and  2  and  may  be 
evaluated  as  follows: 

ft  / e  A  d*  *  I?  pm  V  •  * 

1 

where 
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Hass  Addition  Region  (Continued) 


2 

V  =  f  A  dx,  ln^ 

1 

3 

Pm  =  average  density,  Ibm/ft 

dW  =  :nass  flow  generated  in  the  section,  lb/sec 


and  upon  differentiation: 


dV 
^m  dt 


+ 


Combining  terms  and  using  finite  differences  with  the  perfect 
gas  law  yields  from  the  continuity  equation  the  solution  of  the 
discharge  mass  flow  In  terms  of  the  pressures  and  control  volume 


5* 


+  dW  - 


!2P.(T1+T2) 


V-V1 

At 


V 

12R(TyiT^r 


p  +p  ~  pt.pi 
*1*2  1  2 

At 


Euler's  fluid  acceleration  equation  for  unit  mass  is: 


6. 


dU 

St 


^o  dP 

*  ■ 

P  • 


Multiplying  Euler's  equation  (6)  by  pA  and  the  continuity  equa¬ 
tion  (1)  by  U  gives: 


7. 


+ 


p  U  A 


du 

5S 


and 

8„  U  ££  (pA)  +  U  (f  U  A)  »  0 


Adding  equations  (7)  and  (8),  and  combining  appropriate  terms, 
yields  the  one-d i mens  Iona  1  momentum  equation: 


9. 


:d  (p  U  A)  d  (p 


^  -  -A&  . 


-  *  k  (»)  +  p  I; 
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Hass  Addition  Region  (Continued) 


Then,  Integrating  equation  (9)  with  respect  to  x  between 
station  1  and  2  gives: 


10. 


a 

St 


2 

r> 


J 


(oUA) 


dx  + 


<P2U2VP1U1A1> 


-  <P1A1-P2A2)  +  /  PdA 


The  last  Integral  may  be  evaluated  by  defining  a  mean  pres¬ 
sure,  P^  thus: 


2 

“•  /  P  *  -  VA2  -  V 
1 


The  first  Integral  of  equation  (10)  Is  the  rate  of  change  of 
momentum  from  nonstationary  changes  between  stations  1  and 
2  and  may  be  evaluated  by  defining  a  mean  density,  p^,  and 

a  mean  velocity,  and  Integrating  with  respect  to  x,  thus: 


a  K  um  v> 

St  go 


and  upon  differentiation: 

„  a  (pm  Um  V>  Pm  Um  dV  ,  *m  V  d  Um  ,  Um  V  d  *m 

*'  ^  90  ~  %  dt  9o  dt  So  dt 


14. 


Letting  -  0  and  combining  terms  using  finite  differen¬ 

ces  with  the  perfect  gas  law  yields  from  the  momentum  equation 
the  solution  of  the  discharge  pressure  In  terms  of  the  mass 
flows  and  pressures: 


r 


L 


wi  "i 


u. 


piAi  + 


P  +P 
rl  r2 


(A2  -  Aj) 


(Pl+P2)  (V“V,)  V(Ul+U2)  (P1+P2*P1*PP 

2h  90R(t1+T2)  At  ”  24  S0R(Tj+T2)  At 


The  acceleration  term  to  be  added  to  equation  14  Is  derived  In 
Section  4. 1*3* 
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4.1.2. 1 


Hass  Addition  Region  (Continued) 


Equations  (5J  and  (14)  are  Iterated  In  subroutine  AIBST  to  solve 
for  discharge  temperature,  pressure,  and  flow  rate  at  each  mass 
addition  region  as  follows: 

1.  Estimate  a  starting  value  of  the  discharge  pressure  and 
temperature  using  the  Influence  coefficient  equations  for 
constant  specific  heat  and  molecular  weight,  reference  7. 


BT  =  T, (y-l) 


(1+7  Mj)  Mj  dW 

(1-Mj)  Wj 


a)  T,  =  T.  -  DT 

•  2guess  1 


OP  =  Pj  dW  27  Mj 


1  + 


(7-1)  M 


1 


W1  (1“MP 


b>  P2guess  =  P1  ~  0P 


2.  Determine  the  gas  storage: 


dt 


(Pl+P2quess)(2V~V‘>-V(Pl+P2> 
12  R(T1+T2gu.ss> 


3.  Determine  the  discharge  flow  rate: 


.  •  •  dy 

w  =  y  +  dw  -  — 

2  W1  m  dt 


4.  Determine  the  discharge  velocity: 

u  _  w2 1  haatss 

2  P2  A2 


»<**»«*» 


4.1.2. 1 


Mass  Addition  Region  (Continued) 


1 


Determine  the  discharge  pressure  (equation  14  above): 
W,  U\  W,  U, 

P2  “  "H  '  "H  +  Vl  +  ■— W 

a0  30 

-[v(o1+u2)(P1+P2guess  -  P{  -  p>) 

+  (pl+|,2guess>(Ul+U2)(V-V,)]/[2490R 
<VT2gucss>At)  /A2 


If  an  accelerating  reference  system  is  considered  (see 
Section  4.1.3  for  derivation): 


If  >  0.  determine  acceleration  term: 
9/% 


Tenp  .  !<yf“,guess>  (A^2)  ~  * 


P2  =  ?2  +  Temp 


]/f 


4  R<Tl+T2gu.ss>A2 


IP,«P 

"5.„  l3H£s?  <  #ooi,  go  to  step  8,  otherwise  obtain 
P2guess 

new  value  of  P„  using  method  of  false  position  and  re- 
zguess 

turn  to  step  1(b). 


Determine  U2  based  on  converged  P?: 


U2  * 


W2  R  T2ouess 


p2  a2 


Determine  the  discharge  temperature  using  the  general  energy 
equation: 

(7-1)  U2 

T2  "  To  ‘  2  g0  7* 


4. 1.2.1 


Mass  addition  ftiylon  (Continued) 


10. 


If 


T  -  T 

2quess  2 


’2guesj 


< 


new  value  of  T. 

?guess 

return  to  step  1(a). 


.001,  go  to  step  11,  otherwise  obtain 
using  method  of  false  position  and 


11.  Solution  Is  converged,  determine  discharge  mach  number,  gas 
density,  and  total  pressure: 


4. 1.2.2  Segment  Slot  Mass  Addition 

The  development  of  the  gas  dynamic  equations  for  the  region  be¬ 
tween  grain  segments  of  segmented  motors  (referred  to  as  a  slot) 
Is  similar  to  the  non-steady  gas  flow  equation  development  for 
a  mass  addition  region.  Section  4.1.  The  con tro**vol time  for  a 
slot  Is  defined  from  the  forward  slot  Interface  to  the  aft  slot 
Interface  with  mass  addition  occurring  at  each  Interface  and  not 
within  the  control  volume.  The  control  volume  for  a  slot  Is 
shown  in  Figure  4.2. 

The  following  assumptions  are  made  In  subroutine  SLOT  (Including 

the  assumptions  of  Section  4.1): 

I 

! 

a.  The  static  pressure  and  temperature  at  the  slot  Interface  Is 
the  samfc  as  the  port  static  pressure,  -  P2,  P^  =  P^,  T^  = 
Tj,  and  T^  =  T^. 

b.  The  mass  flow  generated  at  the  slot  Interface  Is  a  function 
of  the  port  static  pressure  only  and  Is  determined  from  the 
following  burn  rate  equation: 

t 

!  dW  =  Af  pf  aPn 
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4.1.2.2 


Segment  Slot  Mass  Addition  (Continued) 


where: 


dW  - 
Af  - 

PjT  * 

P  = 


generated  mass  flow,  lb/sec 

2 

burn  area  at  slot  Interface,  In 

3 

sollfil  propellent  density,  lb/I n 

2 

slot.  Interface  static  pressure.  lb/In 
burn  rate  coefficient 


-  burn  rate  equation  pressure  coefficient. 


c.  Static  pressure  at  station  3  Is  a  function  of  the  area  change 
(dA/dx),  and  the  capacitance  effects  (dP/dt  and  dV/dt)  between 
stations  2  and  3#  and  acceleration  of  the  vehicle. 

The  solution  of  the  gas  dynamics  within  a  slot  is  obtained  from 
the  above  assumptions  and  the  equations  developed  In  Section  4.0 
as  follows: 

1.  Determine  the  mass  generation  rate  at  the  forward  slot  Inter¬ 
face  (station  2)  from  the  static  pressure  at  station  1  (port 
cavity  discharge) 

•  a 

dW^.  -  Aj.  Pf  aP 

2.  Determine  the  Inlet  flew  rate,  velocity,  and  mach  number: 

W ^  +  dWj. 

W  R  Tj 

u2  ■  -rrf 


Hj  = 


kk  R  Ti] 


3.  Determine  the  current  slot  volume  and  rate  of  change  of  volume: 


Aj  +  A- 

v  =  ?  (z3  -  z2) 


where 
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4.1. 2. 2 


Segment  Slot  Hass  Addition  (Continued) 


Zj  =  aft  Interface  station  location.  In. 

Z^  -  forward  Interface  station  location,  In. 


dV  it 
dt  *  2At 


VRf2 


D  2) 

SLOTf' 


+  VRr3  +  RSL0Ta' 


wher^ 

T.  =  current  value  of  forward  slot  Interface  distance 
burned  In. 

T  =  current  value  of  aft  slot  Interface  distance 
burned  In. 


4.  Guess  the  value  of  the  slot  aft  Interface  static  pressure  and 
temperature  using  the  Influence  coefficient  equations  for 
constant  specific  heat  and  molecular  weight,  reference  7. 


dT  = 


Tjtr-l) (l+y  M*)  M*  dwf 

*2 


dP  = 


P1  4  90  U*  *2> 


W2(l  -  K*) 


P,  =  P,  -  dP 

3guess  2 

T,  =  T,  -  dT 

3guess  1 


5.  Determine  the  mass  generation  rate  at  the  aft  slot  Interface 
(station  3)  from  the  guessed  static  pressure  at  station  3j 


dW_ 


=  A  a  P" 

a  Sguess 


6.  Determine  the  stored  gas  In  the  slot 

*t  ^P2*P3auess^  dV  V _  (P2*P3guess~P2 

dt  =  ^‘Y^ues^  dt  «  d<YT3gues7 


H.I.Z.2 


Segment  Slot  Mess  Addition  (Continued) 


8. 


where 


dW 

dt 


a  (£JL_) 

St  '12  RT; 


7.  fietsrm* 


•v*  discharge  flow  rate  and  velocity: 


i  +  dw  -Si 
2  a  dt 


U„  = 


W,  R  T., 

3  .  jquess 


P„ 

5guess  3 


Determine  the  aft  slot  Interface  static  pressure  (see  Section 
4.1. 2.1,  equation  14). 


W-U,  W,U,  P.+P, 

-LI  m  JlI  +  p  a  +  -L-iaass*  u  .  A  \ 
g  g  +  P2A2  +  2  <*3  V 


*P2*P3<juess^  (V'V  dV 


2I%  R<T39u.«+T1>  At 


/: 


!  24g  R (rl+f, ) 

3o  '  2  3guess7 


if 


.P  «*p 

3  3quess 


<  CRP,  go  to  step  10,  otherwise  obtain  new 


value  of  PjgUes,  using  method  of  false  position  arid  return 

to  step  5»  If  CRP  is  not  input  the  program  will  set  CRP 
equal  to  ,001. 


10.  Determine  the  slot  aft  interface  static  temperature  frcm 
the  general  energy  equation: 


SM.TCT 


IT  -T 

11.  If  j  3,  .v3u,en?.s  <  CRT,  go  to  step  12,  otherwise  obtain 

1  3 

new  value  of  T^gUCSS  using  method  of  false  position  and 

return  to  step  7.  If  CRT  Is  not  Input  the  program  will  set 
CRT  equal  to  *001. 

12.  Determine  the  slot  discharge  velocity,  mach  number^  density, 
and  total  pressure: 


Non-Steady  Flow  Gas  Dynamics  with  Acceleration 

The  effects  of  longitudinal  acceleration  of  the  vehicle  on  the 
Internal  ballistic  solution  Is  considered  In  this  section.  In 
conventional  gas  dynamic  studies,  the  effects  of  gravitational 
forces  are  not  considered  because,  for  compressible  fluids,  gra¬ 
vitational  forces  are  significantly  less  tiwn  surface  forces. 
Recent  development  of  missiles  for  low  level  ICBH  Intercept  may 
require  boost  accelerations  that  are  of  significant  magnitude  to 
affect  motor  Internal  pressures  and  temperatures.  In  an  accele¬ 
rating  reference  system,  the  force  field  which  results  from  the 
acceleration  Is  equivalent  to  a  gravitational  force  field  In  a 
nonaccelerating  rr*erence  system. 


Non-Steady  Flow  Gas  Dynamics  with  Accalaratlon 


The  acceleration  effects  are  considered  only  on  the  gas  dynamic 
equations  In  the  port  cavfty  of  the  motor  and  not  tn  the  nozzle. 
Effect  of  acceleration  on  a  nozzle  Is  to  move  the  sonic  point 
upstream  of  the  throat.  The  effect  of  acceleration  on  the  motor 
Internal  ballistic*  Is  to  reduce  the  pressure  drop  along  the  pro¬ 
pellant  grain  as  well  as  fore-head  pressure. 

The  resulting  acceleration  term  Is  c^ded  to  the  momentum  equation 
as  follows; 


a  f  Cota)  J;[ ,  (p2  j  a2  -  pi  U1  V 


Ss  J 


L 

(Pj  A,  -  P2  A2>  +  /  P  dA  + 


(p  A  a  iSL) 

m  ii 


where 

2 

a  ~  vehicle  longitudinal  acceleration,  ft/sec 
AZ  «  length  between  Incre*nents,  ft 

p  «  average  gas  density  tn  Increment,  Ib/ft 

m  2 
=  average  cross  sectional  area  In  increment,  ft 
m 

=  conversion  constant,  32.174  lbm/slug 


The  acceleration  tenn  Is  developed  as  follows: 


1.  From  Newton's  Second  Lav/  of  motion: 


Fbf 

Jo? 

II 

where 

Fbf 

»  body  force,  lbf 

w 

»  weight,  Ibm 

a 

*  acceleration,  ft/ sec' 
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Non-Steady  Flow  Gas  Dynamics  with  Acceleration 


2.  The  body  force  exerted  on  the  gas  within  a  mass  addition 
region  Is? 


p  A  AZ 
m  m 


a 


3.  Assuming  perfect  gas  relationships,  the  body  force  may  be 
written  as: 

(PX+P2>  (Ai+V  ^  a 

Fbf  =  ~ 551^  r(Tj+t2) 


4.  a:id  from  P  =  F/A; 

(P.+P„)  (A.-*,)  AZ 
*  =  % 
where 

2 

AP  =  pressure  change  resulting  from  acceleration,  Ib/In 


The  pressure  change  resulting  from  vehicle  acceleration  Is  added 
to  the  discharge  pressure  of  a  mass  addition  region  as  shown  In 
Section  4. 1.2.1,  stop  6,  for  the  Iteration  of  the  discharge  pres¬ 
sure  and  flow  rate,, 
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Complete  Motor  Gas  Dynamics  . 


Solution  of  overall  motor  gas  dynamics  or  Interna?  ballistics 
from  fore-head  to  nozzle  exit  Is  described  In  this  section. 

1*.2.1  Fore-Head  Pressure  Convergence 

Solution  of  motor  internal  ballistics  for  each  time  point  Is 
obtained  by  an  Iteration  process  which  converges  on  fore-head 
pressure,  PH,  (Figure  3*1).  An  Initial  estimate  Is  made  for 
fore-head  pressure  either  from  the  Input  value  PHI  at  time  * 

0,  or  from  the  previous  time  solution  of  PH  at  time  >  0.  The 
fore-head  section  mass  balance  is  obtained  in  subroutine  MNCHN4 
(flow  chart  No.  10)  after  geometry  values  have  been  determined. 
Then  parameters  necessary  to  solve  cylindrical  section  mass 
addition  regions  In  subroutine  SEGSUB  (f lav  chart  number  11)  are 
determined.  When  the  cylindrical  section  Is  complete-  the  aft- 
head  section  mass  balance  Is  obtained  In  subroutine  MNCHN4  and 
the  fore- '.lead  pressure  Is  checked  for  convergence  in  subroutine 
SETPH  (flow  chart  number  12) . 

The  fore-head  pref»rure  convergence  check  In  subroutine  SETPH 
Is  made  as  follows* 

1.  Determine  the  throat  critical  pressure  ratio: 


2.  Determine  the  nozzle  total  pressure: 


-JL 


where 

P  =  aft-head  section  discharge  static  pressure, 
lbs/ in2 

T  =  aft-head  section  discharge  sMtlc  temperature,  *R 
P 

3*  If  PftU  Cst)  >  the  nozzle  Is  choked.  Determine  the 
on  pq  -  a 

sonic  nozzle  flow  rate: 


L 
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4.2.1 


Fore-Head  Pressure  Convergence  (Continued) 


%  NN  At  P0N 
OIS  »  - - xS-— 

C 

where 

NN  -  number  of  nozz?es 
Go  to  5* 

/P.\ 


ft) 


If  I  j  <  Pfi,  the  nozzle  Is  not  choked.  Determine  the 


subsonic  nozzle  flow  ratet 


SDIS  =  NN  PQN  AE£ 


ft)' 


r 


2zL ,  1  1/2 


ToRtr-0 


(-  fell 


where 

=  nozzle  exit  plane  area. 
DIS  «  SDIS 


W  -  DIS 

If  p-| g""1  ■  >  CRW,  estimate  new  fore-head  pressure  as 

follows  (W  -  grain  discharge  flow  rate);  If  CRW  Is  riot  Input 
the  program  *111  set  CRW  equal  to  .005. 


a)  WD  = 


W  -  WD8 

P  -P 
H  Hold 


where  P^j^  =  previous  Iterative  value  of  P^,  psl 

WDB  ==  previous  Iterative  value  of  WD,  lb/sec 
W  =  grain  discharge  flow  rate,  lb/sec 

Dl S  -  MSB 


b)  DEED  = 


?H  “  PHold 


where  DISB  =  previous  Iterative  value  of  DiS 
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Fore-Head  Pressure  Convergence  (Continued) 


c)  If  DEED  =  WD  or  ff  DEED  =  0: 

1.4 

\i  u  \ 

*  ■  PH 

d)  If  DEED  4  WO: 

• 

y  . 

AP  = 


[/ jL\ 

-il 

[  \  DIS  / 

*1 

DiS 

OEED-WD 


e)  WDj  =  V 
DISB  =  DES 


f)  If  «  <  015,  P-In  -  PH 

g)  If  W  >  OIS,  PMX_  -  P, 


h>  PHgu«s  ’  V* 


I)  If  P 


<  P  .  and  P 


Hguess  —  min 


min 


PHguess  *  5‘° 


J)  If  K  <  P  .  and  P  .  /  0 
J  Hguess  —  min  min 


Pu  =  2.0  Pu 

Hguess  H 


k)  If  P  .  <  p  .  >  P  : 

'  min  'fcuess  —  max 


PHguess  *^Pmax~Pmln^  +  Pmin 
Return  to  fore-head  section  with 


PHguess* 


i*  n  1 
*T#A»  A 


For c~nccu  rressure  Convergence  (Continued^ 


1)  If  P  ,  <  Pu  <  P  ,  return  to  rore-head  sect! 

'  min  Hguess  max' 


wIth  PH  =  PHguess' 


4.2.2 


6.  If  F-r-rr--  <  CRW.  convergence  har  been  attained.  If  CRW 
I  D1S  — 

Is  not  Input  the  program  will  set  CftW  equal  to  .005. 

Nozzle  Gas  Dynamics 

After  the  fore-head  pressure  convergence  criterion  has  been 
satisfied  (step  5  and  6  of  the  previous  section),  the  nozzle  gas 
dynamics  In  subroutine  SETPH  are  determined  as  follows: 


1.  Determine  the  nozzle  exit  area: 


Determine  the  nozzle  expansion  ratio: 


€G  A, 


3.  iterate  the  following  equation  for  —  using  the  method  of 
fslse  position:  ON 

i-lr-sn*)*/- J/l-ftfV 


Determine  the  momentum  portion  of  the  thrust  coefficient: 


Zti  ZZ 

*-[ts(*r  -ft)' 


W!>. 
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Nozzle  Gas  Dynamics  (Continued! 


5.  Determine  the  delivered  vacuum  thrust  coefficient: 


C-  .  *  \  “  €r)  C 

foL  '  fo  n  G'  m 


6.  Determine  the  delivered  thrust  at  P  for  sonic  flow: 

a 


F  88  ^CfoL  PQN  At  “  Pa  eG  VNN 


Determine  the  delivered  thrust  for  subsonic  flow: 


Vi*  L/i.!  ;  | 

'-1  I  \?0N/  | 


VE  W  C  NN 
E  N  m 


8.  Determine  the  fore-head  pressure-time  integral: 


/'rf*  * 


Hold  +  PH>  f 


5-  Determine  the  nozzle  total  pressure-time  Integra? 


/P0Ndt  -  /P0Ndt+<B; 


+  P  )  — 

No  Id  r0N'  2 


10.  Determine  the  nozzle  discharge  flow-time  integral: 

A*  =  +  fytold  +  “T 


11.  Determine  the  total  Impulse: 


iT  *  IT  +  (F  +  FoW)  *§ 


I.  A  A 

tf.Z.Z 


4.3 


4.3.1 


nozzle  Gas  Dynamics 


(Continued) 


The  above  iterative  procedure  is  allowed  to  continue  for  no  more 
than  11  iterations,  if  the  fore-head  pressure  h-s  not  converged 
within  11  Iterations,  a  summary  of  the  last  Iteration  Is  printed 
followed  by  a  complete  program  data  dump  and  the  next  case  Is 
set  up.  When  convergence  has  been  attained,  the  program  output 
Is  printed  and  the  next  time  Increment  solution  Is  initiated  by 
setting  P^guess  equal  to  the  previous  time  Increment  solution 

of  P^.  Program  execution  for  each  case  continues  until  *■  -’na¬ 

tion  options  are  satisfied. 


Propellant  Characteristics  and  Burning  Rste  Model 

The  character i sties  of  the  propellant  are  represented  by  a  mathe¬ 
matical  model  of  the  burning  rate  and  properties  of  the  pro¬ 
ducts  of  combustion.  Basic  assumptions  of  a  perfect  gas,  constant 
combustion  temperature  and  constant  specific  heat  allow  the  pro¬ 
pellant  gas  properties  to  be  described  by  the  characteristic 
velocity,  C*,  the  combustion  temperature,  T  ,  the  specific  heat 
ratio,  7,  ar.d  the  gas  constant,  R.  0 

The  propellant  burning  rate  model  allows  either  Isotropic  or 
anisotropic  burning  of  the  propellant  surface.  Isotropic  burning 
Is  defined  as  uniform  combustion  occurring  normal  to  the  propellant 
surface.  In  anisotropic  burning,  the  burn  rate  varies  with 
distance  burned  as  well  as  with  conditions  at  the  gas-propellant 
Interface.  It  results  from  non-homogeneous  dispersion  of  pro¬ 
pellant  additives  near  the  case  wall  and  core  during  propellant 
casting. 

Propellant  Gas  Properties 

The  propellant  gas  properties  may  be  held  constant  or  may  be 
varied  as  a  function  of  the  static  pressure  In  the  port  cavity. 

If  the  gas  properties  are  to  be  held  constant,  the  parameters 
V  C*,  7,  and  R  are  Input*  If  the  gas  properties  are  to  be 

varied,  tables  of  the  combustion  temperature  (TCOMB),  th  roole- 
cul?'*  weight  (AMWG),  the  specific  heat  ratio  (GAMAG),  e  the 
characteristic  velocity  (CSTR),  are  input  as  a  function  of  static 
pressure  (PRESS). 

When  the  gas  tables  are  Input,  a  spline  Interpolation  procedure 
Is  used  In  subroutine  CSTRSB  to  obtain  the  gas  properties  for  the 
static  pressure  at  the  Increment  dividing  plane.  The  spline 
interpolation  procedure  sets  up  a  series  of  piecewise  cubics 
between  the  table  values  to  obtain  an  Interpolated  value  corres¬ 
ponding  no  this  static  pressure. 
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Propellant  Burning  Rate  Model 


The  propellant  turning  rate  model  may  Include  the  effects  of 
erosive  burning*  Erosive  burning  Is  defined  as  the  change  In 
the  local  burning  rate  resulting  from  gas  velocity  parallel  to 
the  burning  surface. 

The  propellant  burning  rate,  Rg,  at  any  Increment  dividing  plane 
Is  determined  from  the  following  parameters  at  the  adjacent 


upstream  Increment  dividing  plane: 

1.  Static  pressure  P 

2.  Gas  velocity  U 

3.  Mass  velocity  per  unit  area  G 

4.  Distance  from  stagnation  point  •v. 

5.  Burning  rate  Rg,jl 

6.  Solid  propellant  density 


With  the  exception  of  p^,  these  values  are  not  Input  but  aro 
calculated  within  the  program. 

Fifty-one  constants  are  available  to  define  the  burning  rate 
equation.  Only  the  constants  that  are  required  for  the  particular 
burning  rate  equation  to  be  used  are  tnput.  These  constants  are 
as  follows: 

KG!  through  KG5 

KUi  through  KU5 

KR1  through  KR39 

KSL0T1  through  KSL0T2 

prior  to  calculating  burning  rate,  critical  values  of  velocity, 
UCR,  and  mass  velocity  per  unit  area,  GCR,  are  obtained  as 
follows: 

UCR  =  KU1  +  (KU2)  P^KU3^  +•  (KU4)  P*KU5* 

GCR  *  KG1  +  (KG2)  P^K63'  +  (KG4)  P^6^ 

Propellant  burning  rate  Is  then  calculated  In  one  of  two  ways* 

If  G  I 3  greater  than  or  equal  to  GCR,  and  U  is  greater  than  or 
equal  to  UCR,  the  following  relation  is  used: 

RB  *  KR1  +  (KR2)  P^3'  +  (KR4)  P(KR5>  +  (KR6)  U^7* 

+  (KR8)  U<«»  -i  (KR 10)  G^11^  +'  (KR12)  G^13* 


4.3.2 


Propellant  Burning  Rate  Mod>' 


♦  (KR 14)  P*1®^  U(iaU6>  *  (KM?)  P(KR!8)  U*1®19* 

+  (KR20)  G (KR22)  +  (KR2H)  P^2^ 

_ KR26 _ 

(KR27)  P  ^2^  +  (KR29)  P^30^ 

^  RBH1  °f' 

,.{^32)  "  G 

*  {KR31)  771^337  e 

If  the  value  of  G  Is  less  than  GCR  or  if  the  value  of  U  is  less 
than  UCR,  then  the  following  relation  is  used: 

RB  =  KR35  +  (KR36)  p^37^  +  (KR38)  P^1®39) 


To  prevent  the  values  of  GCR  and  UCR  from  being  used  simultaneous¬ 
ly  for  choosing  the  burning  rate  model,  one  of  these  values  must 
always  be  equal  to  zero.  This  Is  accomplished  by  setting  the 
values  of  constants  KG1,  KG2,  and  KG4  or  the  values  of  constants 
KU1,  KU2,  and  hU4  equal  to  zero.  The  program  will  automatically 
stop  If  any  of  the  terms  XGl,  KG2,  or  KG4  are  not  equal  to  zero 
when  any  of  the  terms  Kill,  KU2,  or  KU4  are  also  not  equal  to 
zero.  In  such  cases,  the  program  will  print-out  a  statement  that 
the  GCR  or  UCR  coefficients  are  invalid.  If  KR26  Is  not  zero, 
then  neither  KR2/  nor  KR29  can  be  negative  or  simultaneously  equal 
to  zero;  If  this  restriction  Is  exceeded,  the  p-ogram  will  stop 
and  prlnt-out  a  statement  of  Invalid  KR27  or  KR29. 

The  burning  rate  at  the  segment  slot  face  Is  calculated  by 

RBSLOT  =  (XSL0T1)  p(KSLCT2) 

where  P  Is  the  static  pressure  In  the  port  at  the.  segment 
Interface. 


4.3*3 


Anisotropic  Propellant  Burning 


Anisotropic  propellant  burning  capability,  where  burn  late  de¬ 
pends  on  distance  burned,  was  added  to  the  program  because  of 
Boeing’s  experience  wtth  the  HIBEX  motor.  Anisotropic  burning 
occurred  during  both  Ignition  and  tall-off.  It  appeared  to  be 
the  result  of  two  effects;  1)  variation  In  the  alignment  of  the 
staples  between  the  bulk  of  the  propellant  at  the  cesa  wall  and 
the  core  Interface.,-  and  2)  the  burning  distance  required  to  de¬ 
velop  “coning"  about  the  staples.  Anisotropic  burning  Is  most 
easily  represented  by  variation  of  the  constant  "a"  *-,s  a  function 
of  distance  burned  In  the  burn  rate  equation,  r  =  aP?i.  During 
Ignition,  mass  generation  is  determined  by  muitlj.lylrg  the  port 
perimeter  by  the  anisotropic  burning  rate.  During  tnii-off, 
regions  exist  where  both  Isotropic  and  anf  sotropl*.-  burning  occur. 
The  port  perimeter  Is  subdivided  accordingly.  The  total  nass 
generation  Is  then  the  sum  of  the  Individual  mass  csneratlon 
rates. 

The  following  assumptions  have  been  made  In  developing  the  mathe¬ 
matical  model: 

1.  The  anisotropic  region  at  both  the  co'-e  Interface  and  the  case 
wall  Is  of  uniform  thickness  along  the  motor  length. 

2.  The  thickness  of  the  anisotropic  region  Is  the  same  at  both 
the  core  Interface  and  the  case  wall. 

3.  The  burn  rate  variation  through  the  anisotropic  region  Is  a 
function  only  of  distance  burned  and  local  static  pressure, 

r  =  e(T)Pn. 

4.  The  anisotropic  burn  rate  Increases  from  the  core  Interface 
toward  the  Isotropic  propellant  and  decreases  from  the  Iso¬ 
tropic  propellant  toward  the  case  wall. 

5.  The  vore-head  and  aft-head  burning  rate  during  motor  tail- 
off  Is  the  same  as  the  adjacent  tangent  plane  Isotropic 
burning  rate. 

The  following  limitations  apply  to  the  program: 

1.  Anisotropic  burning  cannot  be  considered  for  propellants  with 
wagon  wheel  grain  configurations  during  tall-off. 

2.  Anisotropic  burning  may  be  considered  for  nsn-steady  flow 
options  only  during  Ignition  and  tall-off. 

Program  simulation  of  anisotropic  burning  Is  accomplished  by 
altering  normal  program  solution  during  the  Ignition  transient 
Interval  to  solve  for  the  burn  rate  coefficient  wit'  a  fixed 
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4.3.3 


1 

Anisotropic  Propellant  Burning  (Continued) 


value  of  the  fore-head  pressure  at  each  time  Increment.  The 
burn  rate  coefficient  Is  stored  In  a  table  as  a  function  of 
distance  burned  at  a  desired  location  within  the  cylindrical 
section  (Input  NiNCPL).  A  program  option  Is  available  to  input 
this  anisotropic  burn  rate  coefficient  table  and  solve  for  fore¬ 
head  pressure  as  d*-tcus.*ed  Jn  Section  4.2.1.  The  burning  rate 
within  the  fors-he«.n  and  aft-head  sections  may  be  specified  by 
Inputs  KRH  and  KRN  or  by  the  anisotropic  burning  rate  table 
Inputs.  During  the  tall -off  Interval,  when  the  burning  surface 
Is  within  the  anisotropic  propellant'  region,  the  burning  rate  at 
each  Increment  dividing  plane  in  sectors  6,  ?,  and  8  becomes  a 
function  of  the  distance  from  the  case  wall  within  the  sector  as 
shown  in  Figure;;  4.4,  4.5,  and  4.6. 


The  program  method  of  solution  for  anisotropic  burning  during  the 
Ignition  transient  interval,  when  the  fore-head  pressure  trace 
Is  Input,  is  altered  to  converge  on  the  burn  rate  coefficient, 
KRST.  At  each  time  Increment,  an  estimate  of  the  burn  rate  coef¬ 
ficient  Is  determined  In  subroutine  RRSTSB  from  the  fore-head 
pressure  rise  rate  and  the  motor  configuration  as  follows: 


From  the  perfect  gas  law  using  finite  differences: 


1. 


"ft 

dt 


a 

St 


12  VH  R  T  ‘ 

V 


12  R  T 
V 


and. 


*  • 

2.  =  W.  -  W  . 

dt  In  out 


where 

Wjn  =  generated  weight  flow  rate,  lb/sec 
WQut  =  nozzle  discharge  flow  rate,  lb/sec 
V  =  free  volume,  in^ 


The  nozzle  discharge  flow  rate  is  determined  from  the  nozzle 
geometry: 


9o  At 


where 


4.S-.3  Anisotropic  Propellant  Burning  (Continued) 

g  =  gravitational  constant,  ft/sec* 

°  2 
A+  *  nozzle  throat  area.  In 

Pnu  =  nozzle  total  pressure,  lfe/fn2 

ON 

if  2 

C  =  characteristic  velocity,  ft/sec 

and  the  generated  weight  flow  rate  Is  determined  from  the 
motor  configuration: 


4.  WJn  =  Rfa  Ab  pf 

where 

R,  =  burn  rate.  In/sec 

®  2 
A,  =  total  turn  area.  In 

“  3 

Pf  =  propellant  density,  lb/ In 


5*  RT 
where 

r2 


P?  c  ^2 


21 I 


7(^)' 


r-3 


Therefore: 
dP. 


6. 


H 


dt 


=  P..  = 


12 


r*  c*2  ,  12  r2  c* 


vgr 


(ftbAbpf^  "  V 


^P0N  V 


Combining  and  arranging  terms  with  P^  =  (TPR)P^,  where  ~pr  I 
an  Input  estimate  of  the  port  pressure  drop,  we  have: 

•  (l2 1,  tl  K  m  ph) 

V 


7. 


PH  + 


12  I*  C*2  A 


b('f 


SoV 


and 


AnlsotropH  Propellant  Burning  (Continued) 


After  the  Initial  estimate  of  KRST  Is  made,  the  ballistic  solu¬ 
tion  Is  converged  for  the  fixed  fore-head  pressure  obtained  from 
the  Input  pressure  trac*.  using  the  convergence  procedure  outlined 
In  Section  4.2.I.  When  the  anisotropic  burn  rate  coefficient 
table  Is  Input,  the  method  of  solution  remain  unchar)  ged  except 
that  the  burn  rate  coefficient  which  depends  •  distance  burned 
at  location  N1NCPL  frcfii  the  forward  tangent  pi.  )e  is  determined 
from  the  Input  table  at  each  time  increment. 

The  Ignition  transient  Interval  Is  terminated  when  the  value  of 
time  exceeds  the  last  table  valu'i  of  the  Input  fore-head  pressure 
trace  Independent  variable  TlMEPh:(NPH),  or  when  the  burn  rate  table 
Is  Input,  the  termination  option  TST.  The  steady  state  Interval 
will  then  continue  with  the  last  table  value  of  the  burn  rate 
coefficient  dependent  variable  AKRTAU(NAKRST)  or  the  burn  rate 
coefficient  Inputs  KR{2)  and  KR (3&)  (depending  on  choice  of 
Inputs)  In  the  general  burning  rate  equation. 

As  the  burning  surface  progresses  toward  the  case  wall,  the  ani¬ 
sotropic  region  Is  first  exposed  In  the  region  of  sector  8  as 
shown  In  Figure  4.4.  This  results  In  a  non-uniform  burning  rate 
along  the  propellant  burning  surface  during  the  tall-off  Interval. 

Three  turning  rates  are  determined  for  the  burning  surface:  R^g 

which  Is  determined  from  the  anisotropic  burning  distance  In  sec¬ 
tor  8,  R^  which  Is  determined  fY.-xn  an  Integration  along  the 

anisotropic  perimeter  of  sector  7  between  the  Isotropic  and  ani¬ 
sotropic  burning  distances  using  the  arlsotroplc  burn  rate  coef¬ 
ficient  table,  and  the  normal  Isotropic  burn  rate,  R^. 

Two  separate  burn  distances  are  defined:  the  Isotropic  burn 
distance  In  sectors  1  through  7#  TAUZ(III),  arid  the  anisotropic 
burn  distance  In  sectors  7  and  8,  7AUZTC (III).  The  anisotropic 
burn  rate  becomes  progressively  less  during  motor  tall-off,  re¬ 
sulting  In  an  anisotropic  burn  distance  less  than  the  Isotropic 
burn  distance  and  producing  burning  that  Is  not  normal  to  the 
grain  surface. 

Figures  4.4,  4.5,  and  4.6  show  the  configurations  of  the  aniso¬ 
tropic  propellant  region  that  can  exist  at  the  case  wall  for  a 
standard  star  configuration  wtth  an  Inert  sliver.  The  angles 
,  1)22*  ar,<*  "ANGLE"  are  used  to  determine  the  anisotropic  pro- 
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pellant  perimeter  lengths,  AL7  and  ALg,  during  motor  tall-of 

Angle  tj2  Is  subtended  from  the  R^  r?d{«  .  x  and  locates  the 

Intersection  of  the  anisotropic  propellent  region  with  ths  Iso¬ 
tropic  propellant  region.  Angle  r)22  Is  subtended  from  the  motor 

axis  to  the  same  Intersection  point  for  r,2.  "ANGLE1*  Is  sub¬ 
tended  from  the  motor  axis  and  Identifies  the  point  of  Intersec¬ 
tion  cf  tf.e  Isotropic  burn  distance  vector,  swung  from  R_,  with 
the  case  wall.  ' 

Subroutine  LPTO  contains  the  geometry  calculations  to  determine 
the  sector  perimeter  length  of  the  anisotropic  propellant  for 
each  reference  plane  during  motor  tall-off.  The  anisotropic  pro¬ 
pellant  perimeter  length  In  sector  8  Is  Identified  as  AL8,  and  In 
sector  7  Is  Identified  as  AL7*  The  perimeter  length  c-f  sector  7 
anisotropic  propellant  (AL7)  Is  assumed  to  be  a  straight  line 
between  the  point.,  determined  by  the  Intersection  of  the  Isotro¬ 
pic  propellant  with  the  anisotropic  propellant  and  the  anisotropic 
burning  distance  with  the  case  wall  or  with  sector  8. 

The  distance  burned  of  the  anisotropic  propellant  In  sectors  7 
and  8  Is  computed  from  the  burn  rate  of  the  anisotropic  propellant 
In  Sectdr  8  (RB8)  In  Subroutine  SEGSUB.  Once  sector  8  has  burned 
out,  RB8  Is  determined  from  the  first  table  value  of  the  aniso¬ 
tropic  burn  rate  coefficient  (the  minimum  value).  The  progres¬ 
sion  of  the  Intersection  of  the  propellant  with  the  case  well  Is 
assumed  to  proceed  at  the  minimum  burn  rate. 

The  mass  flow  generated  Is  determined  In  subroutine  SEGSUB  from 
the  perimeter  lengths  of  the  anisotropic  burning  propellant 
(AL7  and  Al8),  and  the  Isotropic  burning  propellant  (ALP)  and 
their  cor  respond  trig  burn  rates  as  follows: 

dHj  -  NO  pf  £2.  (ALy  Rgy  +  AL^j  Rg^j) 

dWg  =  NO  pj.  CZ  (ALg  Rgg  +  Alggi  rb8hP 
OW  -  (ALPK|  -  AL^,,)  ABH|  -  (ALP  -  A4C>  A6 

+  dWy  +  dWg 

where 

ALto  -  (AL7  +  ALg)2N0 
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and  subscripted  "H!“  values  are  at  'he  Inlet  to  the  mass  add!" 
tlon  region.  The  non-subscrlpted  ones  represent  the  outlet  to 
the  mass  addition  region. 


Figure  HA.  Mathematical  Model  of  Mass  Addition  Region 
Control  Volume 


©  Slot  foreword  Interface 

O  Slot  forward  Interface  mass  addition  re9lon 
©  S>«  .ft  lnt.,f.c.  m„  .dd,tlori  ft3lon 

®  Slot  aft  Interface 


Figure 


Of  Slot  Cr.|n 
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i.  Input  Reference  Planes 

I 


Figure  4*3  Reference  Plane  and  IncrcmcntDlvIdlng  Plane 
Identification 


5.0 


GEOiiKTRICAL  DEFINITION  OF  PROPELLANT  GRAIN 


For  geometrical  analysis,  the  motor  is  divided  longitudinally 
into  th>'ee  sections:  fore-heas  section,  cylindrical  section,  and 
aft-head  section.  The  '^re^head  section  may  Incorporate  a  head- 
end  with  web  or  it  may  he  defined  In  a  manner  similar  to  the 
aft-head  sectJor  which  Is  a  straight  through  grain*  The  follow¬ 
ing  describes  subroutines  and  equations  used  to  calculate  burning 
surface  area,  port  cross  sectional  area,  propellant  volume,  cen¬ 
ter  of  gravity  and  moment  of  inertia  for  the  fore-head,  cylindrical, 
ond  aft-head  sections. 

5.1  Grain  Cross  Section  Geometry 

Propellant  grain  cross  section  options  which  are  programmed  for 
the  computer  are  shown  In  Figure  5.1.  The  grain  design  may  vary 
from  the  more  complicated  forked  wagon  wheel  to  the  standard  star, 
or  to  the  circular  port.  The  slotted-cone  grain  design  ts  a  modi¬ 
fication  of  the  standard  star.  Variations  In  the  Input  parameters 
for  the  forked  wagon  wheel,  which  are  shown  in  Figure  5*2,  will 
produce  the  five -basic  grain  design  patterns.  The  baste  parame¬ 
ters  required  to  describe  the  gratn  cross  section  are  input  for 
reference  planes  located  at  specified  distances  from  the  forward 
tangent  plane  of  the  motor,  Figure  4.3.  The  refeience  planes 
desiribing  the  propellant  cross  section  may  be  placed  at  any  de¬ 
sired  location  within  the  cylindrical  section,  thus  allowing  ac¬ 
curate  descriptions  of  the  propellant  configuration  for  either 
monolithic  or  segmented  motors. 

5.1.1  General  Forked  Wagon  Wheel 

Cross-sectional  geometry  constants  are  required  to  determine  the 
perimeter  length  and  port  area  of  each  Input  reference  plane. 

These  basic  geometry  constants  are  determined  In  the  first  core 
load.  Once  these  plane  constants  have  been  computed  and  stored 
for  each  Input  reference  plane,  the  perimeter  length  for  all 
values  of  distance  burned  and  the  Initial  port  area  are  deter¬ 
mined  for  each  reference  plane  in  the  second  core  load  and  stored 
In  tables  for  use  during  the  internal  ballistic  solution  In  the 
fourth  core  load.  The  reference  plane  constants  for  the  forked 
wagon  wheel  grain  configuration  are  determined  by  subroutine 
PLNCNS  from  the  geometry  inputs.  Figure  5.2.  These  plane  con¬ 
stants  are  shown  In  Figures  5.3  and  5.4.  A  total  of  45  constants 
are  generated  and  stored  for  each  reference  plane: 

Subroutine  PLNCNS  Reference  Plane  Constants 


R 


1 


L 


SI 


R 


3 
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i 


I 


* 


'1  "  NO 


sin  a 


T _  =  R„  +  l 


01 


2max  2  A  cos  a 


01 


sin  a, 


T,  ...  =  T. _ +  L 


02 


4max  ~  2max  B  cos  a, 


02 


cos  a, 


02 


T5max  ^T4max  +  ^  cos  a,.}  "k 


-  R. 


03 


*«  '  (TS™x  +  V  “S  *03 


Yw  =  «!  ♦  RjO-s!*  a01)  *  La  cos  O01  +  t2|mx(sIo  ofo- 


+  Lg  cos  aQ2  +  R^(sln  aQ3  -  sin  aQ2)  -  R5  sin  a( 
K  -  \  -  «5)2-(*45  «*  “o3  ‘  V45  sl"  “o/j 


-  <45  sin  t»03  -  Y45  cos  k03 


^2max  "  V 


cos  +  Tan 


03 

V03  =  R1  +  R2  +  <R3-R2>  sIn 


/(«0i  -  «02) 

\  2 


•jsln  a( 


°0i+  p; 


2s»*-|l3)Tan 


*01 


+  laJ  C0saoi 

03  = 

{  03  Y03  ) 

nc  ~ 

(T^Y  +  W  C0S  “l 

m 
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'm 

1 


i 

5 


i 

i 


y05  ‘  V45  +  (R5  '  V  s,n  a03 


.  f  ¥  2  a.  V  2  1/2 

R05  *  l  X05  YG5 


X07  =  x45  +  LC  s,n  a03 


Y07  =  Y45  +  LC  cos  a03 


07 


=  f  X  2  +  Y  ^  ^ 
1*07  t07  J 


sin  a 


T  -  X  «f  I  —  - 

Smax  5mex  C  cos  a. 


22. 


03 


7max 


-  [  x072  *  (Rf  '  Y07)2  ]  U2  -  aS 


sIn  a05 
Ro  +  Lr  - « 


12max  8  E  cos  a. 


05 


X76  =  *T12max  +  V  cos  ^04 

Y76  =  R9  +  R8^"sln  °W  +  LE  c0S  a05  +  T12max^sln  a05 

-  sin  Ofo)  -  R6  sin 

lD  =  [^Rf  "  Tw  “  *6-  "  ^X76  cos  a04  “  Y76  s,n  j 

-  x76  sln  a04  -  Y76  005  a04 

X0$  =  X76  +  L0  sln  a04 

Y09  “  Y76  +  l0  cos  a04 
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X  ^  +  Y  M  ^ 
*09  t09  I 


a05  ‘  a04 


sin  aQ5  +  cos  aQ5 


®nc  “  aQL 

XQ11  *  ^12max  “  ty  [7an  2  sln  a05  + 

Y011  =  R9  +  R8  +  (R7“R8^sln  a05  +  [^T12tnax  "  V 
Tan  ("to  — j  ♦  4  ]  cos  0^ 

“oil  -  [X0112  +  Y0112] 1/2 


+  4  cos  or05 


T, _  =  T, „  +  L»  s5n  O04 

lOmax  12max  D 

cos 


9max 

= 

[X092 

+  (Rf 

'2  - 

arc 

cos  | 

'!2Z\ 

^  R07  / 

*3  = 

arc 

cos  j 

ft) 

k2  1  1/2 


*4  ^  81  "  02  “  ®3 


<W  ■  arc  C01  !<T7max  *  R5)2  '  (Rf  "  V07)2  j 


(t_  +  r_)  -  a., 

'  7max  5  03 


^?2max  =  2  *  02  ~  P71max  -  C03 


P91max  ~  arc  cos  [^T9max+ V  "*Rf  **  W  j  /  /^Smax+R6' 
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P92max  =  2  +  °3  “  P91max  "  a0k 

Tmax  =  T7max*  or  T9max  •  whlchev*r  Is  larger. 

If  80  Is  Input,  determine  for  slotted-cone  geometry  as  fol¬ 
lows  (Figure  5*7) 2 


L  AJB 

arc  sin 

n*f  -  *s  -  v  »•"(»,)! 

L  T6m  +  R5 

J 

^1A2C 

= 

it  -  82  - 

Zajb 

.^ACB 

arc  sin 

(Rf  -  «5  -  Tw)  sln(ZABC) 

L  Rf 

J 

^:bac 

= 

it  -  ZABC  -ZACB 

GFE 

= 

I "  aoi 

<£AFE 

= 

2  +  °01 

ZAEF 

src  sin 

"(Rj  +  R2)  sln^AFE)" 

-  Rf  ~ j 

^GAE 

= 

it  AFE  AEF 

0  >  <^GAE, 

T  =  F 

max 

lf  ~R1 

If  eQ  <^:GAE, 


BH 


=  fcf  +(Rf  -  *5  -  Tw)  -  2Rf(Rf  -  r5  -  xw)  cos(e0+e2)] 

1  f  (0Q  +  e2)  >  ^BAC, 


Z.4H2  =  arc  s In 


^  abh  =  it  -  Z  ahb  -  eQ  -  e2 


(Rf  "  R5  "  XJ  s5n<8f,  + 
BH 


ZHBC  =  ^ABC  -  ZABH 
BH*  =  BH  cos  UHBC) 
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T  =  (R-  -  R.)  or  (BK1  -  Rc),  whichever  Is  smaller, 

max  T  i  p 

if  (e0  *  e2)  <^bac, 

Tmax  ~  (Rf  “  or  **  R^)>  whichever  Is  smaller. 

The  calculated  plane  constants  appear  on  the  printout  and  allow  the 
user  a  means  of  checking  If  the  proper  grain  design  Is  being  ana¬ 
lyzed.  Although  the  program  Is  designed  to  solve  the  general  con¬ 
figuration  shown  In  Figure  5*2,  there  are  certain  variations  of 
this  configuration  that  exceed  the  mathematical  limits  of  the 
analysis.  To  obtain  a  program  solution,  all  of  the  following  con¬ 
ditions  must  <sxlst  for  each  reference  plane  used: 


R1 

must 

be 

greater  than  or 

f.qual 

to 

zero 

*3 

must 

be 

less  than  or  equal  to 

T2max 

R7 

must 

be 

less  than  or  equal  to 

T12max 

*9 

must 

be 

greater  than  or 

equal 

to 

2ero 

Lc 

must 

be 

greater  than  or 

equal 

to 

zero 

ld 

must 

be 

greater  than  or 

equal 

to 

zero 

P71M 

.mist 

be 

greater  than  - 

equal 

to 

zero 

P91M 

must 

be 

greater  than  or 

Q 

equal 

to 

zero 

“oi 

must 

be 

less  than  90 

must 

be 

greater  than  or 

equal 

to 

zero 

“03 

must 

be 

less  than  SO0 

°bk 

must 

be 

greater  than  or 

equal 

to 

zero 

0Ua 

must 

be 

less  than  90° 

*u 

must 

be 

greater  than  or 

equal 

to 

zero 

If  any  of  these  restrictions  are  exceeded  for  any  of  the  Input 
reference  planes,  the  program  will  automatically  stop  and  print 
the  reference  nlane  dimensions  along  with  a  statement  of  the 
exceeded  restriction.  Although  the  program  Is  self-checking  for 
these  mathematically  Invalid  configurations,  there  are  some  phy¬ 
sically  invalid  configurations  for  which  there  Is  no  such  check. 

Special  attention  should  be  given  to  the  manner  In  which  input 
lengths  LA}  LB'  and  Lj.  are  defined.  The  length  Is  measured 

along  a  common  tangent  to  the  arcs  generated  by  radii  R2  and  R_ 

(Figure  5»2).  One  of  the  points  that  defines  this  length  l-  the 
point  of  tangency  with  the  arc  generated  by  the  radius  R2«  How- 
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ever,  the  other  point  that  defines  length  Is  not  the  point  of 
tangency  with  the  arc  generated  by  radius  R^.  This  end  limit  of 
length  L^  Is  determined  by  bisecting  ?*he  Include!  angle  between 
lengths  L^  and  Lg,  end  extending  this  bisector  Mntll  It  Inter¬ 
sects  the  axis  of  symmetry  of  the  grain  sector.  F/om  this  Inter¬ 
section  point,  a  line  perpendicular  to  L^  Is  drawn.  The  Inter¬ 
section  of  this  perpendicular  and  the  line  Is  the  point  which 

defines  the  end  of  length  L*.  Line  lengths  Lg  and  Lr  are  defined 
In  a  similar  manner. 

The  perimeter  length  and  Initial  port  area  of  the  grain  cross 
section  for  each  reference  plane  are  determined  In  the  second  core 
load  for  Incremental  distances  burned  before  and  after  web  burn¬ 
out  specified  by  Input.  The  geometry  plane  constants  of  each 
reference  plane  are  moved  from  storage  Into  working  locations  In 
subroutine  LPOAPS  and  the  perimeter  length  and  Intttal  port  area 
are  determined  for  each  sector  of  the  cross  section  (shown  In 
Figures  5.2,  5-3,  and  5*4)  at  the  specified  DTAU  Increments  by 
subroutine  AFP5UB. 


Subroutine  LPDAPS 

Subroutine  LPDAPS  sets  up  the  parameters  required  to  determine 
the  perimeter  length,  Lp,  and  the  cross-sectional  propellant 

area,  App,  for  both  the  prtmary  and  secondary  propellant  tips  of 

the  forked  wagon  sheel.  Initially,  the  parameters  required  to 
determine  the  secondary  propellant  tip  are  set  and  subroutine 
AFPSUB  Is  used  to  evaluate  the  perimeter  length  and  the  area  of 
propellant  for  each  sector,  and  then  the  parameters  for  the  primary 
propellant  tip  are  set.  Figure  5.5  shows  the  boundaries  of  the 
sectors  for  the  forked  wagon  wheel  grain. 

The  perimeter  length  and  propellant  area  of  sector  8  are  then  de¬ 
termined  and  summed  with  the  values  determined  from  sectors  1-7 
and  9-13.  The  total  perimeter  of  the  grain  configuration  !s  ob^ 
talned  by  multiplying  the  sum  by  2N0.  The  port  area  is  obtained 
by  summing  the  cross  sectional  area  of  the  propellant  sectors,  mul¬ 
tiplying  by  NO,  and  subtracting  the  result  from  the  area  of  a 
circle  of  radius  Rp. 

Subroutine  AFPSUB 

Subroutine  AFPSUB  determines  the  port  perimeter  length  of  one- 
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half  of  a  symmetrical  section  and  propellant  cross  sectional  area 
of  a  symmetrical  section  of  all  sectors  except  sector  8.  The 
plane  constants  required  for  each  sector  arc  set  by  subroutines 
MNCHN2  and  LPDAPS. 

The  perimeter  and  area  of  the  cross  section  for  all  grain  confi¬ 
gurations  arc  calculated  from  basic  trigonometric  formulas.  Re¬ 
quired  angles  and  line  lengths  are  computed  from  known  dimensions 
using  the  law  of  sines  and  the  law  of  cosines.  Propellant  cross 
sectional  nrea  Is  determined  by  adding  and  subtracting  areas  of 
circular  sectors  and  triangles.  All  configurations  contain  only 
straight  lines  and  circular  arcs. 

Initially,  the  constants  to  determine  the  Initial  area  and  peri¬ 
meter  length  for  sectors  9  through  13,  Figure  5«5,  are  set  by 
subroutine  LPDAPS,  and  then  the  sectors  J  through  7  are  set  such 
that  L^j  —  Lj,  =  L2,  =  =  ar,d  =  L^.  The 

sector  perimeter  length,  L,  and  Initial  propellant  cross  sec¬ 
tional  areas,  AFP,  are  determined  as  follows: 


1.  If  T  <  R2, 

4  -  (R2  -  T) 

AFPj  •*  (R2  -  TjLj 

2.  If  T  >  R2, 

4  =  ° 

AFPj  =  0 

3.  If  T  <  Rv 

*4  ‘  <T2»x  ’  R3!  Ta"  (°°L2--') 
AF0  *  (R3  +  T2*ax  *  2T' 
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4.  If  t>«3  and  lft<tte 

/ <xn . 

^3  "  <t2m*  '  T>  (-£i 
AF0  -  <T2max  *  T>  ^3 

5.  If  *  >  *3  and  If  r  > 

l2"° 

AU  =  0 

3 

ArD  =  0 


6,  If  r  <  r2, 

4  =  la  +  *3 

_  AFP,  =  (R2  -  2T  +  T2max)  La 


7. 


If  T  >  R2, 


f  <  WT>  ' 

^T2max"R2/\ 
^2max  "  ^  ^L2 


+  AL 


8.  Sf  r  <  Ry 

4  «  (R3  -  t) (a01  -  aQ2) 
afp3  -  U(R3  -  t) 


9.  If  T  >  Ry 

Se0 

*V° 
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10. 


11. 


12. 


13. 


14. 

JM 1  i 


4  =  Lg  +  ^2 

AFP4  *  (T2max  “  2T  +  W  4  +  Ai'D 

,f  T2k*x  ^  T  <  TW 

L4  =  'T4max  "  4^  ^T4max  ~  4max^ 

AFP4  =  <T4max  “  T>  4 


|lfT2r«x-<^W 


4  =  0 


AFP4  =  0 


If  T  <  T, 


4msx 


l5  =  (r4  +  T) (a03  -  a02) 


AFFS  ■  +  T5n«x>  tR4  •>  WS,"<°b3  •  a02>-LA  +  T) 


* F  TAmax  ”  T  <  T5max' 


L5  =  (R4  +  T)  ao3  "  arc  cos 


VW  “S<V  ! 

- !5^5 - j 


r  u  i 


AFP5  .  (R„  +  r5mx)  (RJ(  ♦  T)  sin 


Sk  +  T, 


-  »4  +  t)  Ls 
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5.1.1 


l6-  lc 


*"6  '  (TS»ax  -  2T  +  W  LC 


lf  V  <  T  <  W  lC 


~6  ^c^6max  ~  T^'T6max  “  T5max' 

afp6  '  <V^,  ’  T>  L6 


18.  If  -r_  <  t  >  t£  - 

5max  —  &na;r 


!  6  "  0 


AFP6  -  0 


19*  »fT^W 


p71  a  P71max 


*  W(?‘5  *  W^l*  “  ?71(R5  +  T)i 


20.  If  T,  <  T  <  T_  - 
6msx  7«nax' 


0_,  =  +  a.,  -  arc  cos 

71  7 ’max  03 


\  (T  +  Rg)  / 


AFP  =  (R?  +  T7!wx){R5  +  T)stn(P71)  -  P71(R5  +  T)' 


21.  If  TSIVR  >  0, 


rR.,s!n(e2) 

P  4 

'*  ,  A 


V  ’  *r;  slr  i  r5« 

^  5  x 


-  arc  sin 


R75ln(y*SLy) 

(Rc  +  ^slvrT 


Si.V'1 


Tesno 


V  +  -R5  +  W'  “  2  Tw(R5  +  W 


cos 


arc  sin 


/R~  sin  9,  \ 

\v^:) 


5*1*1 


General  Forked  Wagon  Wheel  (Continued) 


Temp  =  arc 


R7  stn  8. 

i,n  vrr-1- 
5  7fnax 


Temp.  =  k  -  arc  sin 


Temp  =  arc  sin 


K s!n  \ 

I  Temp 

|~Temp  sln(Temp( ) 


SLVR  +  R? 


PVXX  =  Tempc  '  Tempa 

ASLVR  =  92Rf2  “  <Rf  “  Tw)R7  stn  e2  *  (02  "  eSLV} 

+  (TSLVR  +  R5)R7  sfn  jrf  -  02  +  6$LV 

arc  sin  [  V^VW'^LV*  *  V] ! 
“  TemP(TSLVR  +  R$)  sIn(*  “  Ternpc  -  Tempb) 

^vxx  +  py x'  ^tslvr  +  Rs^Z 


22.  !f  T<T. 

—  vr 

P72  =  P72max 

AFP?  =  AFP  +  e2Rf2  -  RfR7  sin  02  -  372(r5  +  T)2 
L7  =  ^1  +  *72!  ^R5  +  T) 


23.  If  T  >  T 


=  -*’’c  C2S 


f  (T+R5)2  +  (Rf  -  Tw  -  R.)2  -  Rf2 
(  2(T  +  R5) (Rf  -  Tw  -  R5) 


"  P71max  "  a03  '  2  +  *2 


5*1.1 


General  Forked  Wagon  Wheel  (Continued) 


AFP7  ,  AFP  ♦  02Rf2  -  «fP7  tin  *.  *  A^  +  f>„) 


-  R^  arc  cos  { 1- 


-  P72(R5  +  T)Z 

L?  =  (P?1  +  $n)  (Rg  +  v) 


^P,.«(S!n(P72M||>  P„)  j2)‘/2 


24.  If  T  <  T, 


SLVR' 

L  =  Lj  +  L2  +  L3  +  l.4  +  L5  +  L6  +  L; 

AFP  =  AFPj  +  AFP2  *  AFP^  +  AFP^  +  AFP5  +  AFPg  +  AFPy  -  ASLVR 


25.  If  T  >  *1^ 

L  »  0 
AFP  =  G 


5.1.2  Slotted-Cone 

The  clotted-cone  configuration  Is  an  addition  to  the  standard 
star  In  which  a  segment  of  propellant,  represented  by  angle  0^ 

has  been  Inserted  as  shown  Jn  Figure  5*6.  The  Insert  of  propel¬ 
lant  results  In  additions  to  the  perimeter  and  areo  calculations 
for  a  standard  star  which  requires  the  determination  of  three 
bade  angles,  QAG,  RAC,  and  SAG,  shown  In  Figures  5*9#  5.10,  and 
5.11,  respectively.  Nine  points.  A,  B,  C,  C',  E,  F,  6,  J,  and  H 
are  defined  In  Figure  5.7*  and  are  used  In  the  following  analysis 
to  determine  line  lengths  and  angles.  Four  moving  points,  P,  Q, 

R,  and  S  on  the  burning  perimeter  are  Identified  In  Figures  5.8 
through  5*13. 

Point  A  lies  on  the  motor  axis  (Figure  5*7).  Point  B  coincides 
with  the  center  of  the  Input  radius  R^  and  point  F  coincides  with 

the  center  of  the  Input  radius  R2>  Point  G  Is  located  on  the  case 
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Slof.ted-Cone  (Continued) 


wail  at  the  point  where  a  line  through  points  A  and  F  Intersect 
the  case.  Points  C  and  E  are  located  oh  the  case  wall  where  a  line 
perpendicular  to  side  at  the  Intersection  of  side  with 

radius  R^  Intersects  the  case  and  where  a  line  perpendicular  to 

side  at  the  Intersection  of  side  Lc  with  radius  R2  intersects 

the  case.  Point  H  lies  on  the  case  wall  at  the  angle  0Q  from  point 

G.  Point  C 1  defines  the  point  where  side  Lc  of  the  standard  star 

disappears  as  a  result  of  the  progression  of  the  burning  surface 
of  the  slotted  cone  Insert. 

Perimeter  length  QR,  Figure  5.9,  defines  the  addition  to  side  t-c 

of  the  standard  star,  resulting  from  the  addition  of  the  central 
angle  8g  of  the  slotted-cone.  Perimeter  length  PQ  defines  the 

length  of  the  arc  subtended  from  the  radius,  (R^  +  t),  and  Inter¬ 
secting  with  length  QR.  When  the  burning  distance,  t,  Is  less 
than  or  equal  to  the  Input  radius  R?,  only  the  perimeter  length 

PQ  Is  present.  When  the  burning  distance  Is  greater  than  R2,  but 
less  than  or  equal  to  the  geometry  constant  T*  x,  as  shewn  In 
Figure  5*9,  both  perimeter  lengths  PQ  and  QR  are  present. 


Perimeter  length  RS  defines  the  addition  to  the  ere  subtended  from 
the  point  of  the  Input  radius  R,.  (geometry  point  B).  Whenever 


Is  greater  than  T7max,  and  T7max  Is  greater  than  (Figure 

5.10),  and  T  Is  greater  than  but  less  than  or  equal  to  T7m{|X, 

geometry  point  S  lies  on  the  line  AG  and  perimeter  lengths  PQ, 

QR,  ar.d  RS  are  present. 


Whenever  T  Is  greater  than  T7(wx  but  less  than  T^, 
length  RS  Ts  defined  as  In  Figure  5-11. 


perimeter 


Whenever  T  Is  greater  than  T^,  as  shown  In  Figures  5.12  and  5.13, 

geometry  point  R  and  perimeter  length  QR  do  not  exist  and  only 
perimeter  lengths  PQ  and  QS  remain. 

The  angle  8q  Is  allowed  to  vary  on  Input  between  0°  and  90°. 

Should  0Q  be. less  than  the  angle  QAG  shown  In  Figure  5.9,  then 

the  perimeter  length  PQ  does  not  exist  and  point  Q  then  lies  on 
the  line  AH.  Should  8^  be  less  than  the  angles  QAG  and  RAG,  then 

the  perimeter  lengths  PQ  and  QR  and  the  points  P  and  Q  do  not 
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5.1.2 


Slotted-Cone  (Continued) 


exist  and  the  point  R  lies  on  the  Itne  AH.  Wh*n  angle  SAG  ex- 
ceeds  0^  complete  burnout  has  occurred. 

The  following  are  the  port  perimeter  and  propellant  cross  sectioned 
area  equations  for  this  addition  to  the  standard  Jtar  donflgun- 
tlon.  Figures  5*8  through  5.13  show  the  progression  of  the  burn¬ 
ing  surface  and  the  geometric  figures  required  for  each  equation. 


1.  If  T  <  R2  (Flsure  5.8), 
Perimeter  =  ®q0*1  +  T) 


Area  =  0Q  JRf 


-  (Rx  +  *)2] 


2.  If  r  >  R2  (Figure  5.9 ), 


L  AJB  =  arc 


...  ...  J 


L  ^6max  +  R5 


L  ABC  ■  *  -  02  -  L  AJB 


L  ACB  *  arc  sin 


(».|.  **  P.j 

“  T  ) 
vr 

Rf 

sin  L  ABC 


L  BAC  »  x  -  L  ABC  -  L  ACB 
sin  L  BAC 


BC  =  R 


Tfe 


TV7 


f  sin  L  ABC 
BC  -  Rr 


(Rf-R5-Tw)2-R12>^52-2(Rf^s-Tw)Rs  cos  L  ABC 


2Rj  -  2Rg  +  2(Rf-R5-Twy  cos  /  ABC 
Tempyg  «  minimum  of  Tfo  and  Ife 


5*1*2 


Slotted-Ctne  (Continued) 


3*  If  R2  <  Tenp^,  and  If  T  <  Tgmgx  (Figure  5*9)* 


FR  = 


(T  -  R2) 


s»n  a, 


01 


AR  =  Rj  +  R2  +  FR 


ZRQA  =  it  -  arc  sin 


AR  sin  a 


01 


Rj  +  T 


Z  QAR  =  it  “  L  RQA  -  a. 


01 


OR 


L  QAG  =  arc  sin 


OR  sin  a? 

[liipnf 


a.  and  L  QAG  >  0, 


O' 


AR  sin  0Q 

Per,metw  =  »'!nCW  -  «01  -"V 


2 

Area  =  0Q  R^  -  (Perimeter)  AR  sin  a01 

b.  Or  Z  QAG  <  0^ 

PQ  -  (Rj  +  T)  (90  -  L  QAG) 

Perimeter  =  PQ  +  QR 

=  0oRf2  -  (0O  -  ZQAGJCRj+T)2^.  AR  sin  0^ 


Area 


Slotted-Cone  (Continued) 


•5.1,2 


1>.  If  R2  <  T  <  Temp76  and  If  <  T  <  (Figure  5.10), 


ZASB 

Z  SBA 
Area  j 


arc  sin 


(W'W>stn  92 


(R5  +  T) 


*  -  &2  -  ZASB 

«0R?2  +  (Rf-R5"TW)(R5+T)s!n  *-SBA 


a.  And  ZRAG  >  6Q 


Z  BRA 
Z  RBA 


arc  sin 


(Rf-R^Tw)stn(a24»0) 


rt  -  ZBRA  "  ®2  “  ®0 


Z  R6S  =  Z  SBA  -  /.  RBA 


RS  =  ZRBS(R5  +  T) 

Perimeter  =  RS 

Area  =  Areai-(Rf-R5-Tw)  (R5+T)sln  ZRBA  -  ZRBS(R5+t)2 


b.  Or  Z  RAG  <  eQ 

Area2  *  Areaj  ~Z  RBS  (Rg+T)  2-  (R^-t^)  (R5+t)  s>:f)x  ABC 

(1)  And  If  ZQAG  >  8Q 

Z  QAR  =  eQ  -  Z  RAG 
L  AQR  -  x/2  +  Z  ARB  -  ZQAR 


OR  * 


AH  sin  ZQAR 

~7TirZm 


?er  I  *ter 


OR  +  RS 


Slotted-Cone  (Continued) 


Area  =  Area2  ~  OR  AR  sln(jt/2  -Z_ARB) 

(2)  Or  If  Z_QAG  <  8q 

Z_PAQ  =  80  -  Z_QAG 

PQ  =  Z.  PAtKRj  +  T) 

(Rj  +  t)  sin  ZlRAQ 
**  =  s!n(x/2  -Z.  ARB) 

Perimeter  =  PQ.  +  QR  +  RS 

Area  =  Area2-QR  AR  sln(*/2  -Z.  ARB)  -  ZlPAO.ftj+T)2 


5.  If  R,  <  T  <  Temp.,*,  and  If  <  T  >  T 


76' 

Z.  SAB  =  arc  cos 
Z.  SAG  =  L.  SAB  -  0 


7max 


(Figure  5*11)» 


«,2+(«f-»5-Tw)2  '  <Rr,+T>2  1 


AR  = 


[(R5+t)2  +  (Rf-R5-Tw)2-2(R5+T)(Rf-R5-Tw)cos  Z.ABc] 


1/2 


zLRAB  =  arc  sin 

L  RAG  =  Z_  RAB  - 

L-  ARB  =  *  -  Z.  ABC  -  L  RAB 
L  ARQ  =  £  -Z.ARB 


(Rc+x)  sin  Z  ABC 

_ 2 _ . _ _ _ 

AR 


Z.  AQR  =  *  -  arc  s 
L.  QAR  -  *  -  L  AQR  -  L  ARO. 


.  Tar  s  l  n  Z.  ARO 

foj  +  t) 


L.  QAG  =  Z.  QAR  +  L  RAG 


Slotted-Cone  (Continued) 


a.  And  SAG  >  6Q 

Complete  burnout  has  occurred. 

Perimeter  =  0.0 
Area  =  0.0 

b.  Or  SAG  <  9q 

4  ASB  ,  .rc  s.n  **-•*)] 

SBA  =  JC  -  e2  -  Z.SAG  -L  ASB 
Areaj  =  (0Q  -Z_SAG)Rf1 2 3  +  Rf (R^-Rg-T^)sln(92+  ZSAG) 

L  RQA  =  jc  -  arc  sin 

L  mi  =  |  +Z  ARB  -ZRQA 
L  QAG  *  Z  RAQ  ZRAG 
L-  RBS  =  Z.  SBA  -Zabc 
RS  =  Z-RBS(R5  +  t) 

(1)  and  If  Z  RAG  >  eQ  go  to  4a 

(2)  or  If  L. RAG  <  8Q  <  ZiJAG  go  to  4.bp) 

(3)  or  If  Z  RAG  <  3Q  >  ZQAG  go  to  4b(2) 

6.  If  R2  <  x  >  Tempyg,  and  If  Tfo  « 


Complete  burnout  has  occurred. 

Peclrseter  =  0.0 
Area  =0.0 


5.1.2  .Slotted-Cone  (Continued) 


7.  If  R,  <  t  >  TEMP.,*,  and  if  Tu*‘  >  T,<7',  and  t  <  t  (Figure 


5.12), 


r?6,  a[,u  Iy6  ^  ,v? 
>2 


f(Rv,  +  *0*  +  (Rf  -  Re  -  o2  “  (R5  +  T)2' 
■ are  cos  [— — — ; 


ZAQB  =  arc  sin 


r»f  -  rs  •  v 

L — 


ZABQ,  =  *  -  ZTaQB  -  ZlQAB 


ZLAS8  =  arc  sin 


"(R*  -  Rc  -  T 
'  f  5  w 

_  (v^ 


sln(ZQAB) 


)  sin(82) 


zab3  =  it  -  e2  -  Z'  ASB 
-Z.QAG  =  ZQAB  -  6. 


a.  and  Z'QAG  >  6Q 


Zf  AQB  =  arc  sin 


(Rf  -  Rc  -  \)  s>n(e2  +  e0)' 

- rt7R 7+H - 


S  A  AA 

4.  «ov< 


=  ?r  -  S2  -  8q  *•  ZAQ.B 


ZlQBS  =  ZlABS  -ZlABU 
PER1M  =  ^QBS(R5  T) 


b.  or  ZQAG  <  8- 


PERiM  =  (ZABS  -  ZABQ)(R_  *  T)  +  (0  -ZlQAGMR,  +  T) 


SHZk 

73,: 


S lotted-Cone  {Cor<: ! nued) 


8*  If  R,  <  7  >  Temp7fi,  and  <f  and  T  >  Tw  (FJgur*  5.12) 


'2  *  ' 

Z  SAB  *  arc  cos 


rRf2  +  («,-R.-tw)2-(r5+t)2'1 

~«7VVV 


L  SAG  =  Z_SA8  -  8. 


a.  and  ZSAG  >  6g 


Complete  burnout  has  occurred. 

Perimeter  =  0.0 
Area  =  0.0 


b.  Or  Z  SAG  <  eQ 

Z  QAB  =  arc  cos 
Z  QAG  =  Z  QAB  -  0 

Z  ABS  -  arc  cos 


(R^T) 2 4 (Rf-R  -Tw)2- (RS4T) 


2-(R  ~'2^ 
WfV  (Rf-R5-Tw) 


(Rf«R$>Tw)24(R^T)2-Rf2 


(1)  and  If  ZQAG  >  0Q 


Z  AQB  =  arc  s!r 


!(Rf-R5-T5(R>T) 


(R*’R5*Tw)aln(#24a0)l 
(rT+^-vT - 1 


Z  ABd  *  *  -  e2  -  80  -  Z  AQB 
Z  dBS  =  Z  ABS  -  Z  ABQ 
QS  =  Z  C.B$(R5  +  T) 

Perimeter  =  OS 


»  (8 2+»0  -  Z  $AB}Rf2*Rf  (Rf“R5-Tw)  *  In  Z  SAB 
- <R5+t)  2  Z  ass-  (Rf-R5- V  (R5+t)  *  l  n  z  ABd 


Area 


’  «**  '****••  *■*»«<  *<  >fca. .taut* 


Slotted-Cone  (Continual) 


■325 iK 
1 


I 

1 


i. 


(2)  Or  If  L  Q.AG  <  ®Q 

L  ABQ,  =  arc  cos 

L  OBS  »  L  ABS  -  L  ABQl 
Q3  =  L  OBS  (R5  -f  T) 

L  PAii  =  02  +  60  -  L  QAB 
PQ.  »  L  PAQ.  (Rj  +  T) 

Perimeter  =  PO  +  QS 

Area  =  (&2+«0  -  L  SAB)Af2+Af (Rf-K5“Tw)sin  L  SAB 

-(R5+t)2  ZQ11S-(Rj+T)(Rf-R5-.Tw)sIn  L  QAB 
-(e2+e0-  L  qab)(Rj+t)2 


(Rf  “R5"Tw)2+  fR/T)2^  (R1+T) 2  “ 


3.  If  R2  -<*  T  >  Temp76,  and  If  Tyg  <  T^?  (Figure  5.13), 

AT  =  j^(Rf-R5-Tw)  2+  (R5+T)  2~2  (Rf-R5-Tw)  (R^T)  cos  L  ABC  J  1/2 


L  ATB  =  arc  sin 


(R--RC~T  )aln  L  ASC* 
— 1,-2— £ _ _ _ 


AT 


L  STA  =  */2  -  L  ATB 
L  AST  »  x  -  arc  sin 


'AT  sin 

L  Rf 


L  STA 


L  TAB  =  v.  -  L  ABC  «  L  ATB 
L  SAT  =  ic  -  L  AST  «  L  STA 
L  SAG  =  L  SAT  +  L  TAB  -  Q? 

*.  a.  And  L  SA5  >  8. 

«M  ,1 


Complete  burnout  has  occurred 


5.1.2 


Slotted-Cone  (Continued) 


Perimeter  «  0,0 
Area  =  0.0 


b.  Or  L  SAG  <  6, 


L  QSA  =  it  -  Z  AST 


L  AOS  -  jt  -  sin 


R^  5 In  L.  QSnl 


(R.-S-T) 


L  QAS  «  *  -  L  AQS  -  /_  QSA 
ilO^G  =  4  SAG  +  L  QAS 


(1)  And  If  L  QAG  >  9Q 

L  QAS  «=  eQ  -  £  SAG 
L  AQS  =  rf  -  L  QAS  -  L  QSA 


Q$  - 


Rf  sfn  /.  QAS 
Vfn  ~'l  AQS 


Perimeter  =  Q$ 

Are*  =  4  QAS  *f2  -  Rf  QS  sin  L  QSA 

(2)  Or  If  L  QAG  <  (5„ 


L  PAQ  =  9q  -  L  QAS 
PQ  -  (Rj  +  t)  L  PAQ 


CIS  *  R-  - 

f  sin  L  AQS 

Perimeter  =  PQ  *  QS 

Area  =  (  L  PAQ  +  L  QAS)  Rf2  »  QS  Rf  sin  £  QSA 
-  .£  PAQ  (Rj  +  T)2 


»*/ 


5.2  Grain  Longitudinal  Geometry 

Representative  long! tudlnal  configurations  which  can  be  accc-mo- 
dated  by  the  computer  program  have  been  shown  In  Figure  2.2. 

Many  configurations  are  possible:  monolithic  or  segmented  grains 
In  the  cylindrical  section  of  the  motor;  head-end  with  web  or  a 
straight  through  grain  In  the  fore-head  section;  straight  through 
grain  In  the  aft-head;  external  and  Internal  taper;  and  elliptical 
contour  fore-head  and  aft-head  motor  sections  are  examples. 

Grain  cross  sectional  geometry  can  be  varied  within  the  cylindrical 
section  the  motor. 

5.2.1  Head-Fnd  wt th  Web,  Fore-Head  Section 

This  section  describes  the  analysis  of  the  head-end  with  web 
which  Is  solved  In  the  third  core  load  of  the  computer  program. 
Description  of  the  head-end  with  web  calculations  are  based  on 
the  motor  fore-head  section  shown  In  Figure  5*14.  Grain  geo¬ 
metry  within  the  fore-head  section  Is  based  on  the  Input  geome¬ 
trical  cross  section  at  the  forward  tangent  plane. 

The  burning  surface  area  for  all  distances  burned  and  the  tnttlal 
propellant  volume  are  obtained  by  separating  the  analysis  Into 
four  blocks.  Each  block  performs  the  following  function: 

Block  1  calculates  the  surface  area  versus  distance  burned  and 
the  InJtlal  volume  of  the  propellant  tip.  The  volume  and  area 
are  obtained  by  Integrating  elemental  areas  and  volumes  consist¬ 
ing  of  trapezoids  and  parallelograms. 

Block  2A  calculates  the  surface  area  of  the  pseudoelllpsotd 
minus  the  Igniter  opening  using  the  Theorem  of  Pappus.  Block 
23  calculates  the  volume  of  propellant  between  the  Inner  and  outer 
ellipsoids  and  the  surface  area  on  the  pseudo  ellipsoid  which  Is 
covered  by  the  propellant  tip.  The  volume  and  area  are  obtained, 
by  Integrating  elemental  areas  and  volumes  using  the  Theorem  of 
Pappus. 

Block  3  calculates  the  Initial  total  propellant  volume.  The 
volume  between  the  Inner  and  outer  ellipsoids  Is  determined  from 
the  difference  In  yolume  of  two  oblate  spheroids.  The  volume  of 
the  tip  u  obtained  by  adding  the  volumes  of  cylindrical  segment 
elements. 

Total  surface  area,  A,  versus  distance  burned  Is  obtained  by  com- 
blnlng  the  blocks: 

A  =  ABlock  1  +  A8iock  2A  “  ABlock  2B 

The  Initial  propellant  volume  Is  calculated  In  two  ways,  one 
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5.2.1  Head-End  with  Web,  Fore-Head  Section  (Continued) 


method  as  In  Slock  3  and  the  other  by: 

Volume  =  vBiock  i  +  VBlock  2B 

The  Initial  volumes  are  compared  In  subroutine  VOLSUB  and  a  corrcc 
tlon  factor,  AR,  Is  obtained  which  Is  added  to  the  total  burning 

surface  area.  This  Is  due  to  the  assumption  that  the  difference 
is  caused  by  the  Integration  method  tn  the  Block  1  analysis. 

Subroutine  HDHSUB 

Subroutine  HDNSUB  Is  the  control  routine  which  sets  up  the  cor¬ 
rect  variables  and  equations  to  perform  the  block  1,  2A,  2B,  and 
3  analysis. 

5.2. l.l  Block  1  Analysis 

The  generalized  forked  wagon  wheel  grain  configuration  Is  divi¬ 
ded  Into  13  sectors  (1,  2,  3,  3A,  3B,  4,  5,  6,  7,  9,  10,  11, 

11A,  11B,  12,  13)  as  shown  In  Figure  5.15.  Sectors  (1,  2,  4, 

5,  6,  7,  9,  10,  12  and  13)  share  a  common  analysis,  as  do  sec¬ 
tors  (3  and  11)  and  (3A,  3B,  11A,  and  11B);  however,  some  of 
the  sectors  have  special  equations  for  line  segments  and  angles. 
The  sector  boundaries  used  for  the  head-end  with  web  analysis. 
Figure  5*15,  are  different  than  sector  boundaries  used  for  the 
straight  through  grain  analysis.  Section  5.2.3  and  Figure  5*32  . 

oubr^utlne  SCI  Is  the  control  routine  to  determine  the  surface 
area  versus  distance  burned  and  initial  volume  of  the  propel¬ 
lant  tip  for  the  block  1  analysis  and  Is  called  from  subroutine 
HDNSUB  with  an  argument  L  to  indicate  the  sector  to  be  computed. 
The  analysis  for  sectors  1  and  3  will  be  explained  in  detail 
since  the  analysis  for  the  other  sectors  Is  similar. 

Figure  5.16  shows  2*ometrlc  constructions  and  calculation  con¬ 
trol  planes  used  to  obtain  elemental  burn  surface  areas  and 
propellant  volumes  tn  the  propellant  tip.  Also  shown  Is  part 
of  a  forked  wagon  wheel  grain  and  the  required  constructions. 

For  clarity,  only  one  tip  Is  shown.  Two  planes,  "A"  end  "B", 
control  volume  and  ares  calculations  for  the  propellant  tip. 

These  two  planes  are  generated  as  follows:  The  Intersection 
of  the  surface  of  the  propellant  tip  with  the  forward  tangent 
plane,  l.e.,  the  tangent  plane  perimeter.  Figure  5.16b,  Is  di¬ 
vided  into  Incremental  lengths,  &0.  Lines,  spaced  ALO  apart, 
perpendicular  to  the  tangent  pl.ne,  are  run  from  the  tangent 
plane  along  the  side  of  the  propellant  tip  to  Intersect  the 
Inner  ellipsoid.  These  Intersections  locate  points  P  and  P  . . 


Block  1  Analysts  (Continued) 


Perpendiculars  to  the  side  of  the  propellant  tip  from  points  PQa 
and  Po!j  locate  points  P^a  and  Plb  on  the  plane  of  s  /mmetry  or  on 

the  outer  ellipsoid  depending  on  the  initial  location  of  each 
At?  along  the  propellant  tip  perimeter  in  the  tangent  plane. 
Figure  5.i$a  shows  P^a  arid  P^  on  the  plane  of  symmetry.  Per¬ 
pendiculars  to  the  Inner  ellipsoid  at  points  PQa  and  PQb  locate 
points  P,a  and  P?b  on  the  outer  ellipsoid.  Points  P^,  Pla,  P2a, 
PQb,  P^b,  and  P2fa  define  the  A  and  B  planes.  Points  and 
lie  In  the  A  and  B  planes  and  are  located  on  the  outer  ellipsoid. 

The  trapezoid  formed  by  ALO  on  the  perimeter  of  the  tip  at  the 
tangent  plane  and  points  P  and  P  .  Is  the  area  calculation 
element. 

Perpendiculars,  In  "x-z:i  planes,  for  these  same  four  points 
form  the,  volume  calculation  elements  which  are  bounded  by  the 
traoezotdal  area  elements  as  shown  in  Figure  5 • 16.  The  ALO 
spacing  on  the  propellant  tip  perimeter  which  governs  plane 
placement  varies  In  value  along  the  perimeter  as  follows: 

The  ALO  spacing.  Figure  5.16b,  Is  calculated  first  as: 

a)  ALG.  of  each  sector  Is: 


m 

where  r—  Is  an  Input  parameter. 
Rf 


b)  ALO  Is  then  modified  by; 

,  n  J 

a-0  "  "(o-  .  [S~y~ 

'  pr  ps' 


where 

0  and  0  are  the  distances  between  planes  A  and  B, 
pr  ps 

Figure  5.17. 


The  starting  value  of  ALO  previous  Is  ALOj. 


5.2..1.1  Slock  1  Analysis  (Continued) 


The  factor  K,  an  Input  variable,  Is  the  major  parameter  In 
determining  the  distance  between  the  two  planes.  The  Incremen¬ 
tal  distance  ALO  Is  stepped  along  the  sector  perimeters  from 

L»  0  to  L=  L  . 
x  x  p 

5. 2. 1.1.1  Sector  1  (Figure  5.15) 

Subroutines  ASU8C,  BSUBC,  RASUBB,  XRSUB3,  THETAR,  GAM  SUB,.  GAMA2S, 
P0SU8,  P1SUB,  P3SUB,  ROPSB,  and  VSTRSB  are  used  to  obtain  the 
surface  area  and  Initial  volume  for  sector  1. 

The  distance  D  between  the  points  PQ  and  P^  on  the  Inner  el¬ 
lipsoid  are  shewn  tn  Figure  5.18  for  planes  that  are  located  at 
Increments  of  one-tenth  the  perimeter  length  of  a  sector.  A 
minimum  perimeter- length,  HOLDR,  Is  set  equal  to  the  perlme  er 
length  In  which  the  distance  for  the  plane  Is  greater  tnan 

the  burn  distance,  T.  The  first  plane  Is  located  at  the  peri¬ 
meter  length  HOLDR  and  subsequent  planes  are  spaced  the  distance 
ALO  apart.  H0L0R  Is  Initialized  to  the  perimeter  length  for 
plane  B  on  subsequent  Iterations. 

After  the  minimum  perimeter  length,  HOLDR,  Is  determined,  the 
distance  for  a  plane  located  at  the  total  perimeter  AL(1) 

for  the  sector  Is  determined.  The  parameter  TDMAX  hi  Initia¬ 
lized  to  the  maximum  value  of  TDMAX  or  for  the  sector  and 

then  checked  with  the  maximum  permissible  burn  distance,  TAUM. 

If  TDMAX  equals  or  exceeds  TAUM,  an  error  condition  exists  and 
execution  Is  terminated.  If  all  Dqj's  are  less  than  TAUM,  the 

radial  burning  surface  area  between  the  planes  and  the  burning 
surface  area  on  the  propellant  tip  are  determined. 

Radial  burning  surface  area.  Figure  5.17: 

as  i  =  lb**  V|(V  *  V?. 


Propellant  tip  burning  surface  area: 


AS  l 


(Lp  -  HOLDR)  (R2  -  T)  (Y(y,  +  YQg) 


whore  Yq  Is  defined  In  Figure  5.21 


5. 2. 1.1.1  Sector  1  (Figure  5.15)  (Continued) 


tp  Is  the  perimeter  length  along  the  Initial  grain  perimeter 
p  measured  from  the  beginning  of  a  sector  to  a  general  point. 
HOLOR  Is  the  prior  Incremental  value  of  Lp, 

The  Initial  volume  Is  determined  from  the  following  equation: 

<V0A  +  W  (lP  ‘  H0LDK)R2 
VSTR  *  - 5 - ! - 

Subroutine  ASUBC 

Subroutine  ASUBC  sets  up  the  correct  variables  and  equations  to 
det^rmlnr  the  coordinates  (X.  Y,  and  Z)  of  the  points  PQa,  P^a, 

and  P,  for  planes  located  at  Increments  of  one-tenth  L  along 
3a  p 

the  perimeter  of  the  sector.  Subroutines  RASUBS,  XRSUBB,  THETAR, 

GAHSUB,  GAMA2S,  POSUB,  P1SUB,  P3SUB,  and  TRAN  are  called  to 

determine  the  coordinates. 

Subroutine  RASUBB 

Subroutine  RASUBB  determines  the  length  of  the  radius  vector,  R  _, 

from  the  motor  axis  to  a  general  point  In  a  sector.  The  perimeter 
length  from  the  beginning  of  the  sector  to  a  general  point  Is 
required  In  calculating  a  value  of  RgT  for  each  sector.  Figure  5*19. 

Subroutine  XRSUBB 

Subroutine  XRSUBB  determines  the  X-coordlnate,  X  ,  for  a  point 

IO 

located  on  the  perimeter  of  a  sector.  The  parameter  RgT  from  the 

RASUBB  subroutine  Is  required  to  obtain  the  coordinate.  A  se¬ 
parate  equation  Is  required  for  each  sector.  Figure  5.15. 

Subroutine  THETAR 

Subroutine  THETAR  determines  the  angle,  9r,  between  the  Z-axls 
and  the  line  segment  R  The  parameters  X  and  R  _  are  re- 

o  i  r3  d  I 

qulred  to  obtain  0r,  Figure  5.15. 

Subroutine  GAHSUB 

Subroutine  GAHSUB  determines  the  angle  7j  between  the  line  normal 

to  the  perimeter  In  the  XZ-plsne,  and  a  line  normal  to  the  line 
segment  R#T,  Figure  5.19.  The  perimeter  length  of  the  sector  and 

the  angle  @r  are  required  to  obtain  7j.  A  separate  equation  for 

each  of  the  16  Individual  sectors  is  required.  7j  for  sec to/  3A 

Is  not  equal  to  7.  for  sector  3B  (also  true  for  sectors  11A  and 
11B).  1 


5.2, 1.1.1  Sector  1  (Figure  5. 1-5)  (Continued) 


| 


Subroutine  GAMA2S 

Subroutine  GAM2S  determines  the  angle  yn  between  the  Y-axis  and 


a  line  normal  to  the  ellipse 


-•  1,  which  Is  de- 


(ij 

,  at  the  point 

^2  I*  In  the  rotated  YZ'-plane  and  not  In  the  YZ-plane,  because 


fined  by  the  ellipse  ratio  P  ,  at  the  point  Z*  =  RAT.  Although 

Oc 


of  the  symmetry  of  the  ellipsoid  about  the  Y-axls,  computation 
can  occur  In  the  YZ-plane,  Figure  5.20. 


Subroutine  POSUB 

Subroutine  POSUB  determines  the  coordinates  X,  Y,  and  Z  of  the 
point  Pq  which  Is  located  at  the  Intersection  of  the  Inner  el¬ 
lipse  and  a  line  extended  along  the  propellant  surface  parallel 
to  the  tip  from  a  general  point  on  the  perimeter,  Figures  5.16 
and  5*21. 

The  parameters  RgT  and  0^  are  required  to  determine  the  coordi¬ 
nates.  The  coordinates  XQ  and  ZQ  are  obtained  by  trigonometry 
from  Rgl-  and  0r,  Figure  5*21.  The  equation  of  the  ellipsoid: 


yields  the  coordinate  Yp. 

The  remaining  points  defining  the  plane  (Pj,  P^j  and  P^)  are  ob¬ 
tained  from  Pq.  Radial  burning  In  the  plane  originates  at  PQ. 

Subroutine  P1SUB 

Su£>routlne  PlSUS  determines  the  coordinates  (X^,  Yj,  and  Zj)  of 
the  point  which  is  located  on  the  Y-Z  plane  or  on  the  outer 
ellipse  along  a  line  through  point  PQ  and  normal  to  the  sector 
perimeter  as  shown  In  Figure  5.22.  P^  Is  normally  located  on 
the  Y-Z  plane;  however,  for  sectors  7  and  9,  when  the  angle 
exceeds  In  sector  7  or  Pg^max  In  sector  9,  the  point  Pj 

Is  located  on  the  outer  ellipse  and  is  coincident  with  P y  The 
line  segment  ZlaT  shown  In  Figure  5.19  Is  used  to  flag  subroutine 
P3SUB  that  P,  Is  coincident  with  P,  whenever  Z,__  Is  non-zero. 


Sector  1  (Figure  5*15)  (Continued) 


P91max  and  P71max 
Section  5*1*1* 


are  geometry  plane  constants  discussed 


in 


Subroutine  P3SUB 

Subroutine  P3SUB  determines  the  coordinates  (X^,  Y y  and  Z^)  of 

the  point  that  Is  located  on  the  outer  ellipse.  If  Zja-jv  de~ 

termlned  In  subroutine  P1SUB,  Is  not  equal  to  zero,  P^  Is  colncl 

dent  with  Pj  and  the  coordinates  of  P^  are  set  equal  to  the  coor 

dlnates  of  P^.  If  2jaT  Is  equal  to  zero,  then  P^  Is  In  the  Y-Z 

plane  and  the  coordinates  of  P,  are  determined  from  the  angle  4> 
shown  In  Figure  5*22.  ^ 

Subroutine  ROPSB 


Subroutine  ROPSB  sums  the  values  of  the  Y  axis  coordinate  of  P 

(VQA  and  V  °r  V  (W  and  V0A  °r  W  and  Y0B}  l"  planeS  A 

and  B  which  are  used  In  subroutine  SCI  to  compute  the  surface 
area  along  the  grain  face.  If  the  burn  distance  x  Is  greater 
than  the  maximum  permissible  burr,  distance  for  the  sector,  then 
the  value  of  the  sum  Is  set  to  zero. 

Subroutine  BSUBC 


Subroutine  BSUBC  sets  up  the  correct  variables  and  equations  to 
determine  the  coordinates  (X,  Y,  tfnd  Z)  of  the  points  PQb,  Plfa, 

and  located  at  Increments  of  ztf.0  along  the  perimeter  of  sec¬ 
tor  1,  Figure  5.15*  The  A  plane  coordinates  ere  Initialized 
to  the  prior  B  plane  coordinates  during  the  Integration  along 
the  perimeter. 

Subroutine  VSTRSB 


Subroutine  VSTRSB  determines  the  Initial  volume  of  the  sectors. 
The  Incremental  volume  for  each  sector  Is  determined  from  the 
Incremental  cross-sectional  areas  of  the  sectors  and  the  average 
of  the  PQ  point  Y-coordJnatss  of  the  A  and  8  planes. 

Sector  3 

Subroutines  HASU8C,  HBSUBC,  HAPSBC,  HBPSBC,  STUPRS,  STUPPS, 
OPRASB,  RASUBB,  XRSUBB,  THETAR,  GAMSUB,  GAHA2S,  POSUB,  P1SUB 
P3SUB,  and  ROPSB  are  used  to  obtain  the  surface  area  and  initial 


"'-f  ■  ’  *  ** 
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5.2. \. 1.2  Sector  3  (Continued) 


volume  for  sectors  3,  3A,  and  3B.  As  l \  sector  1,  two  planes 
are  used  to  obtain  the  surface  area.  The  distance  between  the 
planes  A  and  8  Is  determined  by  ALQ  In  the  same  manner  as  for 
sector  1.  The  analysis  for  the  surface  area  In  sector  3  Is 
similar  to  the  analysis  for  the  surface  area  In  sector  1  ex¬ 
cept  that  the  A  and  B  planes  will  cross  as  shown  In  Figure  5.23 
for  sectors  3  and  11  and  In  Figure  5.24  for  sectors  3A,  3B, 

11A,  and  11B. 


The  perimeter  length  of  sector  3A  (AL3A)  Is  determined  fr^  the 
geometry  constants;  nnd  the  coordinates  of  points  PQ,  Pj,  and 

Since  the  surface  area  of  sector 


P^  for  both  planes  A  and  B. 

3B  Is  Identical  with  the  surface  area  of  sector  3A,  the  grain 
surface  areas  of  both  sectors  between  the  planes  are  determined 
simultaneously  by  stepping  the  B  plane  an  Increment  aLO  along 
the  sector  5A  perimeter.  The  surface  area  along  the  grain  be¬ 
tween  the  planes  Is  determined  from  the  following  equation: 


AS  I 


(*„  +  V) 


[<*o- V)2+  <zo- V)2i 


1/2 


The  surface  area  along  the  tip  for  sector  3  Is  determined  from 
the  following  equation  after  the  Intsgratloi  for  sector  3A  Is 
complete: 


AS  I 


(R2  -  T)  (Y^  +  Yob) 


whore  L,  Is  the  initial  sector  3  perimeter  length. 

The  Iv<Ittal  volume  of  sector  3;  and  sectors  3A  and  38,  Is  de¬ 
termined  at  the  beginning  of  subroutine  SCTQR1  after  completion 
of  the  block  1  analysis  as  follows: 


'STR 


^OA*Y3A't~YOB*Y3B^  fooi"0^  ^Y0B*Y3B*R3 


Subroutines  RASUB?,  XRSUBB,  THETAR,  GAMSUB,  CAHA2S,  POSUB,  P1SUB, 
PSSUB,  and  ROPSB  have  been  explained  In  Section  5*2. 1.1.1. for 
sector  1.  Therefore,  only  the  subroutines  unique  to  sector  3 
are  explained  In  this  section. 
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5. 2. 1.1.2  Sector  3  (Continued) 


Subroutines  HASUBC,  HBSUBC,  HAPSUBC,  and  HBPSUBC 

The  function  of  these  subroutines  Is  the  same  as  the  function  of 
subroutines  ASUBC  and  BSUBC  for  sector  1;  the  correct  variable 
and  equations  are  set  up  to  determine  the  coordinates  (X,  Y, 
and  Z)  of  the  points  PQ,  Pj,  and  for  planes  located  at  In¬ 
crements  of  ALO  along  the  perimeter  of  the  sector.  Subroutine 
HASUBC  deals  with  plane  A,  subroutine  HBSUBC  deals  with  plane  B, 
subroutine  HAPSBC  deals  with  plane  A*,  and  subroutine  HBPSBC 
deals  with  plane  8',  Figure  5.24. 

Subroutines  STUPRS  and  STUPPS 

These  subroutines  store  the  variables  that  define  the  planes  pro¬ 
duced  In  sector  3A.  Subroutine  STUPRS  Is  called  In  subroutine 
SCI  after  the  A  and  B  plnnes  have  been  defined  by  subroutines 
HASUBC  and  HBSUBC.  Subroutine  STUPPS  Is  called  In  subroutine 
SCI  after  the  A*  and  B1  planes  have  been  defined  by  subroutines 
HAPSBC  and  HBPSdC.  The  stored  variables  are  then  used  In  sub¬ 
routine  SCI  to  determine  the  area  along  the  grain  face. 

Subroutine  DPRASB 

Subroutine  DPRASB  determines  the  distance  between  the  points 

P  and  P1  and  between  P  an^  P1  that  lie  on  the  planes  pro- 
ra  ra  sa  sa 

duced  In  sector  3A  as  shown  In  Figure  5*24.  These  distances 

DP  ana  DP  are  usod  in  subioutlre  SCI  tc  alter  the  distance 
ra  sa 

between  the  planes  (£LQ)  for  each  successive  Iteration  as  ex¬ 
plained  In  Section  5*2. 1.1  for  sector  1. 


5-2. 1.2  Block  2A  Analysis 

Subroutine  SCT0R1  Is  the  control  routine  to  determine  the  sur¬ 
face  area  of  the  pseudoellipsoid  for  the  block  2A  analysis  and 
Is  called  from  subroutine  HDNSUB.  The  Initial  volume  ~f  , 
pellant  for  sectors  3A,  3B,  11A,  and  1  IB  Is  determined  pr  ..  to 
the  Integration  for  the  pseudoelllpsold  surface  area. 

Points  and  P,  that  are  used  In  Block  2A  are  not  the  same 
o  1 

plane  points  used  In  the  Block  1  analysis.  Points  PQ  and  Pj 

lie  on  the  Inner  ellipsoid  and  are  spaced  a  distance  DS  along 

the  ellipsoid.  Figure  5*25#  starting  at  the  Igniter  opening. 

The  spacing  of  P  and  P,  are  determined  as  follows:  J 

oi  j 


w 
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5. 2.1. 2 


Bloc  k  2A  Analysis  (Continued) 


First 


where 

r-*~  Is  an  Input  parameter 
Rf 

R^  grain  radius  at  the  fon  id  tangent  plane. 
Then  each  successive  Increment  Is 

■  f  « 


where  Is  an  Input. 

Initially,  If  the  motor  has  an  Igniter  hole  In  the  head  end, 

Zpg,  which  Is  the  Z-coordlnate  of  the  point,  PQ,  Is  set  equal 

to  R Otherwise,  Zpg  Is  set  equal  to  zero.  AR  Is  then  added 

to  Zro  to  obtain  ZP1'  which  Is  the  Z  coordinate  of  the  point  PI. 

Subroutine  YPSUB  Is  then  called  to  obtain  the  Y  coordinate  of 
points  PQ  and  Pj.  The  radius  of  curvature  at  the  point  Pj 

on  the  pseudoelllpsold,  the  length,  b,  of  the  Y-!ntercept  of  the 
normal  line  to  the  ellipse  at  PQ,  and  the  Z  coordinate  Zj  are 

then  determined  to  obtain  the  arc  length,  Lq.  Using  these 

values,  the  Incremental  strip  surface  area  is  determined.  This 
procedure  Is  then  repeated  by  setting  Z^  =  Zpl,  and  Zpl  =  Zpg 

+  £R  until  Zp^  =  R^  -  T^.  The  area  adjacent  to  the  Igniter  hole 

is  then  determined  and  added  to  the  Incremental  sum. 

Subroutines  YPSUB,  P.0E1SB,  LBSUB,  ZISU8,  and  AIGSU8  are  used  to 
detemlr.e  the  surface  area  of  the  pseudoelllpsold.  The  surface 
area  for  an  Incremental  strip,  as  shown  In  Figure  5.25,  Is  detcr» 
mined  by  the  fol lowing  equation: 


All 


■-S  (7 

2 NO  *4P1  + 


s,n  “rel  *Zn*T  5,n  °rcO> 
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Block  2A  Analysis  (Continued) 


The  surface  area  for  the  entire  pseudoelllpsold,  Including  the 
Igniter  hole,  Is  obtained  by  summing  the  Incremental  areas. 

The  surface  area  around  tbs  Igniter  hole  Is  added  to  the  sum  of 
the  Incremental  areas.  The  lateral  surface  area  of  the  Igniter 
hole  Is  assumed  to  be  a  non-burning  Inhibited  surface. 

Subroutine  YPSU5 


Subroutine  YPSUB  determines  the  Y  coordinate  of  a  point  on  the 
pseudoel Unsold.  The  angles  a  ^  and  Of  ^  of  the  points  P0a 

and  PQb  as  shown  In  Figure  5.25  are  also  determined. 


Subroutine  RQE1SB 


Subroutine  R0E1SB  determines  the  radius  of  curvature  p,  at  the 
po’nt  PQ  on  the  Inner  ellipsoid.  The  radius  of  curvature  Is 
determined  by  the  standard  equation: 


pi 


)».] 3/2 


where 


Z* 


dZ 

dy  * 


Subroutine  LBSUB 


Subroutine  LBSUB  determines  the  length  b  of  the  Y-Intercept  of 
the  line  normal  to  the  Inner  ellipse  at  PQ,  Figure  5.25-  The 

length  Is  used  In  subroutine  ZISUB  to  obtain  the  coefficients 
of  the  ellipse  equation  which  defines  the  pseudoelllpsold. 

Subroutine  ZISUB 


Subroutine  ZISUB  determines  the  Z-c(>ordInate  produced  by  the 
Intersection  of  the  outer  ellipse: 


and  the  normal  line  to  the  ellipse  at  PI. 
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5. 2.1, 2 


Slock  2A  Analysts  (Cortlnued) 


Subroutine  AIGSUB 

Subroutine  A?S>UB  determines  the  surface  area  around  the  Igniter 
opening.  The  area  Is  obtained  by  finding  the  surface  area  of 
revolution. which  requires  an  angle  it/NO#  an  arm  length 
t  Sln  a 

. .  •  2 . -  +  and  an  arc  length  T  a.^#  Figure  5*26. 


Initially#  the  Y-lntercept  of  line  normal  to  the  Inner  ellipse 
at  the  Igniter  radius  Is  obtained# 


Y 


NO 


B  2  _1LS_ 

S0E  2 

P0E 


1/2 


from  which  the  angle  between  the  Y-axis  and  the  Igniter  opening 
on  the  Inner  ellipse  Is  obtained# 


and  then  the  surface  area  of  revolution  Is  determined  (Theorem 
of  Pappus) 


f  a 


Alg  ‘  <»lg  *  T  sl!>  “rN>  TT 


rN 


5.2. 1*3  Block  2B  Analysis 

Subroutine  SCT0R2  Is  the  control  routine  to  determine  the  sur¬ 
face  srea  on  the  pseudoell Jpsold  that  Is  covered  by  the  pro¬ 
pellant  tips  for  the  block  2B  analysis.  The  surface  area  Is 
determined  by  dividing  the  cross-sectional  grain  configuration 
1 v to  12  sectors  In  subroutine  HDNSUB,  as  shewn  In  Figure  5*27. 
The  surface  area  on  the  pseudoelU  ps&Id  "vs*  each  sector  Is  then 
determined  by  subroutine  S2SK.  Summing  the  surface  areas  for 
all  sectors#  and  subtracting  this  area  from  the  area  obtained 
In  block*  1  and  2A,  yields  the  total  surface  area  of  the  head- 
end  with  web  for  any  value  of  t. 


31ock  2B  Analysis 


Subroutine  S2SK 


Subroutine  S2SK  determines  the  sectot  urface  area  on  the 
pseudoelllpsold  and  Is  called  from  subroutine  SCT0R2  with  an 
argument  to  Indicate  the  current  sector.  Each  sector  Is  set 
up  to  determine  the  Incremental  surface  area  contributed  by  that 
sector.  Thus,  In  sector  3,  ZpQ  Is  set  equal  to  Rp2,  then  Rp^ 

Is  obtained  arJ  9  Is  set  equal  to  9  ZjR  < s  then  added  to 
ro  r! 

ZpQ  to  obtain  Zpj  as  shown  In  Figure  £1.28.  YpQ,  Ypl  and  DS  are 

then  determined.  The  spactrig  for  £R  Is  calculated  In  the  same 
manner  as  was  done  In  Block  2A . 


Finally,  the  parameters  Lpp,  XrJ,  9r0,  9rJ,  a.c(J,  a^,  b, 
Zp  and  are  determined  and  the  incremental  surface  area  for 
the  sector  Is  obtained  from  the  equation.  Figure  5.25: 


-  fci 


and 


<erl  *  6r0> 


AS  -  +  T  arcl  +  2M  +  r  ,l„  arc0)  - 5 


Zpf,  Is  then  set  to  Zpj  and  the  procedure  repeated  until  Zpj 
equals  maximum  radius  of  the  sector. 


The  Initial  vo  isme  between  the  Inner  and  outer  ellipsoids  for 
each  sector  Is  determined  as  follows  (see  Figure  5.25): 


VJT0  '  h  +  V'"  ’i']  <8ri  +  erO><arcl  '  arcO> 
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wncre: 
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s,n  arcl 

pl^arcl  "  arcfP 

<P1  +  °WB)(arcl  -  arcQ> 
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Block  2B  Analysis  (Continued; 


5.2. 1.3 


(2DS2  *  PS1-1  °MB 
\»  ‘  3(HS2  +  #S1) 


(2P0  +  z»i> 


+  Zpl>  a„i  +  a,0 

_21_  +  1<w9  ^ 


Subroutines  YPSUB,  XRTHR,  uOElSB,  LBSUB,  and  ZIS'»B  are  used  to 
obtain  the  surface  area.  All  of  these  subroutines  except  XRTHR 
have  been  explained  in  Section  5*2. 1.2  for  the  block  2A  analysis. 

Subroutine  XRTHR 

Subroutine  XRTitt  is  n  set-up  subroutine  that  uses  subroutine 
XRSUiSB  to  obtain  the  X-coordinate  of  a  point  located  on  the  peri¬ 
meter  of  a  sector  shown  <r.  Figure  5.28.  The  angle  8rj  between 

the  Z  axis  ?nd  a  line  from  the  motor  axis  to  the  point  on  the 
sector  ts  also  determined. 


5- 2. 1.4 


Slock  3  Analysis 


Subroutine  VOLSUB  Is  the  control  routine  which  determir  s  the 
Initial  propellant  volume  for  the  block  3  analysis,  initially, 
the  volume  present  In  the  web  region  Is  determined  from  the 
volume  produced  by  the  difference  of  volunws  of  two  oblate 
spheroids  minus  the  volume  of  the  Igniter  hole  as  follows: 


'EH 


2 

3 


<>..  A.«' 


-  A')-  *  R 


oe  oe 


is 


<Bi.  - 


Bo«> 


Hext,  the  volume  of  the  propellant  tips  is  determined  by  dividing 
the  grain  configuration  into  12  sectors  as  was  done  In  block  2B, 
Figure  5.27*  The  volume  contributed  by  each  sector  (subroutine 
VSEC)  Is  then  added  to  to  obtain  the  total  volume  of  propel¬ 
lant  in  the  head-end,  V^.  Finally,  the  surface  area  generated 

by  radial  burning  from  the  line  of  Intersection  of  the  propel¬ 
lant  tip  with  the  head-end  wob,  A_,  is  approximated  in  the  fol¬ 
lowing  manner: 

The  propellant  volume  in  the  radial  burning  portion  is  deter¬ 
mined. 


VR  =  VfH  "  VEH  +  VSTO  "  2  VSTS  N0 
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5.2. 1.4 


Block  3  Analysis  (Continued) 


The  value  for  thickness  burned  at  which  the  maximum  surface  area 
of  this  radial  burning  portion  occurs  Is  assumed  to  be  £ha  thick"* 
r.ess  at  which  the  Initial  burnout  of  the  web  portion  of  the  heed- 
end  web  begins  (for  most  cases  this  Is  very  nearly  correct). 

The  curve  for  surface  area  versus  thickness  burned  for  this  fuel 
volume  Is  assumed  to  have  the  following  equations  (Figure  5.29); 


A.  = 


(TDMAX  -  tWH)q  T 

—  -  -  -  — 

WH 

(TDMAX  -  'T)q 


for  T  <  T, 


WH 


for  twh  <  T  <  TDMAX 


for  T  >  TDMAX 


V/nere 

TDMAX 

f 

WH 


maximum  burn  distance  In  the  head-end  section 
maximum  web  thickness  tn  the  head-end  section 
calculated  exponent 


The  maximum  thickness  burned  In  the  head-end  web  is  determined, 

lf  bie  -  boe  >  V  ’wh  -  \ 


lf  B!E  -  B0E  2  V  -  B!C  '  B 


OE 


The  radial  burning  volume  Is  matched  to  the  following  Integral 
*VR  =  VRX'  by  an  Iteration  process  to  determine  the  exponent  q, 

nm  (TDMAX  -  ly  )q  TrMAX 

VRX  =  J  - TdT  *  J  (TDMAX  “  T'  dT 

n  WH 

0  *WH 


or 


5.2. 1.4  Block  3  Analysis  (Continued) 

(TOMAX  *  twh)^  twh  (TOHAX  - 
VRX  =  1  “  +  q+I 

The  above  Integra  represents  the  area  beneath  the  curve  shown  In 
Figure 

The  radial  burning  area  AR  Is  added  to  the  total  surface  area, 

AS!,  at  the  completion  of  the  block  3  analysis  tn  subroutine  VOLSUB. 

Subroutine  VSEC 

Subroutine  VSEC  determines  the  Initial  volume  of  the  tips  that  are 
present  in  the  fore-head  and  Is  called  from  subroutine  VOLSUB  with 
an  argument  to  Indicate  the  sector  to  be  computed.  The  same  subrou¬ 
tines  used  In  the  block  2B  analysis  to  obtain  the  Incremental  surface 
area  of  a  pseudoellipsoid  are  used  In  this  subroutine  to  obtain  the 
Incremental  sector  volumes.  An  arc  length  LQ  Is  obtained  from  the 
expression  ^ 

Ld  =  ZP1  “  ZP0 

and  the  Incremental  volume  within  a  sector  Is  obtained  from  the 
expression 

m  =  so  ua(zp,  +  zpoi  !«.|  ♦  er0) (ypo  +  ypi>"1 


5.2.2  Cylindrical  Section 

The  Icngltudlnal  cylindrical  section  Is  that  portion  of  the  motor  be¬ 
tween  the  forward  and  *>ft  tangent  planes.  Figure  4.3.  its  length  Is 
designated  as  hco  and  the  radius  as  Rf  as  shown  in  Figure  5.30. 

It  may  be  either  straight  or  tapered  and  may  contain  either  a  mono¬ 
lithic  or  segmented  grain  with  or  without  v  tapered  poi't.- 

The  longitudinal  cylindrical  section  Is  first  divided  into  Incremen¬ 
tal  mass  addition  regions  by  locating  Increment  dividing  planes  every 
delta  Z  distance,  starting  at  the  forward  tangent  plane  and  proceed¬ 
ing  to  the  aft  tangent  plane.  A  maximum  of  100  Increment  dividing 
planes  is  allowed. 

The  longitudinal  cylindrical  section  is  then  divided  by  a  number  of 
reference  planes.  A  minimum  of  two  Is  required,  one  at  each  of 
the  tangent  planes.  A  maximum  of  eleven,  A  through  K,  Is  allowed 
to  indicate  changes  In  the  cross-sect lonal  grain  geometry  at  spe¬ 
cified  locations  within  the  section  as  shown  In  Figure  4.3.  The 
cross  sectional  grain  geometry  Is  described  at  each  reference 
plane  used,  by  the  port  perimeter,  Lp,  the  case  radius,  the 

port  area,  Ap,  the  propellant  area,  A^p,  the  ra-Uus  of  gyration, 

Kgy,  and  the  distance  to  the  Inert  silver,  T^^,  If  the  Inert 


Cylindrical  Section  (Continued) 
silver  option  l*  used. 

In  segmenting  propellant  grains,  the  location  of  the  slots,  which 
separates  the  grain  segments,  Is  determined  by  defining  the  dis¬ 
tance  of  the  slot  forward  and  aft  Interfaces  from  the  forward 
tangent  plane  of  the  motor.  This  Is  shown  by  tha  values  SH,  SIB, 
S2A,  S2B,  etc..  Figure  4.3.  The  Increment  dividing  planes.  In  this 
case,  arki  determined  as  above/  but  with  each  segment  treated  as 
having  a  forward  and  aft  tangent  plane.  The  reference  planes  nviy 
be  located  within  a  slot,  within  a  grain  segment,  or  on  a  slot 
Interface.  The  restrictions  on  the  number  of  Increment  dividing 
planes  and  reference  planes  Is  also  100  ar.d  11  respectively. 

The  following  describes  the  subroutines  which  determine  tha  longi¬ 
tudinal  geometry  for  the  cylindrical  sections 

Subroutine  HNCHN4 


Subroutine  MNCHN4  contains  the  program  control  logic  required  to 
obtain  the  internal  ballistic  solution  and  the  control  logic  used 
to  Initialize  .ihe  working  reference  planes  for  subroutine  SSGSUB. 
During  burning,  the  distance  burned  at  each  '•eference  plane  for 
each  time  Increment  Is  determined  by  linear  Interpolation  between 
adjacent  increment  dividing  planes.  When  the  reference  plane  Is 
located  within  a  slot,  the  Increment  dividing  planes  lossted  within 
the  grain  segment  are  used  for  extrapolation  to  obtain  the  reference 
plane  distance  burned.  When : the  distance  burned  has  been  determined 
for  each  Input  reference  'lane,  the  perimeter  length,  port  area,  and 
radius  of  gyration  for  e*  h  reference  plane,  and  the  burn  area  and 
CG  location  for  each  head  section  Is  determined  from  a  table  look¬ 
up  procedure  In  tha  geometry  tables. 

Ourlng  the  Internal  ballistic  solution  for  each  time  point  the 
cylindrical  section  working  reference  planes  (X  and  Y)  are  set  up 
for  successive  Input  reference  planes  (A-B,  B-C,  etc.). 

Subroutine  SEGSUB 


Subroutine  SEGSUB  contains  the  program  control  logic  which  determines 
the  perimeter  length,  cross  sectional  propellant  and  port  area, 
propellant  volume,  and  mass  generation  at  each  Increment  dividing 
plane  or  mass  addition  region.  The  control  logic  Is  set  to  check 
for  the  existence  of  a  slot  forward  Interface  between  adjacent 
Increment  dividing  planes,  to  check  the  location  of  the  current 
upstream  (X)  and  downstream  (Y)  working  reference  planes,  and  to  sst 
the  parameters  required  In  performing  the  gas  dynamic  solution  for 
the  slots  and  mass  addition  regions.  Each  increment  dividing  plane 
case  radius,  perimetar  length,  sliver  radius,  fuel  area  and  port 
area  and  radius  of  gyration  are  obtained  by  linear  interpolation 
between  each  reference  plane. 
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I  5*2*3  Straight  Through  Grain,  Motor  End  Sections 

The  following  Is  applicable  to  motor  aft-head  or  fore-head  soctlons 
with  straight  through  grains*  A  straight  through  grain  Is  shown 
In.  Figure  5*31.  The  E  subscript  shown  Is  set  to  N  for  analysts  of 
the  aft-head  and  to  H  for  the  fore-head. 

Baste  geometry  constants  required  for  analysts  are  calculated  prtor 
to  determination  of  burn  surface  area  versus  distance  burned  and 
propellant  volume.  These  constants  are  calculated  In  subroutine 
ENDCSB  ani  are  shown  In  Figure  5*30* 

5*2*3. 1  Geometry  Constants 

Subroutine  ENDCSB 


Initially,  the  case  opening  radius  Is  determined  from  the  Input 
parameter  DEI: 


Then  the  angle  between  the  tangent  to  the  ellipse  section  at  the 
radius  and  the  motor  axis  Is  determined: 


»  arc  cos  R£l/  [(R^2  -  Rgj)  3*  +  R^  3 


1/2 


where  R^  and  are  Input  parameters. 

If  ct^  Is  less  than  or  equal  to  the  maximum  allowable  angle, 

^Emax*  d#f tned  bV  ,nPut>  th,n  °b£max  U  s,t  to  °ER'  RE2  set 
to  Rgj  and  h^j  Is  set  to  zero,  if,  however,  Is  greater  than 

°i>iw  th" 

co»  Rf  h _ 

“  ’  '“*<W>2  4  -  -  ^W*  +  1-0)1/2 

end 

sin 

hEl  “  *RE2  "  RE1^  coT 

Next,  the  length  of  the  head  elliptical  section  Is  determined: 

1/2 

hE2  “  lRf2  "  RE23  7  PE 
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5.2.3. 1 


Geometry  Constants,  Subrc;»i!ne  ENDCSB  (Continued) 
and  the  length  of  the  end  section  Is  computed! 

hE0  =  hEl  +  hE2* 

If  hgQ  Is  less  than  the  maximum  burn  distance  at  the  adjacent 
tangent  plane,  'tmax>  determined  In  Section  5.1,  Geometry  Constants, 
the  end  section  Is  lengthened  to  the  maximum  burned  distance: 

h_  »  t 
ER  max 

If,  however,  h£0  Is  greaser  than  or  equal  to  t  ,  then  =  hEQ. 

Finally,  the  maximum  burn  distance  tn  the  conical  section  and  the 
end  section  Is  determined: 

2  2  1/2 

TE1  -  E<RE2  *  RE1>  +  hEl] 

2  2  2  1/2 
TE0  =  EERf  “  RE1*  +  hE0^ 

The  complete  end  section  case  volume  Is  determined  from: 

R  2 

VCE  =  EhER  "  *W  *  Rf  +  3  E3  "  hE2^  E*  hE2PE^ 

+  (R£2  R£1  +  Re1  +  RE2)  *  y 

The  following  coefficients  are  calculated  for  use  In  subroutine 
RCSUB  when  the  burning  distance  T  >  They  are  the  coefficients 

of  a  fourth  degree  equation  of  the  Intersection  of  an  ellipse  and 
a  circle  and  are  used  to  compute  the  location  of  Intersection  of 
the  burning  surface  and  the  case  wall. 

CAE  -  (P£2  -  l)2 

CBE  =  -(eE2  -  1)  PE2  RE1  4 

CCCE  -  0E4  Re2  3  +  (PE4  +  PE2)  hE0  +  R2  (p2  -  1)  -  R2  P2  ]  2 
CCVE  -  2(PE4  -  pE2) 

COCE  *  [  (h20  +  R2j)  p2  *  Rf2  p2  ]  4Re1 
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Geometry  Constants,  Subroutine  ENDGS6  (Continued) 


COVE  =  4  R  f3E4 

CECE  =  Rf4  +  <Ref  -  h£2)  PE2  Rf2  2  +  (R-2  +  h£2)  p£4 
C«VE  =  (h£2  +  Re2  3e2  +  RfZ)  2  PE2 


5. 2. 3« 2  Calculation  Sectors  and  Zones 

The  cross-sect! on*. i  grain  geometry  In  motor  end  sections  Is  based 
on  the  geometry  a-:  the  forward  and  aft  tangent  planes.  Figures 
5.32  and  5*33  define  sectors  and  zones  used  for  burn-surface  area 
and  propellant  volume  calculations.  The  general  forked  wagon  wheel 
grain  In  Figure  5*32  Is  divided  Into  13  sectors.  Depending  on  the 
exact  grain  geometry,  more  than  one  sector  may  exist  In  a  zone. 

This  Is  the  case  for  £one  A,  Figure  5*33  and  Initial  sectors  1, 

2,  3,  and  4  In  Figure  5-*32>  The  opposite  may  exist  where  only  a 
portion  of  a  sector  Is  bounded  by  a  zone  boundary. 

Four  zones  can  exist  In  the  end  section.  Their  boundaries  are 
defined  as: 

Zone  A  Is  that  region  between  the  minimum  propellant  radius, 

Rj,  and  the  case  opening,  R£j.  Burning  will  occur  along  the 

side  of  the  propellant  tip  and  on  the  end  face. 

Zone  B  Is  that  region  between  the  case  opening,  R_j,  and 

where  the  end-burning  surface  Intersects  the  case,  R  .  Burn- 

c 

Ing  will  occur  along  the  propellant  tip  and  on  the  toroidal 
end  face. 

Zone  C  Is  that  region  between  where  the  end-burning  surface 

Intersects  the  case,  R  ,  and  the  web,  t  -  ,r.  Burning  will 

c  w 

occur  only  along  the  side  of  the  propellant  tip  or  In  the 
valley. 

Web  Zone  Is  that  region  between  the  radius  R,.  -  t  +  t  and 
•*  t  w 

the  radius  of  the  case  Burning  will  occur  along  the  peri¬ 

meter  of  sector  8  and  along  the  end  face  tf  the  case  opening, 
R^j,  Is  greater  than  the  radius  to  the  web,  R^  -  +  T. 

The  sectors  that  are  In  zones  B,  C,  and  web  are  subdivided  Into 
smaller  elements.  The  sectors  In  zone  A  are  not  subdivided.  Sec¬ 
tor  boundaries  and  zones  which  the  sectors  occupy  era  recalculated 
as  the  burning  surface  regresses. 
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Calculation  Sectors  and  Zones  (Continued) 


Subroutine  ASESUB 

Subroutine  ASESUB  Is  the  control  routine  to  determine  the  total 
burning  surface  area  and  Initial  volume  of  the  end  sections.  The 
grain  cross-section  at  the  adjacent  tangent  plane  Is  divided  Into 
sectors.  Figure  5*32.  The  correct  equations  from  subroutines 
XRSUB  and  RASUB  are  set  up  and  the  proper  values  for  the  coordlnat 
of  the  origin  of  the  circular  arc,  X^,.  the  radius  of  curva¬ 
ture  of  the  sector,  RT,  the  angle  between  the  bisector  of  the  pro¬ 
pellant  tip  and  the  straight  side  sectors,  0^,  and  the  perimeter 
length,  L^,  are  assigned  for  each  sector.  After  the  required 

parameters  have  been  set  for  a  sector,  the  sector  burning  area 
and  volume  are  determined  In  subroutine  AESUB. 

Subroutine  AEPSUB 


Subroutine  AEPSUB  tests  for  the  existence  of  a  sector  and  's 
called  from  subroutine  ASFSUB,  If  a  sector  has  burned  ou"  or 
was  not  present  in  the  Initial  grain  configuration,  control  Is 
returned  to  subroutine  ASESUB;  otherwise,  control  proceeds  to 
subroutine  AESUB.  After  the  sector  burning  area  and  volume  have 
been  determined  In  subroutine  AESUB,  the  sector  values  are  added 
to  the  sum  of  the  values  for  the  previous  sector  and  control  Is 
returned  to  subroutine  ASES',3. 

Subrout  I ne  AESUB 


Subroutine  AESUB  determines  the  total  burning  surface  area  and 
Initial  volume  of  the  sectors.  A  test  Is  made  at  the  beginning 
of  each  zone  to  determine  If  the  sector  exists  In  that  zone. 


The  parameters  XRQ,  RAQ, 


and  Yaq  are  associated  with  the  beginning 


coordinates  of  a  sector  perimeter,  and  the  parameters  X^,  RA^, 
and  are  associated  with  the  end  coordinates  of  a  sector 


perimeter  *5  shown  In  Figure  5.32.  These  parameters  are  deter¬ 
mined  by  subroutines  XRSUB  and  RASUB  from  the  Pythagorean  Theorem. 
The  beginning  and  end  points  of  a  sector  may  not  be  the  beginning 
and  end  points  of  the  area  to  be  computed  should  a  sector  exist 
In  more  than  one  zone.  The  beginning  coordinates  of  an  area  are 
then  subscripted  min  and  the  end  coordinates  are  subscripted  max 
Instead  of  0  and  X,  respectively.  The  perimeter  length  Is  then 
defined  as  L^.  Figure  5 shows  an  example  of  the  min  and  max 


coordinates  that  are  used  In  the  following  zone  A  analysts.  When 
an  Integration  scheme  ts  employed,  zones  B,  C,  and  Web,  the  min 
and  max  coordinates  of  each  Increment  are  determined  to  compute 
sector  surface  area. 
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5.2. 3*2  Calculation  Sectors  and  Zones  (Continued) 


Subroutine  RCSUB 

Subroutine  RCSUB  obtains  the  radius  vector  from  the  motor  axis 
to  the  intersection  of  the  aft  burning  surface  with  the  case  wall, 

Rc,  Figure  5*33.  The  Intersection  Is  determined  from  the  equation 

of  a  straight  line. when  t  <  te1>  and  from  the  equation  for  the 

aft-dome  configuration  (circle  or  ellipse)  when  t  >  V  where: 

TE1  =  ^RE2  "  ReP  “  hEl  ^ 

Subroutine  AR3SUB 

Subroutine  ARSSUB  determines  the  chord  length,  tRS>  between  the 

minimum  point  of  a  sector  and  a  general  point  along  the  perimeter 
of  a  sector  as  shown  In  Figure  5*35*  The  required  parameters 

are  RA0’  XR0’  V  a"d  V 

Subroutine  ALRSUB 

Subroutine  ALRSUB  determines  the  arc  length  of  a  sector,  LR,  from 

the  minimum  point  of  a  sector  to  a  general  point  along  the  peri¬ 
meter  as  shown  In  Figure  5-35-  The  required  oarfmeters  are  L,- 
and  Rr  ™ 

Subroutine  XRSU-'i 

Subroutine  XRSUB  determines  the  X-coordlnate  of  a  general  point 
on  the  perimeter  of  a  sector.  A  separate  equation  I  s  used  for 
each  sector.  The  required  parameters  are  R^  and  T, 

Subroutine  RASUB 

Subroutine  RASUB  determines  the  length  of  a  radius  vector  from 
the  motor  axis  to  a  general  point  on  the  perimeter  of  a  sector. 

A  separate  equation  Is  used  for  each  sector.  The  perimeter 
length  along  the  sector  and  tho  distance  burned  t  are  required. 

Subroutine  HESUB 

Subroutine  HESUB  determines  the  length  of  the  trapezoidal  elements, 
hg,  used  to  determine  the  Incremental  volumes  and  areas.  Figure  5-3*>. 
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I 

;*)«2,3.2.1  Zone  A  Calculations 


When  Raq  <  Rg^,  the  surface  area  Is  computed  in  Zone  A. 
>  Rgj,  the  analysis  proceeds  to  the  next  zcne  and  R^^ 
the  smaller  value  of  R^  and  hi' 


When  Raq 
Is  set  to 


The  following  Zone  A  analysts  applies  to  .ectors  1  through  7  and 
9  through  13.  Figure  5*35  Is  used  as  an  example  to  define  cel" 
culatlon  parameters.  It  Illustrates  sector  >. 


The  burning  surface  aree  In  zone  A  Is  determined  from  an  algebraic 
composition  of  simple  geometric  figures  such  as  shown  In  Figure 
5.35. 


The  value  of  y^,  LR,  and  are  determined  «*s  follows: 

i2  ,  2,1/2 


7j  -  2  arc  cos 


Si  -  h 


“RS 


(IfajY 


V 


>  ■  Wi"  \ 

2Rt  ) 


R  4  .  y  S 1/2  /r2  y  S 1/2  . 

R*— ■  X“—  }  {RAmIn  XRmIn  }  + 


Amax  ‘Rmax 

(X  “X  1  i/^ 

l*Rmax  *Rm.!n;  j 

where  Ry  =  R^  -  T  for  sector  3 


AT0  fs  the  area  between  the  chord  LRS  and  the  circular  a;c  LR 

and  Is  determined  by  subtracting  the  area  of  the  Inscribed  tri¬ 
angle  from  the  area  of  the  circular  sector,  abg: 


Aro  -  Irt'  lr  *  Sc?  I"/  '  j 


2,12 


A^j.  Is  uqual  In  magnitude  to  A^.q  and  Is  positive  If  R^.  Is  positive, 
negative  If  R  Is  negative,  and  zero  IF  R^  Is  zero.  A?p  Is  the 
area  of  the  trapezoid,  aboF: 


y2  ,1/2  r„2  Jl  ,1/2  |  *XRmax"XRmtn* 

ARmax*  "  l*Aroln  ”  RmIr,J  !  2 
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5.2.3.2<.l  Zone  A  Calculations  (Continued) 

AR  I n  the  area  bed  determined  by  subtracting  the  triangular  ares 
obd  from  the  circular  sector  obc: 

7  ..  x  fR2  -  x3 

.  wnax  7R  .  *Rmax  lRAm*x  xmaxJ 


where 


7r  =  arc  sin 


(  XRmax  \ 
\  RAmax  / 


Arq  Is  the  area  aef  determined  by  subtracting  the  triangular  area 
oaf  from  the  area  of  circular  factor  oae: 

R^  y  -  X  fR2  -  X2  iX/^ 

.  RAmln  7R0  *Rmln  LKAmIn  xmlnj 


where 


r„0  -  •«  s 


In  ) 

\RAmIn  / 


The  burning  area  along  the  side  of  the  propellant  tip  Is  h^CU) 
and  tne  total  sector  area  Acc  Is:  fe  R 

Aee  -  Z  HO  [LRhE  *  An  *  AFf  ♦  Ar  -  Aro] 

The  sector  volume  Is; 

DV  "  <AFF  +  AA  -  aro  +  V  hE  1(0  2 

If  the  maximum  pol.it  of  the  sector,  R^max  Is  within  zone  A,  then 

control  Is  returned  to  subroutine  AEPsUB;  otherwise,  computation 
will  proceed  to  zone  3. 


o.2. 3*2. 2  Zone  8  Calculations 


When  Rac  for  a  sector  Is  less  than  R^.,  surface  area  Is  computed  In 
zone  B.  When  R^Q  Is  greater  than  or  equal  to  R£,  computation 
proceeds  to  Zone  C. 

Initially,  R^m#x  Is  set  to  the  smallest  value  of  R^  and  Rc  from 
which  Lb__„  Is  determined  by  subroutine  ARSSUB  and  AIRSUR.  L,,* 


Zone  B  Calculations  (Continued) 


is  determined  In  the  same  manner  from  RAmjn*  The  burning  surface 
area  Is  then  computed  In  Increments  of  £L  where  AL  -  R^.  (DLRF) 
by  Integrating  over  the  perimeter  length,  LRmax»  DLRF  Is  an 
Input  parameter.  For  each  Increment,  RAmax  Is  determined  from 
L  =  +  AL,  and  R^j  Is  deteiii.'~"d  from  L  =  Lrj;,  Figure 

The  X-coordlnate,  XR,  of  the  centroid,  Figure  5.36;  Is  obtained 

from  subroutine  XRSUB  and  the  radius  vector  R^  Is  obtained  from 
subroutine  RASUB.  The  angle  7R  Is:  a 


arc  sin 


The  cross  sectional  area  of  the  end  face  Is  equal  to  the  product 

of  the  Incremental  arc  length  (\  -  \  ,  It,  and  the  arc  length, 

min 

w  through  which  the  centroid  Is  rotated. 

The  surface  area  along  the  side  of  the  propellant  tips  are 
approximated  by  trapezoidal  increments  and  Is  added  to  the  area 
of  the  end  face.  The  burning  surface  areas  and  volumes  of  each 
Increment  are  added  *-o  the  sum  of  the  previous  Increment  values. 
They  are  determine.*  > or  an  Increment  as  follows: 


EE 


^  2  NC 


hp  +  hi 

<L-W  -h-1  +  ra  >r  fc  -  WT 


where 

hr 


‘  1 f  '  (Ra™»  •  Rn)2l‘/2| 


'Amax  "El" 


4 


I^ER  ’  &  -  <RAmln  “ 


h  +  h« 

OV  =  2  NO  (R.  -  R.  ,  )  ~~ — -  R. 

'  Amax  Amin'  2  A 


The  next  increment  In  the  sector  Is  determined  by  setting  i 

hE  =  hE"  RAm!n  =  RAmax>  and  LRi  =  LR1  +  ^  =  L*  L.!s  set 
LRmax  for  the  last  Iteration  for  a  sector. 


R1 

to 


=  L 
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Zone  B  Calculations  (Continued) 


The  parameters  Rfi,  Rfl  ,  L,  and  h_  are  determined  for  each. In¬ 
crement.  *  Amax  & 


The  tofoldal  end  burning  arra  Is  formed  by  revolving  an  arc  about 
the  motor  axis,  and  (s  equal  to  the  product  of  the  length  of  the 
arc  increment  and  the  arc  length  through  which  the  centroid  of 
the  arc  Increment  Is  rotated  (Theorem  of  Pappus).  The  area  along 
the  side  of  the  propellant  dp  Is  determined  from  the  trapezoids. 


The  radial  vector,  RA,  to  the  centroid  of  the  arc  Increment,  as 
shown  In  Figure  5.36*  Is: 


where 


1% 


X  -  X 


dn 


nm 


Is  the  distance  from  the  centroid 


of  the  arc  Increment  to  the  origin  of  the  circle  about  which  the 
arc  Increment  Is  revolved. 


The  angles  X  and  X^  are  determined  as  follows: 


X  =  arc  s 


m  (SaasLil) 


\in  *  arc  sl 


If  R^  Is  less  thsn  R^,,  computation  will 
sector.  If  R^  Is  greater  than  or  equal 
proceed  to  zone  C. 


proceed  to  the  next 
to  RC,  computation  will 


Zone  C  Calculations 


The  burning  surface  area  In  Zone  C  consists  only  of  the  surface 
along  the  side  of  the  propellant  tip,  and  Is  determined  from 
trapezoidal  elements,  if  a  sector  exists  in  zones  B  and  C, 
is  set  equal  to  R^^  of  xone  B;  otherwise,  R^  *  RAQ.  With 


is  determined  from  subroutines  ARSSUB  and  ALRSUB.  The 


length  of  the  Initial  edge,  h£,  of  the  trapezoidal  element  Is 
determined  In  subroutine  HESUB  as  follows; 


SHEET 

ICS 


5.2. 3. 2. 3  Zone*  C  Calculations  (Continued) 


IfRA<RE2, 


hE  ■  hER  "  hEl 


/RA  ~  RE1  \ 
\R£2  ‘  RE1  / 


Or  If  R,  >  R„, 

(Rf 2  -  »2) 1/2 

hE  “  h£R  +  p6  nE0 


An  Incremental  length  &L  Is  added  to  and  a  corresponding  RA 

Is  determined  from  subroutine  RASUB.  The  length  of  the  top  edge 
h£  of  the  trapezoidal  element  Is  determined  from  the  above 

equations  In  subroutine  NESUB  and  the  elemental  trapezoidal  area 
and  volume  are  added  to  the  previous  sectors  as  follows: 


/h  +  h|\ 

Are  -  2  *  <L  *  W  (-V1) 


where 


»  RAmax  *  ^Amln 

‘A  s'  2 


(hF  +  h*) 

DV  =  2  NO  (R.  -R.  ,  )  — v 

'  Amax  Amin'  2  R 


5.2, 3. 2.^  Web  Zone  Calculations 

Subroutine  AWESUS  determines  burning  surface  area  and  Initial 
volume  of  the  web  zone. 

The  burning  surface  area  of  sector  8  Is  determined  first  from 

the  trapezoidal  element  of  length  hf>  which  Is  dett,  ..Ined  from 

R.rt  =  -  t ,  +  t.  as  follows: 

AO  f  W 

AWE  =  2  NO  h£  Lg 

Lg  Is  the,  perimeter  length  of  sector  8. 

The  volume  of  propellant  In  the  web  zone  Is  determined  as  follows 
(geometf  ic  symbols  are  given  In  Figures  5.33  end  5*36) 


5*2. 3*2.4  Wftb  Zone  Calculations  (Continued) 

If  >  (Rf  “  ca^culat«  the  parameter, 

(R„  -  Rf  +  T  -  T) 

h  a  -—ST. . —L. - Z  ■  .  .  h' 

VC  R£2  -  Re1  hE 


and  the  volume  as: 


0V  -  j[(^)  hE2  ’  ^-f-]  h2  *  <hER  -  hE0>  Rf2 

+  [<Rf  -  \  *  T)2  +  RE2<Rf  ’  Tw  +  T)  +  RE22] 
"  <bEFC  +  bER  ‘  hElJ <Rf  "  T„  +  T>2  * 

J 

If  re2  —  (Rf  "  V"  ca^cu^at®  the  parameter, 

K2-  <Rf- •„  +  t'2I  1/2 


zi  - 


and  the  volume  as: 


■  itef  ■.  •  vj  v  •  * 


ER  “  hElP  Rt 


“  *Z1  +  hER  "  "  Tw  +  T)2]  * 

The  additional  end  burning  area  and  Initial  propellant  volume, 
when  Rgj  >  Rf  -  4  xy  are  determined  as  follows: 


yhen  Re1  <  Rc 

'*«■  -  2  *>  81  <REl2  -  RA02' 


OV  -  *  «E  (re12  -  rA02) 


The  toro.'tel  end  burning  surface  exists  when  >  R»,.  The  re- 

V  t  * 


,1 


5.2. 3.2.4  W«b  Zone  Calculations  (Continued) 


volyed  area  Is  a  product  of  the  arc  length  and  the  circumference 
of  the  circle  described  by  the  centroid  of  the  arc  length.,  the 
radial  vector  to  the  centroid  of  the  arc  Increment,  Figure  5*36, 


9  _  _ 2T_ 

h  ~  \  -  \ 


The  toroidal  surface  area  Is  determined  from  the  product  of  the 
arc  length,  -  \mjn)T,  and  the  circumference  of  the  circle, 

2  NO  8j  R^,  described  by  the  centroid  of  the  curvet 
\E  -  -  W  T  “  *A  8! 


Monwits  of  Inertia  and  CQ  Location 


Pitch  and  roll  moments  of  Inertia  (HOI)  and  the  center  of  gravity 
(CQ)  location  during  motor  burning  can  be  calculated.  The  roll 
moment  (J-ROLL)  la  taken  about  the  longitudinal  axis  of  the  motor. 
The  pitch  moment  (J-X-Y)  Is  taken  about  an  axis  passing  through 
the  motor  CQ  and  centerline.  The  center  of  gravity  Is  measured 
from  the  eft  tangent  plane.  The  value  Is  positive  when  the  CQ  Is 
forward  of  the  aft  tangent  plane.  The  moments  about  the  pitch  and 
yaw  axes  are  assumed  to  be  equal.  This  assumption  Is  valid  for  any 
conflgtratlon  with  an  even  multiple  of  4  propellant  tips. 

The  *&l  and  CQ  location  of  the  motor  are  based  on  the  combined 
values  of  each  section^  fore-head,  cylinder,  and  aft-head. 


The  pitch  HOI  of  the  fore-head  Is  Initially  determined  about  the 
forward  tangent  plane  and  then  transferred  to  the  aft  tangent  plane, 
the  cylindrical  add  aft-head  section  HOI*s  are  Initially  determined 
about  the  aft  tangent  plane.  These  values  about  the  aft-tangent 
p’ane  are  then  transferred  to  the  motor  CG.  The  transfer  formula 
ins  , 

JCG  ~  ^aft  tangent  pla  e  +  ^ 


where  d  Is  the  distance  between  the  motor  CQ  and  eft  tangent  plane. 


The  HOI  of  a  body  with  respect  to  e  given  axis  Is  defined  as  the 
product  of  the  mass  and  the  square  of  the  distance  from  the  axis. 
If  dm  =  dW/go  represents  an  elemental  mass  end  Y  Its  distance  from 
an  axis,  the  HOI,  J,of  the  object  about  this  axis  will  be  equal  to 

//  dW/g,,,  Reference  8. 


The  CG  Is  that  point  at  which  the  mass  of  an  object  Is  concentrated 
so  that  the  moment  of  the  concentrated  mass  about  any  axis  or  plane 
Is  the  same  as  the  sum  of  the  moments  of  all  the  elements  of  the 
mass  about  the  same  axis  or  plane.  The  sun  cf  the  moments  from  a 

O1*™  15  /kdW 
J  So 

and  the  CG  ts  defined  as:  x  = 


where 


x  3  distance  to  the  CG  location  from  a  plane. 


The  HQI  end  CG  location  .or  each  section  ate  determined  from  a 
summation  of  Incramental  volumes  and  areas  about  the  desired  axis. 

Tha  Incremental  volumes  are  hollow  circular  cylinders  for  the 
cylindrical  section  end  thin  shells  fo~  the  heed-end  sections. 

The  radius  of  gyration  for  the  cyllndr  _sl  section  Is  determ. ned 
from  the  roll  HOI  of  the  grain  cross-sectional  area.  The  Incremental 
rectangular  HOI  s  are  first  taken  about  the  CG  of  the  lncrementt.1 
volume  and  then  transferred  to  the  desired  axis.  The  transfer 
formula  Is  I  «  5cG  +  where  ?cG  the  HOI  about  the  CG  end 
d  Is  tha  distance  to  the  reference  axis. 


5.3.1 


5.3.1. 1 


Roll  Moffcent  of  Inertia 

The  roll  HOI  of  an  Incremental  thin  shell  for  the  heed  sections 
ere  determined  from  the  fundamental  equation;  Figure  5.37: 

dp  -  7-  r2,  slug- In2 
P  9  o 

The  thin  shells  are  based  on  a  subdivision  of  sector  boundaries. 

The  roll  HOI  of  the  cross  sectional  area  of  on  Incremental  hollow 
circular  cylinder.  Figure  5*3§>  Is  determined  from  the  fundamental 
aquation: 

Jp  -  <f.4  -  '  I4)  HO,  In4 

where 

rQ  =  outer  Incremental  radius.  In 
r(  =  Inner  Incremental  radius.  In 
The  radtus  of  gyration  Is  determined  from: 

in 

K"  >  ,n 

2 

where  AFP  Is  the  cross  sectional  area  of  propellant.  In  . 

Motor  Eno  Sections 

The  roll  MO  I  s  for  the  motor  end  sections  are  determined  In  sub¬ 
routines  PT1AA  and  SD1D13  from  a  summation  of  thin  shells  of  each 
sector  and  th<?  wab  region  as  follows: 

t 

1;  9999  (hp+hp)  p# 

AJPP  -  ^  ,x]  ♦  (R20  ♦  Ra2)  (Ra2  -  RaJ>  1C 2^ 

where,  Jpx  ■  A  +  «A2} 

*So 

Ws  -  (9rI  ♦.  ®r3>  NO  pf  (R2  -  R  2)  <hE  ♦  h^) 

- % - 

Rao  Is  the  minimum  radius  of  the  Increment 


Ra  Is  the  maximum  radius  of  the  Increment 


5.3.1»2 


5.3.2 


Cylindrical  Suction 

The  roll  HOI  for  tha  cylindrical  suction  l«  determined  In  subrout  in* 
3EGSUB  using  the  rad’us  of  gyration  of  th«  propallant  cross  section. 
The  radius  of  g> ration  of  the  cross  ruction  is  determined  In  sub¬ 
routines  PT1AA  and  SD1013  from  a  s«;metIor>  of  the  roll  HOI  for  the 
elementary  hollow  circular  cylinders  of  each  sector  and  the  web 
region  as  follcwst 


Z  [  (ra4-ra^>  no  i.f  (O-r^) 

1 


WT  ‘  f  I  CA\oi  »  )  +  *  (*A2-»Ao> 

.  .i-r 

Nsy  "  *■  WT  3 

The  radius  of  gyration  is  calculated  orly  et  the  Input  reference 
planes  and  linearly  Interpolated  at  each  increment  dividing  plane. 

The  cylindrical  section  total  roll  MO?  Is  determined  from  a  summation 
of  individual  mass  addition  region  roll  M9I  values  as  follows. 

Figure  4.5* 

Nl 

!pcyl  -  2  t  («  »w?  -  APH|)  ^4  *  <*  r/  -  V  ]  4l£i 

The  subscripted  Hi  parameters  are  the  values  at  the  adjacent  upstream 
increment  dividing  pl-.ne.  HI  Is  the  total  number  of  Increment 
dividing  planes. 

Pitch  MO  I 

The  pitch  HOI  s  of  the  head  sections  are  taken  with  reference  to 
the  adjacent  tangent  plane  such  that: 

.  _  V.  fl!  +  jA  +  SL  P2 
JB  ~  go  (  2  +  12}  +  go  RC0 


where 


r  =  radius  to  CG  of  cross  section.  In 

L  =  Incremental  length  of  shell.  In 

Rrr  =  distance  from  reference  axes  (tangent  plena) 
to  the  CG,  In 
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5.^2 


Pitch  HOI  (Continued) 


For  the  Incremental  hollow  circular  cylinders  of  the  cylindrical 
section  the  calculations  «re  taken  with  reference  to  the  CG  of  the 
Incremental  cylinder  and  then  transferred  to  tin  aft  tangent  plane 
such  thats 

i  W  ,K2  .  lA  .  W  .2 
*  TT  (“?  +  +  7T  d 


where 


5.3.2. 1 


K  =  radius  of  gyration  of  cross  section.  In4 
d  «  distance  from  the  CG  to  ait  tangent  plane.  In 
Motor  End  Sections 


The  pitch  MO  I  for  the  motor  end  sections  are  determined  In  sub¬ 
routines  PT1AA  and  SD1D13  Worn  a  summation  of  M0|  s  for  elemental 
thin  shells  of  each  sector  and  the  web  region  as  follows: 

....  £  1  r  (hE';'hE'  ,2  Wl  JPX  Wl  Rca  1 
AJBB-i.  |t  2  12i^  +  — +  — | 


RC«WI 


(hE+hE) 


+  +  V  +  rao>  ’sfe 


where 


rcg  = 


hE0  + 


hg+hg 


(h£+hE) 


for  the  fore-head  and 


for  the  aft-head 


2  2  Wi 

JPX  =  *RA  +  RA0>  2 go 


5.3>2.2 


Cylindrical  Section 

The  MO  I  Is  determined  In  subroutine  SEGSUB  from  a  summation  of 
individual  jm*?  addition  region  values  as  fallows: 


5lff 

f  - 


—  1  "gWI  (AINCW-A6NCH0  +  .  A|NJ 


5.3. Jr, 2  Cylindrical  Section  (Continued) 


5.3.3 


5.3*3. 1 


£  i  ♦  i  **  ♦  2  (msqjmn  *  „co .  AWa,)*j 
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C6  of  Sections 


The  CQ  for  each  section  Is  determined  from  e  summation  of  the 
moments  of  Incremental  volumes  about  the  desired  axis.  The  CG 
for  the  cylindrical  section  Is  determined  from  a  summation  of  the 
moments  of  Incremental  hollow  circular  cylinders  about  the  aft 
tangent  plane.  Figure  5.3$  and  the  CG  for  the  head  sections  are 
determined  from  a  sanitation  of  the  moments  of  thin  shells  aboit  the 
adjacent  tangent  plane.  Figure  5>37.  Thus: 

MX  *  ~  X  «  Z  X. 

So  1  S0 

therefore 


x  -  rxf  dW/W 

where 

Xg  --  moment  arm  of  Incremental  volume.  In 
=  Incremental  volt  me  weight,  lb 
W  =  total  weioht,  lb 
Head  Sections 

The  CG  Is  determined  In  subroutines  PT1AA  and  SD'^IS  from  a  summation 
of  moments  for  elemental  thin  shells  of  each  senior  and  the  web 
regton  as  follows: 

X|  =  X  [  (WT  Xj  previous  +  *hE  +  M  T  ^Wl  +  WT^ 

13 

where  W-  *  r  W. 

T  1  1 

subscript  I  Is  N  for  aft-haad, 

H  for  fora-head. 


S.3.3.2 


Cylindrical  Section  CG  Location 


The  CG  location  Is  determined  In  subroutine  SEGSUB  from  a  summation 
of  the  Individual  mass  addition  region  moments  as  follows.  Figure 
5.38: 
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5.3.4 


Motor  HO i  end  CG 
Motor  Roll  MO  1 1 

JPP  3  JPHed  *  JPNoz  +  *Pcyl 

The  motor  CG  location  and  pltxh  MOi  ana  determined  by  a  transfer 
of  axes  of  the  MOi  of  each  section  as  follows: 


Motor  CG: 


<*H  +  W  vfH  "t  +  A0WCyl-  -  *N  VfK 


Fore-Head  Section  Pitch  MOI: 

JBHed  3  ^BHedWt  "  [  *H  lXH  +  hC0  "  *1^  3 


Aft-Heed  Section  Pitch  MOI: 

'Soz  "  'JBNozWt  ”  ^  *  H  "  +  XiH^  ^  ■  f^-Pf 

9o 

Cylindrical  Section  Pitch  MOI: 


‘Beyl  =  ^BcylUft 


!BcylWt  ”[(  ^  XIH  Wf  "  +  hC0*  VfH  pf 

f  *H  VfN  °f  Wf  “  ^VfH  +  VfN^  Pf  ■*  )  } 
[  Wf  -  (vfH  +  vfN)  Pf  3  /  So 


Motor  Pitch  MOI: 

J8B  3  JBHed  +  JBNoz  +  'ecyl 


FIflyr*  5.2  One-Half  Fork of  General  Modified  Wagon  Wheel  Configuration 
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HfKPt  5.29  DistritnrttaB  of  Volme  in  IteAUl  Burning  Section 
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5.33 


Figure  5#}8  Cylindrical  Suction  Elements!  Volumes  for 
HOI  and  CG  Calculations 


DETAILED  PROGRAMMING  INFORMATION 


Information  required  to  maintain  and  revise  the  program  Is 
presented  in  this  section.  A  brief  description  of  all  program 
subroutines,  macroscopic  program  logic  flow  charts,,  the  program 
storage  allocation,  a  description  of  the  computing  system, 
program  diagnostic  aids,  and  a  list  of  the  program  nomenclature 
are  included. 

As  stated  in  Section  3.0  on  the  method  of  solution,  the  primary 
purpose  of  this  computer  program  Is  to  obtain  the  performance 
characteristics  of  solid  propellant  rocket  motors,  which  requires 
the  solution  of  the  Internal  ballistics.  In  order  to  obtain 
the  performance  characteristics,  the  computer  solution  is 
separated  into  two  major  sections.  One  section  obtains  the 
propellant  burning  surface  area  and  volume  (geometry  calculations), 
and  the  other  section  obtains  the  solution  of  the  gas  dynamic  equa¬ 
tions  (internal  ballistic  solution).  The  computer  solution  of  the 
geometry  equations  requires  three  separate  core  loads  and  the 
computer  solution  of  the  gas  dynamic  equations  requires  a  fourth 
core  load.  Each  computer  core  load  contains  a  control  subroutine 
which  calls  the  required  subroutines  and  computes  values  of  control 
parameters.  A  common  data  region  named  C0M  and  dimensioned  6600, 
which  resides  in  the  machine  memory  core  at  all  times  during 
program  execution,  contains  all  of  the  computed  variables  so  that 
each  core  load  may  be  linked  to  all  other  core  loads.  Another 
common  section  named  TBLK  and  dimensioned  2580  is  u?«d  to  store  the 
geometry  tables  whHch  are  calculated  In  core  loads  1-3  and  which 
are  used  in  core  load  4.  The  Individual  core  load  control 
subroutines  are  linked  through  the  MAIN  program  by  a  control 
vai  fable  named  ICHN.  The  control  subroutine  for  the  first  computer 
core  load  is  MNCHNl ,  fer  the  second  coTputer  core  load  -  MNCHN2, 
for  the  third  computer  core  load  -  MNCHN3,  and  for  the  fourth 
computer  core  load  -  MNCHN4. 

The  first  computer  cor*  load  (ICHN  =  1)  contains  the  subroutines 
required  to  read  the  input  data,  initialize  the  dot*  coils, 
compute  the  input  reference  plane  constants,  locate  the  incre¬ 
ment  dividing  planes,  check  for  input  data  errors,  and  print 
the  program  inputs  and  computed  constants.  The  second  computer 
core  load  (ICHN  «  2)  contains  the  geometry  subroutines  required 
to  compute  the  tnitiel  propellant  area,  perimeter  length,  end 
radius  of  gyration  for  the  cylindrical  section  reference  planes; 
end  to  compute  the  Initial  propellent  volume,  burn  area,  center  of 
gravity  location,  end  moments  of  Inertlalfor  the  eft-heed  end 
straight  through  grain  fore-head  section.  The  third  computer  core 
load  (ICHN  -  3)  contains  the  geometry  subroutines  required  to 
compute  the  initial  propellant  volume,  burn  area,  center  of  gravity, 
end  moments  of  inertia  for  the  fore-heed  with  web.  The  fourth 
computer  core  load  (ICHN  -  4)  contains  the  subroutines  required 
to  obtain  the  internal  ballistic  solution. 


Subroutine  Description 

This  subsection  is  divided  Into  two  pcrts.  Part  one  (6.1.1) 
contains  a  verbal  description  of  the  purpose  of  each  subroutine 
and  Includes  cross-references  to  other  sections  which  give 
additional,  detailed  explanations  *f  particularly  complex  or 
important  subroutines.  Part  two  (b.1.2)  contains  a  tabular 
description  of  the  subroutine  linkage.  The  table  Identifies 
all  \&/<er  level  subroutines  called  by  a  particular  routine  plus 
all  routines  which  call  that  particular  subroutine. 

Subroutine  Descriptions 

ACOS  -  Subroutine  ACOS  determines  the  arc  cosine  for  argu¬ 
ments  between  0  and  2k  radians. 

AEPSUB  -  Subroutine  AEPSUB  tests  for  the  existence  of  a 

sector  for  the  end  sections.  If  a  sector  Is  burned 
out  or  does  not  exist,  program  control  is  returned 
to  subroutine  ASESUB,  otherwise  computation  proceeds 
to  subroutine  AESUB.  The  sector  areas  are  then 
summed  and  program  control  Is  returned  to  subroutine 
ASESUB  (Section  5. 2. 3. 2). 

AESUB  -  Subroutine  AESUB  determines  the  surface  area  and  Ini¬ 
tial  volume  of  a  sector  for  the  zones  In  the  end 
sections  (Section  5*2. 3. 2). 

AFPSUB  -  Subroutine  AFPSUB  determines  the  perimeter  length 
and  cross-sectional  area  of  all  sectors,  except 
sector  8,  In  the  cylindrical  section  (Section  5  • 1.1) • 

AIBST  -  Subroutine  AIBST  determines  the  pressure  loss  and 
changes  In  the  gas  properties  between  increment 
dividing  planes  due  to  mass  addition  and  area  change 
for  the  non-steady  flow  soluUon  of  the  Internal 
ballistics  (Sectton  4. 1.2.1). 

AIBSUB  -  Subroutine  AlBSUB  determines  the  pressure  loss  and 
changes  In  the  gas  properties  between  Increment  di¬ 
viding  planes  due  to  mass  addition  and  area  change 
far  the  solution  of  the  Internal  balli3ttcs  neglec¬ 
ting  transient  effects  (Section  4.1.1). 

AIBSUB  -  Subroutine  AiGSUB  determines  the  surface  area  around 
the  igniter  opening  in  the  head  end  section  (Section 
5.2. 1.2). 

-  Subroutine  ALRSUB  determines  the  arc  length  of  a 
sector  in  the  end  sections  from  the  minimum  point  of 
a  sector  to  a  general  point  along  the  perimeter  of 
the  s s'. tor  (Section  5. 2. 5.2). 


AUSUB 


6.1.1 


Subroutine  Descriptions  (Continued) 


ARSSUB 

ASESUB 

ASTSUB 

ASUBC 

AWESUB 

BRAKSB 

BSUBC 

COHPSC 

CCNV 


Subroutine  ARSSUB  determines  the  chord  length  between 
the  minimum  point  of  a  sector  and  a  general  point 
along  the  perimeter  of  the  sector  (Section  5*2. 3. 2). 

Subroutine  ASEhUB  sets  up  the  correct  equations  for 
subroutine:  XRSUB  and  RA3UB  and  assigns  for  each 
sector  the  proper  values  for  the  coordinates  of  the 
origin  of  the  circular  arc  (RAO,  XAO),  the  radius  of 
curvature  of  the  sector  (RT),  and  the  perimeter  length 
(AL)  of  the  sector  (Section  5 *2. 3. 2). 

Subroutine  ASTSUB  sets  up  the  correct  equations  for 
subroutine  PT1AA  to  determine  the  moments  of 
Inertia  for  the  block  1  analysis  of  the  head  end 
with  web  (AJSTP  and  AJSTB) . 

Subroutine  ASUBC  sets  up  the  c  Tect  variables  and 
equations  to  determine  the  coc  j I  nates  (X,  Y,  and  Z) 
of  the  points  PDA,  P1A,  and  P3A  for  the  block  1 
analysis  in  subroutine  SCI  (Section  5.2. 1.1). 

Subroutine  AWESUB  determines  the  total  burning  sur¬ 
face  area  of  the  web  zone  at  thickness  TAU  and  the 
burning  surface  area  of  sector  8  for  the  end  sections 
(Section  5.2. 3.2.4) . 

Subroutine  8RAKS8  determines  the  length  of  the  dia¬ 
gonal  of  the  parallelogram  that  Is  formed  by  the 
Intersection  of  two  planes  In  the  block  1  analysis 
(See  Figure  5.2’). 

Subroutine  BSUBC  sets  up  the  correct  variables  and 
equations  to  determine  the  coordinates  (X,  Y,  and  Z) 
of  the  points  POB,  P1B,  and  P3B  on  the  pseudoelllpsold 
for  the  block  1  analysis  In  subroutine  SCI  (Section 
5. 2. 1.1.1). 

Subroutine  COHPSC  determines  the  sine  and  cosine  for 
the  reference  piane  geometry  angles  A(l)  thru  A(5). 

Subroutine  CONV  determines  new  Iterations  value  for 
either  AKRST  or  PCTAB  during  the  start  transient 
calculations  for  the  internal  ballistic  solution. 
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■I  C.l.i 


Subroutine  Description-;  (Continued) 


CSTRSB 


GPRASB 


ENOCSB 


FOGRE 


GAHA2S 


GAHSUB 


Subroutine  CSTRSB  determines  the  propellart  gas  pro¬ 
perties,  '“STAR,  molecular  weight,  specific  heat 
ratio,  gas  constant,  and  combustion  temperature  as 
a  function  of  the  chamber  pressure  during  the  In¬ 
ternal  ballistic  solution  (Section  4.3.3)* 

Subroutine  OFRASB  determines  the  distance  between 

the  points  P  and  P*  ;  P  and  P*  ;  P  .  and  P'.J 
ra  ra  sa  sa  ro  rb 

and  P^  and  P*^  that  lie  on  the  planes  produced  In 

sectors  3A  ar.d  3B  or  11A  and  11B  In  the  block  1 
analysis  of  subroutine  SC!  (Section  5*2. 1.1. 2). 

Subroutine  ENOCSB  determines  the  coefficients  of 
a  fourth  degree  polynomial  equation  obtained  by  the 
Intersection  of  an  ellipse  and  a  circle,  CAE,  C8E, 
CCCE,  CCVE,  COCE,  COVE,  CECE,  and  CEVE  and  deter¬ 
mines  the  constants  RE1,  RE2,  ALFE,  ALFEM,  HE1, 

HE2,  and  HEO  that  define  the  geometry  of  the  fore¬ 
head  and  aft-head  sections  (Section  5*2. 3.1). 

Subroutine  FOGRE  solves  for  the  roots  of  a  fourth 
degree  polynomial  equation  by  reducing  to  the  form 

X4  +  A  X2  +  8X  +  C  =  0 


by  the  substitution 
Y  =  (X  -  P/R) 


and  then  solving  the  resultant  cubic 


T3  +  |  T2  +  ^  yg  4C^-  T 


Subroutine  GAMA2S  determines  the  angle  between  the 
Y-axfs  and  a  line  normal  to  the  ellipse 

(bIe*  +  te)  =  1 

(which  Is  defined  by  the  ellipse  ratio  pn_)  at  the 
point  Z‘  =  RAT  (Section  5*2. 1.1.1).  Ufc 

Subroutine  GAHSUB  determines  the  angle  between  a 

line  normal  to  the  perimeter  In  the  X-Z  plane,  and 
a  line  normal  to  the  line  segment  RAT  which  Is  a 
radial  vector  from  the  motor  axis  to  a  point  on  a 
sector  perimeter  (Section  5*2.1.1»1). 
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6*1.1 


Subroutine  Descriptions  (Continued) 


HAPSBC  -  Subroutine  HAPSBC  sets  ip  the  correct  equations  to 
determine  the  coordinates  (XOP,  YOP,  ZOP)  of  the 
points  PDA',  P1A',  and  P3A'  on  the  pseudoell Ipso  Id 
for  the  block  1  analysis  of  the  Intersecting  planes 
for  sectors  3  and  11  In  subroutine  SCI  (Section 
5*2.1. .1*2) . 

HASUBC  -  Subroutine  HASU8C  sets  up  the  correct  variables  to 
determine  the  coordinates  (XO,  YQ,  ZO)  of  the  points 
POA,  P1A,  and  P3A  that  lie  on  the  pseudoell Ipsold 
produced  In  the  block  1  analysis  cf  the  Intersecting 
plane  for  sectors  3  and  11  in  subroutine  SCI  (Section 
5. 2. 1.1. 2). 

HBPSBC  -  Subroutine  HBPS8C  stts  up  the  correct  equations  and 
variables  to  determine  the  coordinates  (XOP,  YOP, 

ZOP)  of  the  points  PCB',  P1B’,  and  P3B*  that  lie 
on  the  pseudoell Ipso  Id  produced  In  the  block  1 
analysis  of  the  Intersecting  planes  for  sectors  3 
ana  11  In  subroutine  SCI  (Section  5.2. 1,1. 2); 

HDNSUB  -  Subroutine  HDNSUB  sets  up  the  correct  variables  to 
perform  the  block  1,  2A,  2B,  and  3  analysis  of  the 
head  end  with  web  (Section  5*2.1). 

HESUB  -  Subroutine  HESUB  determines  the  length  of  the  zones 
In  the  end  sections  from  which  the  Incremental  sur¬ 
face  ar^as  In  subroutine  AESUB  are  determined  (Sec¬ 
tion  5*2. 3*2). 

LBSUB  -  Subroutine  LBSUB  determines  rhe  length  (b)  of  the  Y 
Intercept  produced  by  the  Intersection  of  the  line 
Y  =  0  and  the  line  normal  to  the  ellipse  at  PG 
(Section  5*2. 1.2) . 

LPDAPS  -  Subroutine  LPDAPS  sets  up  the  correct  variables  to 
determine  the  perimeter  length,  ALP,  and  the  cro^s- 
sectlonal  fuel  area,  AFP,  of  the  propellant  tips  !n 
subroutine  AFPSUB  for  the  cylindrical  section  re¬ 
ference  planes  (Section  5*1*1)* 

-  Subroutine  LPTO  determines  the  perimeter  length  AL7 
and  AL8  of  the  anisotropic  propellant  In  sectors  7 
and  8  during  the  motor  tall-off  Interval  (Section 
4.3*3). 


LPTO 


6.1.1 


Subroutine  Descriptions  (Continued) 


MNCHN1  -  Subroutine  MNCHN1  Is  the  control  routine  for  the 

first  core  load  which  reads  the  Input  data,  computes 
the  constants  that  define  the  geometry  of  the  grain 
cross-section  and  kngltudlnal  conf Ig^-atlon  for  use 
by  the  second  and  tHH  core  loads,  checks  for  data 
errors,  and  prints  the  program  Inputs  and  geometry 
constants. 

MNCHN2  -  Subroutine  MNLHN2  Is  the  control  routine  for  the  se¬ 
cond  core  load  which  determines  the  Initial  fuel 
area  (AFF)  and  port  area  (A')  for  the  cylindrical 
section  reference  planes  and  determines  the  burn 
area  and  Initial  fuel  volume  of  the  straight  through 
grain  fore-head  and  aft-head  sections  as  a  function 
of  distance  burned. 

MmCHN3  -  Subroutine  MNCHN)  Is  the  control  routine  for  the 
third  computer  core  load  which  determines  the  Ini¬ 
tial  fuel  volume  and  burn  area  as  a  function  of  time 
for  the  head  end  with  web. 

MNCHN4  -  Subroutine  MNCHN4  Is  the  control  routine  for  the 

fourth  computer  core  load  which  obtains  the  Internal 
balltstlc  solution  (Section  5.2.2). 

HODTSB  -  Subroutine  HODSTB  modifies  the  value  of  TIME  and 

the  value  of  the  nozzle  throat  diameter  after  con¬ 
vergence  of  the  ballistic  solution  Iteration  and 
Initializes  X  and  Y  reference  planes  for  Interpola¬ 
tion  of  the  Increment  dividing  planes  In  subroutine 
TISUB. 

MSISUB  -  Subroutine  MSISUB  determines  the  location  of  the 
center  of  gravity  and  the  polar  and  rectangular 
moment  of  Inertia  for  the  block  2A  5”"lyr.js  of  the 
head-end  with  web. 

MTfSUB  -  Subroutine  MTISUB  determines  the  location  of  the 
center  of  gravity  and  the  polar  and  rectangular 
moment  of  Inertia  for  the  block  2B  analysis  of  the 
head  end  with  web. 

PCHGJS  -  Subroutine  PCHtfU  prints  the  program  Inputs  for  the 
propellant  properties  and  the  burning  rate  equation 
constants. 

PCHPLI  -  Subroutine  PCHPLl  prints  the  geometry  constants  cal¬ 
culated  for  each  reference  pl^ne. 


Subroutine  Descriptions  (Continued) 


PLNCNS 

PLNt.CS 

POSUB 

PT1AA 

P1SUB 

P3SUB 

RASUB 

RA3UBB 

RBSTSB 

RBSUB 

I5BVSUB 


Subroutine  PLNCNS  calculates  the  geometry  constants 
of  each  cylindrical  section  reference  plane  (Section 
5.1.1). 


Subroutine  PLNLCS  checks  for  reference  plane  Input  data 
errors,  prints  appropriate  diagnosed  comments,  and 
flags  the  program  for  case  termination  If  an  orror  exists.- 

Subrouttne  POSUB  determines  the  coordinates  (XO,  YO, 

ZO)  of  the  point  PO  that  Is  located  on  the  Inter¬ 
section  of  the  Inner  ellipsoid  with  the  core  (Section 
5.2.1. 1.1). 


Subroutine  PT1AA  determine*  the  polar  and  rectangular 
moment  of  Inertia  for  the  cylindrical  section,  straight 
through  grain  end  sections,  and  the  block  1  analysts 
of  the  head  rnd  with  *.  sb  (Section  5.3). 

Subroutine  P1SUB  determines  the  coordinates  (XI,  Yi, 
and  Zl)  of  the  point  PI  that  Is  located  on  the  Z  a? Is 
along  a  line  through  potnt  PO  and  normal  to  the  sector 
perimeter  (Section  5. 2. Ul. 1). 

Subroutine  P3SUB  determines  the  coordinates  (X3,  Y3, 
and  Z3)  of  the  point  P3  ti»t  Is  located  In  the  Y-Z 
plane  and  on  the  outer  ellipse  (Section  5.2. 1.1.1). 

Subroutine  RASUB  determines  the  length  of  a  radius 
vector  from  the  motor  axis  to  a  point  on  the  peri¬ 
meter  for  each  sector  In  the  end  sections  with  a 
straight  through  web  (Section  5*2. 3*2). 

Subroutine  RASUBB  determines  the  length  of  a  radius 
vector  from  the  motor  axis  to  a  general  point  In  a 
sector  for  the  block  1  analysis  of  the  head  end  with 
web  (Section  5*2. 1.1.1). 

Subrouth>-  RBSTSB  determines  the  Initial  estimate  of 
the  burn  urs;e  coefficient,  AKRST,  for  each  time  In¬ 
crement  during  the  start  transient  Interval  end 
performs  the  table  look-ups  for  PH  and  AKKST  (Section 
4*3*3) . 

Subroutine  RBSUB  determines  the  propellant  burning 
rata  at  each  Increment  dividing  plane  (Sectton  4.3.2). 

Subroutine  RBYSUB  checks  the  validity  of  the  burning 
rate  equation  constants  and  prints  appropriate  diag¬ 
nostic  comments. 
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6.1.1 


Subroutine  Descriptions  (Continued) 


RCSUB  -  Subroutine  RCSUB  determine?  the  value  of  the  radial 
distance  from  the  motor  axis,  to  the  Intersection  of 
the  aft-end  burning  surface  and  the  motor  case  for 
any  configuration  (Section  5. 2. 3. 2). 

RGISUB  -  Subroutine  RGISUB  determines  the  radius  of  gyration 
for  the  Incremental  thin  shells  of  the  head  end  with 
web. 


R0E1SB  -  Subroutine  ROEISB  determines  the  radius  of  curvature, 
p^,  at  the  point  on  the  pseudoelllpsoia  for  the 

block  2A  analysis  In  subroutine  SCI  (Section  5-2. 1.2). 

R0P5B  -  Subroutine  R0P5B  sums  the  values  of  the  Y  coordinates 
PC  (YOA  and  YOB)  or  P0l  (YOA'  and  YOB')  for  the  A 
and  B  planes  from  which  the  block  1  surface  area  Is 
obtained  In  subroutine  SCI  (Section  5*2. 1.1-1). 

RSSPLN  -  Subroutine  RSSPLN  determines  the  coefficients  of 
the  piecewise  cublcs  that  are  used  f r, •  the  spline 
Interpolation  procedure. 

SC i  -  Subroutine  SCI  Is  the  control  routine  to  determine 
the  surface  area  and  Initial  volume  of  the  propel¬ 
lant  tip  for  the  block  1  analysis  of  the  head  end 
with  web  (Section  5*2. 1.1). 

SCT0R1  -  Subroutine  SCT0R1  Is  the  control  routine  to  determine 
the  surface  of  the  pseudoelltpsold  for  the  block  2A 
analysis  of  the  head  end  with  web  (Section  5. 2. 1.2). 

SCT0R2  -  Subroutine  SCT0R2  Is  the  control  routine  to  deter¬ 
mine  the  surface  area  on  the  pseudoelllpsold  of  the 
projected  propeliant  core  In  the  block  2B  analysis 
of  the  head-end  with  web  (Section  5 .2. 1.3). 

SD1D13  -  Subroutine  SD 1 D 13  determines  the  center  cf  gravity 
and  moments  of  Inertia  of  the  straight  through  grain 
end  sections  (Section  5.3). 

SEGSUB  -  Subroutine  SsGSUS  Is  the  cylindrical  section  control 
routine  to  determine  the  mass  generation  rate,  port 
area,  perimeter  length*  end  cross-sectional  fuel 
area  for  each  increment  dividing  plane  and  mass  addi¬ 
tion  region  during  the  Internal  ballistic  solution 
(Section  5.2./,). 


SCT0R2  - 
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Subroutine 

SETPH  - 

SLOT 

SPLNiA  - 

SPLN2A  - 

SPLN3A  - 

SORT 

STUPPS  - 

STUPRS  - 

S2SX 

TDGRE  - 


Descriptions  (Continued) 


Subroutine  SETPH  Is  the  Internal  ballistic  solution 
control  routine  to  obtain  convergence  on  the  fore 
head  pressure  by  matching  the  grain  discharge  flow 
with  the  nozzle  flow  determined  from  the  nozzlts 
pressure.  The  performance  calculations  for  tnrust, 
total  Impulse,  thrust  coefficient,  etc.,  are  Inclu¬ 
ded  In  subroutine  SETPH  (Section  4.2). 

Subroutine  SLOT  determines  the  gas  properties  at 
the  discharge  section  of  e  slot  between  grain  seg¬ 
ments  for  the  non-steady  flow  solution  of  the  In¬ 
ternal  ballistics  (Section  4. 1.2.2) . 

Subroutine  SPLNIA  sets  up  subioutint  RSSPLN  to  per¬ 
form  the  initialization  for  the  spline  Interpola¬ 
tion  procedure. 

Subroutine  SPLN2A  sets  up  subroutine  RSSPLN  to 
determine  the  functional  values  and  derivatives 
for  the  arguments  of  the  spline  Interpolation 
functions  after  SPLNIA  has  been  executed. 

Subroutine  SPLN3A  sets  up  subroutine  RSSPLN  to 
determine  the  coefficients  of  the  Interpolating 
function  for  the  spline  Interpolation  procedure. 

Subroutine  SORT  fs  s  modification  of  the  IBSYS-13 
monitor  library  routine  to  determine  the  square 
root  of  negative  numbers.  Only  the  pasltlve  value 
of  the  argument  Is  transferred  to  the  designated 
register. 

Subroutine  STUPPS  stores  the  variable*  defining  the 
plane  A  produced  !n  sectors  jA  or  1L-*.  :ir  the 
block  1  analysis  In  subroutine  SCI  of  toe  head  end 
with  web  (Section  5. 2. 1.1. 2). 

Subroutine  STUPRS  stores  the  variables  defining 
plane  B  produced  in  sectors  3A  or  HA  for  the  block  1 
analysis  In  subroutine  SCI  of  the  head  end  with  web 
(Section  5. 2. 1.1. 2). 

Subroutine  S2SK  determines  the  sector  surface  area 
on  the  pseudoelllpsold  of  the  projected  grain  cross- 
section  for  the  block  23  analysis  lr«  subroutine 
SCT0R2  of  the  head  end  with  web  (Section  5.2. 1.3). 

Subroutine  TDUkE  determines  the  largest  real  root 
of  a  third  degree  polynomial  for  the  argument  X. 
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6.1.1 


Subroutine  Descriptions  (Continued) 


THETAR 

TISUB 

TRAN 

VFPPSB 

V0LSU6 

VSEC 

VSTRSB 

XRSUB 

XRSUBB 

XRTHR 


Subroutine  THETAR  determines  the  angle  between  the 
Z-axls  and  the  line  segment  RAT  for  the  block  1 
analysis  In  subroutine  SCI  of  the  head  end  with  web 
(Section  5- 2. 1.1.1). 

Subroutine  TISUB  determines  the  value  of  DELT  for 
the  steady  state  Internal  ballistic  solution  neg- 
lectlng  transient  effects  and  modifies  the  values 
of  thickness  burned  In  each  Increment  dividing 
plane  after  the  ballistic  solution  Ic  converged  for 
each  time  point. 

Subroutine  TRAN  transfers  the  geometry  constants 
from  the  permanent  common  storage  location  to  the 
working  array  common  storage  location. 

Subroutine  VFPPSB  determines  the  port  volume  of  eacu 
cylindrical  section  segment  and  sums  the  segment 
port  volumes  to  obtain  the  total  cylindrical  section 
port  volume. 

Subroutine  VOLSUB  Is  the  control  routine  which  de¬ 
termines  the  Initial  volume  for  the  block  3  analysis 
of  the  head-end  with  web  (Section  5.2. 1,4). 

Subroutine  VSEC  determines  the  volume  produced  by 
the  difference  of  volumes  of  two  oblate  spheroids, 
minus  the  volume  of  the  Igniter  hole  In  the  block  3 
analysis  of  the  head  end  with  web  (Section  5-2. 1.4). 

Subroutine  VSTRSB  determines  the  Initial  core  volume 
that  Is  present  In  the  head  end  with  web  for  the 
block  3  analysis  (Section  5-2, 1.4). 

Subroutine  XRSUB  determines  the  X-coordtnate  of  a 
general  point  on  the  perimeter  of  a  sector  In  the 
end  sections  (Section  5. 2. 3. 2). 

Subroutine  XRSUBB  determines  the  X-coordlnate  of  the 
RAT  line  segment  which  Is  a  radial  vector  from  the 
motor  axis  to  a  point  on  a  sector  perimeter  for  the 
analysis  of  the  head  end  with  web  (Section  5*2. 1.1.1). 

Subroutine  XRTHR  Is  a  set  up  subroutine  that  uses 
subroutine  XRSUBB  to  obtain  the  X-coordlnate  of  a 
point  located  on  the  perimeter  of  a  sector  In  the 
block  2B  analysis  of  the  head  en-1  with  web.  The 
angle  0^  between  tne  Z  axis  and  a  line  from  the 

motor  «k!s  to  a  general  point  In  a  sector  Is  also 
determined  (Section  5 *2. 1.3). 

_ _ 
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Subroutine  Descriptions  (Continued) 


YPSUB 


ZISUB 


-  Subroutine  YPSUB  determines  the  Y-coordlnate  of  the 
points  PO  and  P3  which  are  located  on  the  surface 
of  the  Inner  and  outer  ellipsoids  respectively  for 
the  block  2  analysis  of  the  head  end  with  web  (Sec* 
tlon  5*2. 1.2). 

«  Subroutine  ZISUB  determines  the  Z-coordJnate  pro¬ 
duced  by  the  Intersection  of  the  outer  ellipse  and 
the  line  normal  to  the  ellipse  at  PI  (Section 
5.2. 1.2). 


6.1.2 


Subroutine  Linkage  Table 


NAME  LOVER-LEVEL  CALLS  CALLED  BY 


AC0S 

SORT 

AESUB 

GAMA2A 

RASUBB 

AFP SUB 

LPT0 

SD ID  13 

AIGSUB 

MNCHN2 

THETAR 

ALRSU8 

PLNCKS 

XRTHR 

AWESUB 

P1SUB 

YPSUB 

ENDCSB 

P3SUB 

FDGRE 

RASUB 

AEPSUB 

AESUB 

ASESUB 

AESUB 

AC0S 

AEPSUB 

ALRSUB 

ARSSUB 

HESUB 

RASUB 

SQRT 

XRSUB 


AFPSUB 

AC0S 

LPDAPS 

SQRT 

A1BST 

AIBSU8 

AIBSUB 

SQRT 

MNCHN4 

AI8SUB 

A1BST 

AIBST 

SQRT 

MNCKN4 

SEGSUB 

AIGSUB 

AC0S 

SCT0R1 

SQRT 

ALRSUB 

AC0S 

AESUB 

SQRT 

ARSSUB 

SQRT 

AESUB 

ASESUB 

AEPSUB 

AWESUB 

ENDCSB 

RCSUB 

MNCHN2 

6.1.2 


Subroutine  Linkage  Table  (Continued) 


NAME 

ASTSUB 

ASUBC 

AWESUB 

BRAKSB 

BSUBC 

C0MPSC 

C0NV 

CSTRSB 

OPRASB 


LOWER  LEVEL  CALLS 


PT1AA 

GAKA2S 

THETAR 

GAHSUB 

TRAN 

P0SUB 

P1SUB 

P3SUB 

RASUBB 

XRSUBB 

AC0S 

HESUB 

SQRT 

SQ.RT 

GAMA2S 

THETAR 

GAMSUB 

TRAN 

P0SUB 

P1SUB 

P3SUB 

RASUBB 

XRSUBB 

CALLED  BY 

MNCHN3 


ASESUD 

SCI 

SCI 


LPT0 

MNCHN1 

MNCHN2 

MNCHN3 


SETPH 


SPLN1A  MNCHN4 

SPLN2A  RBSTSB 


SORT 


SCI 


AC0S 

ASESUB 

LPDAPS 

SQRT 


ENDCSB 


MNCHN1 


6.1.2 


Subroutine  Linkage  Table  (Continued) 


NAME  LOWER- LEVEL  CALLS  CALLED  BY 

FDGRE  A CDS  RCSUB 

SQRT 


5AMA2S  AC0S  ASUBC 

SQRT  BSUBC 

HAPSBC 

HASUBC 

HBPSBC 

H8SUBC 


GAMSUB  ASUBC 

BSUBC 

HABSBC 

HASUBC 

HBPSBC 

HBSUBC 


HAPSBC  GAMA2S  THETAR  SCI 

GAMSUB  TRAN 
P0SUB  XRSUBB 
PI  SUB 
P3SUB 
RASUBB 


HASUBC  GAMA2S  THETAR  SCI 

GAMSUB  TRAN 
P0SUB  XRSUBB 
PI  SUB 
P3SUB 
RASUBB 


HBPSBC  GAMA2S  THETAR  SCI 

GAMSUB  TRAN 
P0SUB  XRSUBB 
PI  SUB 
P3SUB 
RASUBB 


HBSUBC  GAHA2S  THETAR  SCI 

GAMSUB  TRAN 
P0SUB  XRSUB8 
PI  SUB 
P3SUB 
RASUBB 


Subroutine  Linkage  Table  (Continued) 


LOWER- LEVEL  CALLS 


NAME 

HONSUB 

HESUB 

I  MPT 

LBSU6 

LPDAPS 

LPTS 

ivMN 

hNCKNl 


SCI 

SCT0R1 

SCT0R2 

V0LSU8 


SQRT 


AFPSUB 

PT1AA 


AC0S 

C0KPSC 

SQRT 

TRAN 


INPT 

MNCHN1 

MNCHN2 

MNCHN3 

MNCHN4 


MNCHN3 


AESUB 

AWESUB 

PT1AA 

SD1D13 


MAIN 


SCT0R1 

S2jK 


ENDCSB 

MNCHN2 

MNCHN3 


SEGSUB 


C0MPSC  RBVSUS  MAiN 

ENDCSB 

PCHGiS 

PCHPLI 

PLNCMS 

PLNLCS 


CALLED 


Subroutine  Linkage  Table  (Continued) 


NAME 

LOWER-LEVEL  CALLS 

MNCHN2 

AC0S 

ASESUB 

C0MPSC 

LI'OAPS 

PT1AA 

RCSUB 

TRAN 

MAIN 

MNCHN3 

ASTSUB 

C0MPSC 

HONSUB 

LPDAPS 

RCSUB 

RG 1  SUB 

SQRT 

TRAN 

MAIN 

MNCHN4 

AIBST 

AIBSUB 

CSTP.SB 

0UTPUT 

RBSTSB 

RBSUB 

S".v“SUB 

SET  PH 

SQRT 

TRAN 

MAIN 

M0DTSB 

TI'^UB 

0UTPU7 

MSI  SUB 

SQRT 

SCT0R1 

MTISUB 

SQRT 

■  S2SK 

flUTPUT 

M0DTSB 

MNCNN4 

PCHGIS 

MNCHN 1 

PCHPL1 

MHCHN1 

PLNCNS 

AC0S 

SQRT 

MNCHN1 

165 


6.1.2 


Subroutine*  linkage  Table  (Continued) 


ffigsBsmgsm 


LOWER-LEVEL  CALLS 


CALLED  BY 


PLNLCS 

MNCHN1 

PSJSUB 

SQRT 

ASUBC 

BSUBC 

HAP SBC 

HASUEC 

HBPSBC 

HBSUBC 

3CT0R1 

PT1AA 

KESUB 

ASTSUB 

SD ID  13 

LPDAPS 

SQRT 

MNCHN2 

RGISUB 

P1SUB 

AC0S 

ASUBC 

SQRT 

BSUBC 
HAP SBC 
HASUSC 
HBPSBC 
HBSUBC 

P3SUB 

AC0S 

ASUBC 

SQRT 

BSIBC 

TP 'RE 

HAPSBC 

HASUBC 

HBPSBC 

HBSUBC 

RASUB 

AC0S 

AESUB 

SQRT 

SD ID  13 

RASUBB 

AC0S 

ASUBC 

SQRT 

BSUBC 

HAPSBC 

HASUBC 

HBPSBC 

HBSUBC 

SCT0R1 

RBSTSB 

CSTRSB 

;<8SUB 

66 


6.1.2 


Subroutine  Linkage  Table  (Continued) 


NAME 

LOWER-LEVEL  CALLS 

CALLED  BY 

R&SUB 

RBfTSB 

MNCHN4  ' 

SEGSUB 

SL0T 

RBVSU8 

MNCHN1 

RCSUB 

FDGRE 

ASESUB 

TDGRE 

MNCHN2 

MNCHN3 

RG 1  SUB 

PT1AA 

SdRT 

MNCHN3V  ’ 

PJ3E1SB 

SCT0R1 

S2SK 

R0PSB 

SCI 

RSSPLN 

SPLN1A 

SPLN2A 

SPLN3A 

SCI 

ASUBC 

HBSUBC 

HDNSUB 

BRAXSB 

R0PSB 

BSUBC 

SQRT 

DPRASB 

STUPPS 

HAPSBC 

srupRs 

KASUBC 

TRAN 

HBPSBC 

VSTRSB 

SCT0R1 

AlGSUB 

R0E1SB 

HDNSUB 

LBSUB 

SORT 

MSI  SUB 

YPSUB 

P0SUB 

RASUBB 

ZISUB 

! 


SCT0R2 


S2SK 


HDNSUB 


SuhynyttnA  Linkage  T»bl *  (Continued) 


NAME 

LOWER-LEVEL  CALLS 

CALLED  BY 

SD ID  13 

AC0S 

HESUB 

RASUB 

SQ.RT 

XRSUB 

PTiAA 

SEGSUB 

AIBSUB 

LPT0 

RBSUB 

SL0T 

SORT 

MNCKN4 

SETPH 

C0NV 

SORT 

VFPPSB 

MNCHN4 

SL0T 

RBSUB 

SORT 

SEGSUB 

SPLNIh 

RSSPLH 

CjTRSB 

SPLN2A 

SSSPLN 

CSTRSB 

SPLN3A 

RSSPLN 

(For  Diagnostic  Use) 

SORT 

AC0S  LPT0  SDID13 

AE3UB  MNCHN3  SEGSUB 

AFPSUB  MNCHN4  SETPH 

AIBST  MS  I SUB  SL0T 

AIBSUB  MTISUB  S23K 

AIGSUB  PLNCNS  TDGRI: 

ALRSU9  P0SUB  THETAR 

ARSSUB  PT1AA  V0LSUB 

AWSSUB  PI SUB  VSEC 

BRAKSB  P3SUB  YSTRSB 

DPRASB  RASUB  XRSUB 

ENJCSB  RASUBB  XRSCBB 

FOGPiE  RG I  SUB  ilRTHP 

GAMA2S  SCI  YPSUB 

HESU6  SCT0R] 


Subroutine  Linkage  Taole  (Continued) 


NAME 

LOWER-LEVEL  CALLS 

STUPPS 

SC! 

STUPRS 

set 

S2SK 

LBSU3 

MT 1  SUB 

R0E1SB 

SQRT 

XRTHR 

YPSUB 

2!  SUB 

SCT0R2 

TDGRE 

S(iRT 

P3SUB 

RCSUB 

Z  t  SUB 

THETAR 

AC0S 

ASU3C 

SQRT 

BSUBC 
HAP SBC 
HASUBC 
HBPSBC 
HBSUBC 

TISUB 

M03TSB 

TRAN 

ASUBC 

BSUBC 

HAPSBC 

HASUBC 

KBPSBC 

H8r.U8C 

VFPPSB 

SETPH 

V0LSVB 

SQRT 

HDNSUB 

VSEC 

VSEC 

SQRT 

XR‘<  HR 

V0LSUB 

YPSUB 

CALLED  BY 


LPT0 

MNCHN2 

MHCHN3 

MNCHN4 

SCI 


Subroutine  Linkage  fable  (continued) 


-  NAME 

LOWER-LEVEL  CALLS 

VSTRSB 

SORT 

set 

XRSUB 

SORT 

AESUB 

SO  1013 

XRSU5B 

SORT 

AS'JBC 

LJUSC 

HaPSBC 

FASUBC 

MBPSBC 

HBSUBC 

XRTHR 


XTHR 

AC0S 

S2SK 

SORT 

VSEC 

XRSUBB 

YPSUB 

AC0S 

SCT0R1 

SORT 

S2SK 

VSEC 

Z I  SUB 


TDGRE 


SCT0R1 

S2SK 


Flow  Charts 


*!>» I s  section  contains  macro-scoplc  flow  charts  In  schematic  form 
of  the  first,  second  and  third  level  subroutines.  The  Fortran 
listings  themselves  are  considered  to  be  the  flow  charts  for  the 
lower  level  subroutines.  In  order  to  facilitate  their  use  each 
routine  has  been  processed  by  a  special  program  which  re-numbers 
all  statement  numbers  In  steps  of  ten.  These  rout'nns  are  also 
annotated  by  comment  cards  which  define  the  computational  blocks 
and  locate  Important  logic  branches. 

Tne  following  schematic  flow  charts  are  Included: 


Name  Flow  Chart  No. 

MAIN  Program  1 

Subroutine  KNCHN1  2 

Subroutine  MNCHN2  3 

Compute  Plane  Constants  4 

Subrout  li:  AESllB  5 

Subroutine  MNCHN3  6 

Subroutine  SCI  7 

Subroutine  SCT0R1  8 

Subroutine  SCT0R2  9 

Subroutine  MNCHN4  10 

Subroutine  SEGSUB  li 

Sub  rout  iiie  SETPH  12 

Subroutine  TISUB  13 

Check  For  Case  Termination  14 


Frlmry  Uf!c  (ranch  ContnHW  hy  I  CM 


Flow  Chart  #1 
HA  IN  Program 


JnftUM**! 


Flow  Chart  #3 
Subroutine  KNCHN2 


Flow  Chart  #4 
Compute.  Plane  Constants 


Flow  Ch*rt  #6 
Sub rout  In*  MWCHN3 
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Incremental  Strip  t*  b*  Uaad 
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Calc.  Surface  Area  af  th* 
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Flow  Chart  #8 
Sub rout  In*  SCT0R1 


180 


^  Cntrl^ 

_ I 

tSet  U|t  Sedlua  erd  X,l  Coor- 
dlneto*  of  Current  Section 

| 

(*A 

|  Inltlalixe  Integration 

1  . 


call  mu»  > 


(  CAU  YPSU*  ~S 

i”  V, 


I  Correct  Integration  F»cc>~  (Li^) 
|  »nd  Incrunt  1  Coordinate  of  P| 


r«ot  ^ 


H Calculate  tncronenuF] 

I  Perlaetr.r  cf  Sector  I 

in: _ 

!  Calculate  *-«ord!.*ate»  and 
lAnjIei  of  Increaantul  Sector 


Calculate  Coordinated  of 
Increment*!  Stria  for  P 
of  Current  Sector  “ 


Calculate  Coordlrote*  of 
Increeentel  Strl(  for  P. 
of  Current  Sector 


<W*te  P# 
Ccardlnetee 


Calculate  i 


Calculate  Incremental  Area  ri  Peeude 
ClMpto  end  Velmne  of  Pre.olUnt 


/  lnte*^\ 
(ration  Cea- 


Flow  Chart  Hi 
Subroutine  SCTJ8R2 
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Vm4  Print  krn 


Flow  Chert  #11 
Subroutine  SEG'VUB  (2  of  2) 


Flow  Chert  #12 
Subroutine  SETPH 
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Flow  Clwrt  #13 
Subroutine  TISUB 
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6. 3,1.1 
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Stonge  Allocation 


The  program  Is  divided  Into  four  core  loads.  Subsection  6. “.1.1 
consists  of  a  table  showing  which  routines  *,e  required  In  each 
core  load.  Subsection  6. 3. 1.2  contains  a  diagram  explaining  the 
computer  allocation  of  the  program  code  and  geometry  tables.  Sub¬ 
section  6. 3. 1.3  contains  both  a’j  alphabetical  and  a  numerical  listing 
of  COMMON  allocation.  Subsection  6. 3. 1.4  contains  a  table  which 
defines  the  labeled  COMMON  allocation  of  the  geometry  tables. 
Subsection  6. 3. 1.5  consists  of  the  allocation  map  printed  out  by 
the  computer  system. 


Subroutine  Requirements  for  Each  Core  Load 
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X 
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X 
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SPLNIA 
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X 

MTISUB 

X 
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X 
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A 
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X 
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X 
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X 
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X 
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HBSUSC 
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X 

ZISUB 

X 

HDNSfJB 

X 

ROE 1  SB 

X 

HESU8 

X 

X 

X 

ROPSB 

X 

INPT 

X 

RSSPLN 

X 
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MFCHH4,  etc 


CUMMIN  Allocation 

The  following  pages  contain  output  from  a  special  Boeing  Com¬ 
pany  pre-compi lev  which  lists  numerically  and  alphabetically 
those  variables  which  are  assigned  to  unlabled  CUMMIN  through 
equivalence  statements. 

There  are,  in  addition,  two  block*  of  labeled  canon. 

1)  "TBJJK"  contains  the  geometry  tables  and  is 
explained  on  ptge  J$9» 

2)  "CRTIjfcr  contains  the  three  convergence 
criteria  CRP,  CRT  and  CRW  in  cells  1-3 
respectively. 


NUMERIC  listing  OF  PROGP»H  VARIABLES  IN  CO*  ARRAY 
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•i.4  Geometry  Table  Storage 

(Labeled  C0MM0N  Block  "TBLK") 
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End  Table 

Variable  Names 
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PM0IHD 
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Computer  System 


The  program  is  written  In  the  FORTRAN  IV  computer  language; 
howe-'er,  one  subroutine  (SQRT)  is  a  modified  FAP  version  of  the 
standard  library  routine  which  permits  a  negative  argument.  The 
program  is  designed  for  use  on  an  IBM  7094  computer  using  an  IBSYS, 
version  13>  monitor  system.  In  Oi  *er  to  conserve  core  space  the 
abridged  output  package  Is.  used  by  exercising  the  "ALTIO1  system 
option. 

Diagnostic  Tools 

Several  intermediate.,  diagnostic  print  statements  have  been 
permanently  built  into  the  program  to  aid  in  program  and  data- 
error  analysis.  The  complete  COMMON  regions  are  printed  before 
program  termination  in  the  event  of  any  convergence  failure.  In 
addition,  several  user-controlled  dumps  are  available  through  the 
use  of  Input  card  sets  48  and  49  as  explained  below. 

Card  Set  48  (Format  (12,  IX, 11E7.0)) 


Column 

Name 

Location 

Def in?  tion 

1-2 

N 

Card  Set  48 

4-10 

CKDUMP(l) 

MNCHN4 

Time  at  which  KDUMP(71)  will 
start  printing  dumps  I,  2  and 

3  (seconds) 

11-17 

CKDUMP(2) 

MNCHN4 

Time  at  which  KDUMP(71)  jI  1 1 
stop  printing  dumps  I,  2  snd 

3  (seconds) 

18-24 

CKDl’MP(3) 

MNCHN4 

Time  at  which  KDUMP(72)  will 
start  printing  dump  5  (seconds) 

25-31 

CKDUMP(4) 

MNCHN4 

Time  at  which  KDUMP(72)  will 
stop  printing  dump  5  (seconds) 

Card  Set  49  (Format 

(8011)) 

Column 

Name 

Location 

Definition 

1 

KDUMP(l) 

MNCHN4 

Fore*»head  and  Art-head  calcu¬ 
lations  of  T'UZ,  RBZ,  P,  AP,  ALP 
WDOT  and  DW  DOT 

*b 

KDUMP(2) 

SEGSUB 

Cylindrical  section  reference 
plane  calculations  of  TAUZ,  RBZ, 
P,  AP,  ALP,  W  DOT  and  DW  DOT 

3 

KDUMP(3) 

AIBST 

Start  transient  calculations  of 
PD  (!il),  DW  DOT,  WDOTD,  AMACHD 
In  cylindrical  section 

Diagnostic  Tools  {Continued} 


Column 

Homs 

Location 

Definition 

4 

kdump(4) 

SEGSUB 

Anisotropic  burn  perimeters  and 
burn  rates. 

5 

K0UMP(5 ) 

MNCHN4 

COMMON  and  TBLK  Data  regions. 

This  particular  dump  is  turned 
off  after  printing  unless  it  is 
being  controlled  by  KDUMP(72). 

6 

KDUMP(6) 

SEGSUB 

Fuel  area,  port  area  and  port 
perimeter  for  interpolation  of 
each  increment  dividing  plane 

71 

KDl'MP(71) 

MNCHN4 

SEGSUB 

AIBST 

This  flag  overrides  KDUMP(l), 
KDUMP(2)  and  KDUMP(3)  settings 
and  prints  these  dumps  according 
to  the  simulated  burning  times 
specified  in  CKDUMP(l)  and 
CKDUMP(2). 

72 

KDUMP(72) 

MNCHN4 

This  flag  overrides  KDUMP(5) 
and  prints  this  dump  according 
to  the  simulated  burning  t fires 
specified  in  CKDUMP(3)  and 
CKDUMP(4). 

i 

The  subroutine  SPLN3A  has  also  been  fnclu  ad  in  the  program  to 
permit  printing  out  the  coefficients  of  the  spline-fit  curves  used 
in  subroutine  CSTRSB,  If  those  coefficients  are  desired  the 
appropriate  coding  must  be  inserted  following  calls  to  SPLN1A  which 
perform  the  curve  fits. 


Nomenclature  I 

1 

The  following  pages  contain  the  nomenclature  list  which  contains 
ail  input  snd  output  variables,  all  variables  used  in  the  document 
text  and  a  significant  number  of  important,  internal  program 
variables.  In  order  to  facilitate  inclusion  of  additional  variables 
in  this  third  category  the  nomenclature  lists  have  been  punched  on 
data  cards.  A  copy  of  this  nomenclature  "deck!1  is  provided  with  the 
program  d*cks. 
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AfC  IN  INITIAL  LENGTH  OF  SECTOR  6*  DETERMINED  IN 

SUBROUTINE  PLNCNS  (FIGURE  5.24). 


IN  ANY  ZONE  AT  THICKNESS  TAU  MEASURED 
FROM  END  OF  SECTOR  NEAREST  MOTOR  AXIS 
TO  END  OF  SECTOR  OR  ZONE*  WHICH  EVER  IS 
SMALLER.  DETERMINED  IN  SUBROUTINE  AESUB 
(  FIGURE  5.21 > • 
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B72M  RADIANS  GEOMETRY  CONSTANT  DETERMINED  IN  SUBROUTINE 

PLNCNS  (  FIGURE  5.25) 
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XCGN  IN  AFT-HEAD  GEOMETRY  T AC i. E  CENTER  OF  GRAVITY 

LOCATION  FROM  AFT  TANGENT  PLANE  < DEPENDENT 
VARIABLE)*  INPUT  ON  SUBSET  CARD  NUMBER  40-6 


X-COORD I  NATE  OF  A  POINT  ON  THE  PERIMETER  OF 
A  SECTOR  FOR  THE  BLOCK  2B  ANALYSE  OF  THE 
HEAD-END  WITH  WEB  (FIGURE  5.16) 

DETERMINED  IN  SUBROUT  I N  XRSUBB. 
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-COORDINATE  CORRESPONDING  TO  XRX 
FIGURE  5.21) 
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YA5  IN  GEOMETRY  CONSTANT  WHICH  DEFINES  THE 

Y-COORDINATE  OF  THE  POINT  LOCATED  ALONG  LINE 
T 6  MAX  AT  THE  DISTANCE  95  F ROM  SIDE  t.C 
( F I  GORE  5.25)  DETERMINED  IN  SUBROUTINE 
PLNCNS. 
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RESULTS 


The  major  mod  I  float  ions  to  the  original  Thiokoi  Chemical 
Corporation  (TCC)  program  made  by  Boeing  were  done  during  the 
HiBEX  program.  HiBEX  stands  for  "High  6  Boost  Experiment." 

Because  of  its  extreme  flight  environment  and  a  new  stapled 
high  burning  rate  propellant,  it  was  not  known  if  the  internal 
ballistics  of  the  HiPEX  motor  would  be  affected,  i.e.,  ignition 
transient  interval,  maximum  chamber  pressure,  motor  burn  time, 
and  shape  of  the  chamber  pressure-time  trace.  Accordingly, 
program  modifications  were  made  as  discussed  in  Section  1.1. 

Resulting  predictions  using  these  modifications  were  compared 
with  measured  results.  The  dimensionless  fore-head  pressure¬ 
time  traces  are  shown  in  Figure  7*1*  The  prediction  indicates 
good  agreement  with  the  results  of  three  full  scale  motor  firings. 
During  the  HiBEX  program,  various  grain  configurations,  propellant 
formulations,  and  nozzle  throat  sizes  were  used.  The  three 
firings  for  which  data  are  shown  in  Figure  7*1  are  from  identical 
motor  configurations. 

Security  classification  of  the  HiBEX  program  prohibits  the 
discussion  of  specific  numerical  results  In  this  unclassified 
report.  Initial  analysis  that  as  corducted  is  reported  in 
Reference  4,  and  an  updated  analysis  using  an  anisotropic  burning 
rate  model  developed  from  small  scale  "Forty-Pound  Charge"  (FPC) 
motor  firings  to  predict  the  traces  shown  in  Figure  7*1  is 
discussed  In  Reference  5* 

Figure  7.2  shows  the  influence  of  internal  gas  flow  along  the 
prcpeilant  grain  and  the  propellant  burning  rate  model  on 
ballistic  predictions.  Curve  A  is  based  on  steady  internal 
flow  aid  isotropic  burning.  Curve  B  .'s  based  on  non-steady 
internal  flow  and  isotropic  burning.  Curve  C  represents  use  of 
the  complete  program  capability  except  for  an  accelerating 
reference  system.  *t  is  based  on  non-steady,  internal  flow  and 
anisotropic  burning  during  ignition  and  ta!loff.  The  method  of 
solution  for  these  and  other  program  options  are  discussed  in 
Section  8.1. 


